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Viscosity is an important fluid property which permits the 
development of shear stresses and resistance to flow. For static 
fluids, only normal stress or pressure exists. 

The ideal fluid was defined to be inviscid i.e. the viscosity has 
no effect. Thus, there were no frictional forces between moving 
fluid layers or between the fluid and bounding walls. 

The flow of a real fluid is more complex than that of an ideal fluid, 
owing to the phenomena caused by the existence of viscosity. 

Viscosity introduces resistance to motion by causing shear and 
friction forces between fluid particles and boundary walls. 

For flow to take place, work must be done against resistance 
forces, and in the process, energy is converted to heat. The 
inclusion of viscosity also allows the possibility of two physically 
different flow regimes. The effects of viscosity on the velocity 
profile also render invalid the assumption of uniform velocity 
distribution. Although the Euler equations may be altered to 
include shear stresses of a real fluid, the result is a set of partial 
differential equations to which no general solution is known. 

Laminar and Turbulent Flow 

In laminar flow, agitation of fluid particles is of a molecular 
nature only, i.e. at a length scale of the order the mean free path of 
the molecules. On the usual macroscopic scale of observation, fluid 
particles appear to be constrained to motion in essentially parallel 
paths by the action of viscosity, Fig. l.la. 


Fig. l.la 
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Fig. l.ib 


The shearing stress between adjacent layers in the laminar flow 

simple parallel flow) is completely defined by x = ju— , see Fig. 1.2 
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Fig. 1.2 

If the laminar flow is disturbed by wall roughness or some 
other obstacle, the disturbances are rapidly damped by viscous 
action, and downstream the flow is smooth again. Therefore a 
laminar flow is stable against such disturbances, but a turbulent 
flow is not. 

In turbulent flow, fluid particles don 't remain in layers, but 
move in a heterogeneous fashion through the flow, sliding past 
other particles and colliding with some in an entirely haphazard 
manner that results in a rapid and continuous macroscopic mixing of 
the flowing fluid, with length scales which are very greater than 
the molecular, scales in laminar flow, Fig. 1.1 b 

The random motion and the observed eddies in a turbulent flow 
suggest that both the inertia forces, associated with the 
accelerations during motion, and the viscous forces induced by the 
action of the viscosity, may be important. When the viscous forces 
are dominant, the flow might be expected to be laminar. When the 
inertia forces are dominant, the flow might well be turbulent. 

Both laminar and turbulent flows were demonstrated by Reynolds, 
(1883), with an apparatus similar to that of Fig. 1.4. 
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REYNOLD ’S EXPERIMENTS 

The effects of viscosity cause the flow of a real fluid to occur 
under two different conditions or regimes: that of laminar flow 
and that of turbulent flow. Reynolds was the first one 
demonstrated the characteristics of these regimes, with an 
apparatus similar to that of Fig. 1.1 a 




Fig. 1.4 (a) Reynolds Apparatus (b) different flow types 


Water flows from a tank through a bell-mouthed glass pipe, the 
flow being controlled by the valve. A thin tube leading from a 
reservoir of dye has its opening within the entrance of the glass 
pipe. 

Reynolds discovered that, for low velocities of flow in the glass 
pipe, a thin filament of dye issuing from the tube did not diffuse 
but was maintained intact throughout the pipe, forming a thin 
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straight line parallel to the axis of the pipe (Fig. 1.4b). As the 
valve was being opened, however, and greater velocities were 
attained, the dye filament wavered and broke, eventually diffusing 
through the flowing water in glass pipe. 

Since mixing of fluid particles during flow would cause 
diffusion of the dye filament, Reynolds deduced from his 
experiments that at low velocities this mixing was absent and that 
the fluid particles moved in parallel layers, or lamina, sliding past 
adjacent lamina but not mixing with them, this is the regime of 
laminar flow. Since at higher velocities the dye filament diffused 
through the pipe, it was apparent that mixing of fluid particles 
was occurring, or, in other words, the flow was turbulent. 
Laminar flow broke down into the turbulent flow at some critical 
velocity above that at which turbulent flow was restored to the 
laminar condition. 

Reynolds was able to generalize his conclusions from his dye 

stream experiments by the introduction of a dimensionless term 

Re, later called the Reynolds Number, which was defined by: 

n Vdp Vd 
Re = — — = — 
ju v 

In which V is the mean velocity of the fluid in pipe, d is the 
diameter of the pipe, and p, p and v are the specific mass, 
dynamic viscosity and Kinematic viscosity of the fluid flowing 
therein. 

The upper limit of laminar flow is defined by 
2100<Re cr <4000. The lower limit of turbulent flow, defined by the 
lower critical Reynolds number, is of greater engineering 
importance; it defines a condition below which all turbulence 
entering the flow from any source will eventually be damped out 
by viscosity. This lower critical Reynolds number thus sets a limit 
below which laminar flow will always occur; many experiments 
have indicated the lower critical Reynolds number to have a value 
of approximately 2100. Between Reynolds number 2100 and 4000 
a region of uncertainty exits. 

The concept of a critical Reynolds number to the flow of any 
fluid in cylindrical pipes ‘ one may predict that the flow will be 
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laminar if Re<2 100 and turbulent if Re>4000. However < critical 
Reynolds number is very much a function of boundary geometry. 

As conclusion, laminar flows (small Re) are characterized by 
low velocities, small length scales (e.g. small diameter pipes ) and 
fluids with high kinematic viscosity. Turbulent flows (large Re) are 
characterized by high velocities, large length scales, and fluids of 
low kinematic viscosity. 

1. Introduction to Laminar Flow 

The various ways of classifying flows have been met in earlier 
notes. Below is a summary of the broad classifications: 

Flow Description Steady/Unsteady the manner in which the fluid 
velocity varies with time. 

In steady flow the velocity is unchanging at any site with time. If 
you measure the fluid velocity at time T—0 and you measure the 
velocity some time T=tj after, at the same place, then the velocity 
would be the same. 

In unsteady flow, the velocities would not be the same. 
Uniform/Non-uniform the manner in which the fluid velocity 
varies in space. 

In uniform flow, the velocity is unchanging at any place with 
time. If you measure the velocity at a point, and then you measure 
the velocity at another point in the fluid, the velocities would be 
the same for any time at which the experiment is performed. 

In non-uniform flow, the velocities at the two points would be 
different. 

Laminar/T urbulent, the manners in which fluid particles move 
relative to each other. 

In laminar flow, individual particles of fluid follow paths that do 
not cross those of neighboring particles streamlines. Streamlines 
help in our understanding fluid flow. 

Recall from previous notes that a streamline is an imaginary 
curve drawn through a mass of fluid. The streamlines at time t 
form the family of curves everywhere tangent to the velocity field 
at that time (every point on the streamline is tangent to the net 
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velocity vector). When the flow is steady, the fluid particles move 
along the streamlines. There can be no net flow across a 
streamline. 


Flow Regimes 

At low flow rates, fluids move in "layers". The velocity 
components of the individual fluid elements do not cross the 
streamlines. This type of behavior is called laminar flow (or 
streamlines flow). In pressure driven flows laminar flow is 
predominant only at low flow rates. At higher flow rates, the 
streamlines are disrupted by eddies moving in all directions. 
These eddies make the flow turbulent, and although they cross and 
re-cross the streamlines, the net velocity is such that, the net flow 
from eddies is zero. 

Eddies form from contact of the fluid with a solid boundary or 
from two fluid layers moving at different speeds. As these eddies 
grow, the laminar flow becomes unstable and velocities and 
pressures in the flowing fluid no longer have constant or smoothly 
varying values. 

In pipes, relatively large rotational eddy form in regions of high 
shear near the pipe wall. These degenerate into smaller eddies as 
energy is dissipated by action of viscosity. The presence of eddies 
means that the local velocity is not the same as the bulk velocity 
and that there are components of velocity in all directions. 

The transition between laminar and turbulent flow is fuzzy. The 
behavior depends on entrance conditions and distance from the 
inlet. Often it is useful to speak of a transition region for flows 
that are neither laminar nor fully turbulent. 

Distinguishing Features of Laminar Flow 

-as the name suggests, the flow is characterized by layers of fluid 
traveling at different velocities. 

- Individual particles are held in place by molecular forces that 
prevent them from wandering outside their streamline. 

- Viscous forces predominate over inertia forces. 
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— Velocity gradients set up across the flow. 

Occurrences 

Associated with slow moving of viscous fluids; it is relatively rare 
in nature, but one example is flow of water through an aquifer. 
The velocities may be as low as a few meters per year. 

The text treats several examples that have specific interest to 
engineering. We will only be looking at a few of these. We look at 
flows in circular pipes since they are the most commonly used 
shape. 

Viscous Forces: Only Newtonian fluids are being considered in 
this course. Recall from the first few lectures that for such fluids, 
the shear stress developed between layers moving relative to each 

other was given by the following relation: x = pi— 

dy 

Where, x is the shear stress parallel to the fluid motion, du/dy 
is the velocity gradient perpendicular in the transverse direction 
and p is the dynamic viscosity. This relation is fundamental to our 
development of the relation. 

Pressure forces: Also of importance are the forces set up by 
pressure at the ends of the flow. Note that for steady, fully 
developed flow, there is no acceleration and so the total force is 
zero. 


Water at 13°C has a kinematic viscosity of 1.139 X 10" 6 m 2 /$ and flows in a cylindrical 
pipe of 30 mm diameter. Calculate the largest flowrate for which laminar flow can be 
expected. What is the equivalent flowrate for air? 
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Solution: 

Assume that the flow is one-dimensional and the flowrate 

Q = AV 

while the Reynolds number 

R = — 


(4.4) 


(7.1) 


For this calculation, the kinematic viscosity of water and the pipe diameter are given as 
v = 1.139 X 10“ 6 m 2 /s and d = 0.03 m. Taking R c = 2 100 as the conservative upper 
limit for laminar flow, 


2 100 = 


V(0.03) 


0.000 001 139 ’ 


V = 0.080 m/s 


Gwater = (0.080) X - (0.03) 2 = 5.64 X 10“ 5 m 3 /s © 


(7.1) 


(4.4) 


The water has the given viscosity at about 15°C so from Appendix 2 we have for air at 
zero altitude, which is at that temperature, 

v aiT = (/x/p)air = (1-789 X 10“ 5 Pa -s)/( 1.225 kg/m 3 ) = 1.46 X I0“ 5 m 2 /s 


and 


2 100 


(1.022) X 


V(0.03) 


0.000 014 6’ 


V = 1.022 m/s 


(7.1) 

(0.03) 2 = 7.22 X 10 ~ 4 m 3 /s = 13<2 wa[er • (4.4) 


2. Turbulent Flow and Eddy Viscosity: 

The characteristics of turbulent flow can be listed as follows: 

- Irregularity or randomness in time and space. 

-Diffusivity or rapid mixing. 

-High Reynolds number. 

-Three-dimensional vortices fluctuations. 

-Dissipation of the kinetic energy of the turbulence by viscous 

shear stresses. 

- Turbulence is a feature of fluid flow. 

The vorticity fluctuations symbolize the three-dimensional 
nature of the turbulence. The characteristic of dissipation 
makes it clear that there must be a continuous energy supply, 
from the mean flow to the turbulence, or it will decay. 

To analyze turbulence it is useful to focus on fluid particles. 
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These particles are observed to travel in randomly moving fluid 
masses of varying sizes called " eddies " these cause at any point in 
the flow, a rapid and irregular pulsation of velocity about a well 
defined mean value. This may be visualized as in Fig. 1.5, where V 
is the time mean velocity which can be measured directly with a 

t 

small pitot tube and V is the instantaneous velocity, such 
that:V —V', the over bar denotes a time average value. 

The instantaneous velocity V may be considered to be 
composed of vector sum of the mean velocity and the components 
of pulsations V x and V y , which can be measured by -hot-wire 
anemometer and its record is shown in Fig. 1.5, which 'discloses no 
regular period or amplitude because of the random nature of 
turbulence. The shown records allow the definition of certain 
turbulence characteristics: 



V X (or V y) 


'Jh i 


^rms ofy t (or v y ) 

v I'M 


"i: 7 " 


T* V/ W 


■Time 


J w 


Fig. 1.5 

-The root-mean square (r m s) values ) 2 or (y * ) 2 are a 

measure of the violence of turbulent fluctuations, i.e., the 
magnitude of departure of V' from V. 

-The relative intensity of turbulence (T x 2 ) 2 IV 

-The scale of the turbulence, the mean time interval between 
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reversals in the sign ofV x (or V y ), because it is a measure of the size 
of the turbulent eddies passing the point. 


Shearing stress in a simple parallel turbulent flow may be 
visualized by considering two adjacent points in a flow cross 
section, Fig. 1.6, at one of these points the mean velocity is V, at 
the other V + AV, where V & V + AV is taken to be the mean 
velocities of fluid layers. The turbulence velocity V y represents the 
observed transverse motion of small fluid masses between layers 
in both directions. The fluid masses move towards layer with 
small mean velocity is slowed down, while the fluid masses move 
towards layer with higher mean velocity is speeded up, i.e. fluid 
masses momentum is changed during their transfer process 
( tending to speed up the slower layer and slow down the faster 
one), just as if there were a shearing stress between layers. Thus the 
existence of shearing stress in turbulent flow is deducible from 
momentum considerations. 



The turbulent shear stress was given by Bossinesq (1877) who 
followed the pattern of the laminar flow equation and wrote: 


'T = S 


dZJ 


Where s is the eddy viscosity, it is property of the flow (not fluid 
alone) which depends primarily on the structure of the turbulence. 
Prandtl (1926) derived the following relation for the eddy viscosity in 
the form: 

s — pl 2 (dv / dy) 

Where he termed the distance l the mixing length, which is the 
mean distance which small aggregations of fluid particles are 
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transported by turbulence from region of one velocity to region of 
another. 

As near the boundary wall, the turbulence is strongly influenced by 
the wall (V y , V x — 0, fluctuation velocities of fluid particles due to 
turbulence, normal and along the direction of general mean 
motion), this let the mixing length l vary directly with distance from 
the wall y, where 


l = tcy 


Where K is von Karman constant, which has been determined by the use 
of experimental data, its nominal value is K= 0.4, then the turbulent 
shear stress yields: 


7- = ptdy 


2 . 



2 


Show that, if the velocity profile in laminar flow is parabolic, the shear stress profile must 
be a straight line. 


For sl parabolic relation between u and y , i> 

'T = 




y 

du 

dy 


From Eq 


Therefore, because du / dy J 2.a 1 y 7 nr = ^(2^^) and clearly, -/ oc y 


PROBLEM 7.3 


A turbulent flow of water occurs in a pipe of 2 m diameter. The velocity profile is measured 
experimentally and found to be closely approximated by the equation u = 10 + 0.8 In y, 
in which v is in metres per second and y (the distance from the pipe wall) is in metres. 
The shearing stress in the fluid at a point ^ m from the wall is calculated from measurements 
of pressure drop (see Section 7.10) to be 103 Pa. Calculate the eddy viscosity, mixing 
length, and turbulence constant at this point. 


The velocity gradient is needed and from the given profile 


v = 10 + 0.8 In y. 


dv _ 08 _ 0.8 

dy y (6.33) 


at a distance of 0.33 m from the wall. 

From Eq. 7.2, r = e ^ yields 103 Pa = 103 N/m 2 = s(2.4 s _1 ), 

dy 


s = 42.9 Pa-s (Note /x water «= 10 3 Pa-s) ® 
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dv 


From Eq. 7.3, r = - pu x v y = pi 2 ( — I yields (if p = 1 000 kg/m 3 ) 


From Eq. 7.6, r = 


dy 

103 N/m 2 = 1 000/ 2 (2.4) 2 , l = 0.134 m 
px 2 y 2 yields 

103 = 1 000k 2 ( | j (2.4) 2 , k = 0.401 


The magnitude of the eddy viscosity e when compared with the viscosity ji (approx- 
imately 0.001 Pa-s) is of special interest in that it provides a direct comparison between 
the (large) turbulent shear and (small) laminar shear for the same velocity gradient. The 
mixing length, /, when compared with the pipe radius is found to be about 10% of the 
latter dimension; this is a nominal value of correct order of magnitude, as is the turbulence 
constant, k. 


3. Fluid Flow past Solid Boundaries. 

For real fluid, experimental evidence shows that the velocity 
of the layer adjacent to the surface is zero (relative to the surface). 
This means that a velocity profile must show a zero value at the 
boundary. For laminar flow, surface roughness has no effect on the 
flow picture, Fig. 8.6, 




1 

J 
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t 
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/dv 




(a) Smooth boundary 




t 

y 

\ 

/do 




(b) Rough boundary 

Fig. 7.6 Laminar now over 
boundaries that act smooth. 

Fig. 1.7 

In turbulent flow, the presence of the boundary causes reduction of 
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mixing Length to zero value (i.e. the turbulence is completely 
extinguished ) and a film of viscous flow ( viscous sublayer) over the 
boundary results. A boundary surface is said to be smooth if its projection 
of protuberances are so completely submerged in the viscous 
sublayer and then they have no effect on the structure of 
turbulence, Fig.l.8.However, experiments have shown that 
roughness heights larger than about one-third of sublayer 
thickness will augment the turbulence and have some effect on the 
flow. 
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(a) Smooth boundary 


Turbulent / 
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1 — vr — — 

- • ;• f • , 

(b) Rough boundary 

Fig. 7.7 Turbulent flow over 
boundaries that act smooth. 

Fig. 1.8 

EXTERNAL FLOWS 

In external flow cases, energy or work is used typically to 
move the object through the fluid. In external flows, one is 
interested to determine the flow pattern ( including the pattern of 
viscous effect) around an object ( wing or flat plate, etc) immersed 
in the fluid, from which to determine the lift and drag ( resistance to 
motion) on the object. 

For internal flow, the focus is often not on lift and drag, but 
on energy or head losses, pressure drops and cavitation where 
energy is dissipated, because in internal flow, energy or work is 
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used to move fluid through passages. 

4. Characteristics of the Boundary Layer. 

Analysis made on the basis of ideal (inviscid) However 
streamlined bodies produces tow results: 

- The calculated drag on the body is negligible, while the observed 
drag is not. 

- The ideal fluid slips smoothly by the body, real fluid does not. 

For fluid with small viscosity (e.g. air and water), Prandtl 
suggested the boundary layer concept to explain the resistance of 
the streamlined bodies Fig. 8.8, flat plates parallel to the flow 
, airfoil and so forth. The essential point is that the frictional 
aspects of the flow are confined to the boundary layer and 
perhaps a wake behind the body, where the, flow is rotational, 
but outside the boundary layer, ( the viscosity of the fluid is 
essentially inoperative ), the flow is effectively frictionless and 
irrotational. 

Boundary layer phenomena are visualized on a smooth flat plate 
parallel to the oncoming flow, Fig T9. On this plate the boundary 
layer may be either laminar or turbulent. With smooth plate 
leading edge, a laminar boundary layer is to be expected adjacent 
to the upstream portion of the plate since the boundary layer is 
thin, viscous action intense, and turbulence is inhibited. As 
explained in section 8.1, the Reynolds number is used to define the 
character of boundary layer, and we can define two Reynolds 
number depend on the used length scale, Re x (.on the base of 
distance from the plate leading edge x), and Re d (on the base of 
boundary layer thickness S), in the form : 


R, 


V 


and 


R, 




v 
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Fig. 1.9 

Experiments show that the typical critical values for Re x =500 000 
and Res = 3 900, such that below these values, laminar boundary 
layer are to be expected. Whereas Reynolds numbers above these 
values the flow is transitioning to turbulent boundary layer with a 
thin viscous sub-layer. For excessive turbulence in the oncoming 
flow or a rough leading edge of the plate, it may imply a turbulent 
boundary layer beginning at the leading edge of the plate. 

The critical Reynolds number Re s is also very sensitive to pressure 
gradient as follow: 

-A favorable gradient, pressure decreasing in the direction of flow 
(flow in nozzle), stabilizes the laminar boundary layer and Re s > 

3900. 

-Conversely, an unfavorable pressure gradient, pressure increasing 
in the direction of the flow, as in diffuser, causes early breakdown of 
the laminar flow and Res <3900. 


[FLUID MECHANICS (1-B)] 


Page 15-CHI 


INTRODUCTION 


For real fluid flow an engineering approach for solution of such 
flow problem is to solve first the outer problem of ideal fluid 
motion about the body, ignoring viscous effect entirely. Then 
using the "outer” solution values of velocity and pressure at the 
surface of the body as approximate values for the edge of the 
boundary layer, the "inner" viscous flow problem, is solved. For 
streamlined shapes, this procedure gives, from the "outer" solution, 
the pressure distribution (including an accurate estimate of the lift 
force, if any) and, from the "inner" solution, an estimate of the 
friction force or drag on the shape is calculated. 

Now comparison between the streamlined body and flat plate 
boundary layers, Figs. 8.8 & 8.9 gives: 



Fig.1.10 

- Superficial similarity is noted immediately. 

- For the streamlined body: 

* The surface curvature may affect the boundary layer 
development due to either internal effects or induced separation. 

* The velocity just outside the boundary layer ( in the irrotational 
flow) change continuously along the body because of the 
disturbance of the overall flow due to body finite width. 

Flow Establishment- Boundary Layers 

It is considered here the simplest case, flow from a large 
reservoir into a cylindrical passage with a round entrance Fig. 
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3.23. The boundary layer will start at the entrance and will 
grow, with the distance along the tube ( accompanied by a 
diminishing core of irrotational fluid at the center of the 
passage. Thus the boundary layer steadily thickness until they 
meet and so envelop the whole flow, section-B, beyond which all 
velocity profiles will be identical. The flow is, .therefore established 
( developed ) and is everywhere rotational. Downstream of section 
B, the turbulence of wall friction is felt throughout the flow field. 

Of course, the flow in a boundary layer may be either laminar 

or turbulent. If the Reynolds number, Re = — , for the established 

v 

flow is less than 2100, it may be said that the established 
laminar flow has resulted from the growth of laminar 
boundary layers. In this case the zone of establishment has a 

x Rc 

length given by: — = — . For slightly higher Re, the flow may 

d 20 

be laminar up to and past 

^ = with subsequent transition to turbulent flow before the 
flow is truly established. 



In Re » 2100 the boundary layers are ultimately turbulent. For a 
well-shaped and smooth entrance, the boundary layers may be 
laminar at the upstream end of the zone of non-established flow, 
followed by turbulent boundary layers in the downstream portion of 
this region, Fig. 3.24. Between the turbulence (of boundary layers or 
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established flow) and solid boundary, there exists the viscous sub 
layer ( to satisfy the conditions of zero velocity on the walls). 



Note: Comparison of flat plate boundary layers with those of the 
pipe entrance 

Comparing Fig. 3.9 and 3.24, superficial similarity is noted 
immediately, however, for pipe entrance, these differences exist: 

The plate has been rolled into a cylinder, so it is not flat. 

The core velocity steadily increases, while for flat plate case, the 
free stream velocity remains constant. The pressure diminishment in 
the fluid in the downstream direction, where for flat plate there is no 
such pressure variation. 
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THEORY OF BOUNDARY LAYERS 


1-1. Introduction, 1-2.. Boundary layer definitions and characteristics of boundary layer disturbance 
thickness (5)-displacement thickness (5*)-momentum thickness (0)- energy thickness (a), 13-3 
Momentum equation for boundary layer by Vonkarman, 13-4.Laminar boundary layer, 13- 
5.turbulent boundary layer, 13-6.Total drag due to laminar and turbulent layers, 13-7 Boundary 
layer separation and its control - Highlights - Objective Type Questions-Theoretical Questions - Unsolved 
Examples. 


13-1. Introduction 

External flows past object encompass an extremely wide variety of fluid mechanics 
phenomena. Clearly the character of the flow field is a function of the shape of the 
body. For a given shaped object, the characteristics of the flow depend very strongly 
on various parameters such as size, orientation, speed, and fluid properties. 
According to dimensional analysis arguments, the character of the flow should 
depend on the various dimensionless parameters involved. For typical external flows 
the most important of these parameters are the Reynolds number Re —UL/v, where 
F- is characteristic dimension of the body. For many high -Reynolds -number flows 
the flow field may be divided into two regions. 

1. A viscous boundary layer adjacent to the surface of the vehicle 

2. The essentially inviscidflow outside the boundary layer 

W knows that fluids adhere the solid walls and they take the solid wall velocity. 
When the wall does not move also the velocity of fluid on the wall is zero. In region 
near the wall the velocity of fluid particles increases from a value of zero at the wall 
to the value that corresponds to the external ’’frictionless” flow outside the 
boundary layer (see figure). 
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The concept of boundary layer was first introduced by L. Prandtl in 1904 and since then it 
has been applied to several fluid flow problems. 





Ludwig Prandtl 
4 . 2.1875 15 . 8.1953 


When a real fluid (viscous fluid) flows past a stationary solid boundary, a layer of fluid 
which comes in contact with the boundary surface, adheres to it (on account of viscosity) 
and condition of no slip occurs (The no-slip condition implies that the velocity of fluid at 
a solid boundary must be same as that of boundary itself)- Thus the layer of fluid which 
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cannot slip away from the boundary surface undergoes retardation; this retarded layer 
further causes retardation for the adjacent layers of the fluid, thereby developing a 
small region in the immediate vicinity of the boundary surface in which the velocity of the. 
flowing fluid increases rapidly from zero at the boundary surface and approaches the 
velocity of main stream. The layer adjacent to the boundary is known as boundary layer. 
Boundary layer is formed whenever there is relative motion between the boundary and the 


fluid. Since r 0 = ju 


( du' 

aL 


, the fluid exerts a shear stress on the boundary and boundary 


exerts an equal and opposite force on fluid known as the shear resistance. 


According to boundary layer theory the extensive fluid medium around bodies moving 
in fluids can be divided into following two regions. 


(i) a thin layer adjoining the boundary called the boundary layer where the 
viscous shear takes place. 

(ii) A region outside the boundary layer where the flow behavior is quite likes that of 
an ideal fluid and the potential flow theory is applicable. 

1-2. Boundary Layer Definitions and Characteristics 

Consider the boundary layer formed on aflat plate kept parallel to flow of fluid of velocity 
U (Fig. 1-1) (Though the growth of a boundary layer depends upon the body shape, flow 
over a flat plate aligned in the direction of flow is considered, since most of the flow 
surfaces can be approximated to aflat plate and for simplicity). 

— The edge facing the direction of flow is called leading edge. 

— The rear edge is called the trailing edge. 

— Near the leading edge of a flat plate, the boundary layer is wholly laminar. For a 
laminar boundary layer the velocity distribution is parabolic. 
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Turbulent 

layer 

Overlap layer 
Buffer layer 
Viscous sublayer 


4 Laminar boundary 

layer 


. Transition 
region 


Re=5. 10' 


,5 


Turbulent boundary 

layer 

Re=3 10 6 y 


Fig.,1 Boundary layer on aflat plate. 

— The thickness of the boundary layer (S) increases with distance from the leading 
edge x, as more and more fluid is slowed down by the viscous boundary, becomes 
unstable and breaks into turbulent boundary layer over a transition region. 

For a turbulent boundary layer, if the boundary is smooth, the roughness 
projections are covered by a very thin layer which remains laminar, called laminar 
sub layer. The velocity distribution in the turbulent boundary layer is given by Log 
law or Prandtl's one-seventh power law. 

The characteristics of a boundary layer may be summarized as follows: 

(i) 3 (thickness of boundary layer ) increases as distance from leading edge x 
increases. 

(ii) S decreases as U increases. 

(ill) 3 increases as kinematic viscosity (v) increases. 

( 

(iv) T t) = p — ; hence r Q decreases as x increases. However, when boundary 

\3J 

layer becomes turbulent, it shows a sudden increase and then decreases with 
increasing x. 

(v) When U increases in the downward direction, boundary layer growth is 
reduced. 
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( vi ) When U decreases in the downward direction, flow near the boundary is 
further retarded, boundary layer growth is faster and is susceptible to 

separation. 

(vii) The various characteristics of the boundary layer on flat plate (e.g. 
variation of 3, r 0 or force F) are governed by inertial and viscous forces; hence 

they are functions of either — or — . 

v v 

Ux 

(viii) If - — <5xl0 5 ... boundary layer is laminar ( velocity distribution is 
v 

parabolic). 

Ux 

If — >5xl0 5 ... boundary layer is turbulent on that portion (velocity 
v 

distribution follows Log law or a power law). 

Ux 

( ix) Critical value of — at which boundary layer change from laminar to 

v 

turbulent depends on: 

- Turbulence in ambient flow, 

- Surface roughness. 

- Pressure gradient. 

- Plate curvature, and 

- Temperature difference between fluid and boundary. 

(x) Though the velocity distribution would be a parabolic curve in the laminar sub-layer 
zone, but in view of the very small thickness we can reasonably assume that velocity 
distribution is linear and so the velocity gradient can be considered constant. 

1 3-2- 1. Boundary layer thklmess ( S) 

The velocity within the boundary layer increases from zero at the boundary surface to the 
velocity of the main stream asymptotically. Therefore the thickness of the boundary layer is 
arbitrarily defined as that distance from the boundary in which the velocity reaches 99 per 
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cent of the velocity of the free stream (u = Q-99U). It is denoted by the symbol S. This 
definition however gives an approximate value of the boundary layer thickness and 
hence S is generally termed as nominal thickness of the boundary layer. 

The boundary layer thickness for greater accuracy is defined in terms of certain 
mathematical expressions which are the measure of the boundary layer on the flow. 
The commonly adopted definitions of the boundary layer thickness are: 

1. Displacement thickness { 3 *) 

2. Momentum thickness (6) 

3. Energy thickness ( S e ). 

13-2-2. Displacement thickness ( S*) 

The displacement thickness can be defined as follows: 

"It is the distance, measured perpendicular to the boundary, by which the 
main/free stream is displaced on account of formation of boundary layer. " 

Or "It is an additional "wall thickness" that would have to be added to compensate for 
the reduction inflow rate on account of boundary lay er formation” . 

LA I LA Page 6 5/9/2015 

Figure 2: Definition of boundary layer thickness: 

(a) Standard boundary layer (u = 99%U), (b) Boundary layer displacement thickness. 

Far-away free stream “sees ” solid wall as if it was being displaced up by d* (x) at 


x. 



The displacement thickness is denoted^ by S*. 

Let fluid of density p flow past a stationary plate with velocity U as shown in the 
Fig.2. Consider an elementary strip of thickness dy at a distance yfrom the plate. 
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Assuming unit width, the mass flow per second through the elementary strip. = pu.dy 
Mass flow per second through the elementary strip ( unit width) if the plate were not 
there = pUdy Reduction of mass flow rate through the elementary strip = p(U-u)dy 
[The difference (U-u) is called velocity of defect]. 

s 

Total reduction of mass flow rate due to introduction of plate - J p(U -u)dy 

o 

(if the fluid is incompressible ). Let the plate is displaced by a distance 8* and velocity 
of flow for the distance <5* is equal to the main/free stream velocity (i.e. U). Then, loss 
of the mass of the fluid/sec, flowing through the distance 3* 

= pUS* 

Equating Eqn. (ill) and (IV), we get 

s 

pU 8* = J (U - u)dy 
o 



1 - 2 - 3 . Momentum thickness ( 9) 

Momentum thickness is defined as the distance through which the total loss of 
momentum per second be equal to if it were passing a stationary plate. It is denoted 
by 6. 

It may also be defined as the distance, measured perpendicular to the boundary of 
the solid body, by which the boundary should be displaced to compensate for 
reduction in momentum of the flowing fluid on account of boundary layer formation. 

Refer fig. 1.2. Mass of flow per second through the elementary strip = pudy 

Momentum/sec, of this fluid inside the boundary layer = pudyxu = pu dy 

Momentum/sec, of the same mass of fluid before entering boundary layer - puUdy 

Loss of momentum/sec. = puUdy — pu 2 dy = pu(U — u)dy 
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Total loss of momentum/sec. 



Let 6 = distance by which plate is displaced when the fluid is flowing with a 
constant velocity U. 

Then loss of momentum/sec, of fluid flowing through distance 6 with a velocity U 


All three boundary layer thickness definition 8, Si, 0 are use in boundary layer analysis. 

Viscous/friction effect in a boundary layer causes mass and momentum flux deficits (at 
cross section x ) relative to ideal inviscid/frictionless flow. If it were not for the viscous effect, 
the velocity profile would be full/uniform. 

In general, the mass and momentum flux deficits increase with x and the boundary layer grow 
with x. 

1. Momentum flux deficit increases because the fluid in the boundary layer is continuously 
subjected to opposing (net) shear as the BL develops along the wall. Here, we focus mainly in 
the case of no pressure gradient (uniform pressure). 

- If a fluid element within a BL is considered, the net shear is illustrated in Fig. 2. 

- if a semi-infinite CV is considered, the net shear is the wall shear stress as illustrated in Figs. 

2 and 8. 

2. Recall that the concept of flux is associated with a surface, i.e., flux through a surface. Here 
the surface is a semi-infinite plane (0 < y ->°°) at x.] 

Motivation: Due to the force/mechanical interaction between the body and the fluid stream 
and Newton’s 3rd law of motion, the lost of momentum of fluid stream must relate to force on 
the body. 

Question: In order to study the BL a little more conveniently, how should we parameterize 
(/describe) the deficit in mass and momentum fluxes at any cross section x? 

Viscous/friction effect in a boundary layer causes fluid in the boundary layer to be retarded, 
losing its velocity/momentum, and resulting in mass flux deficit relative to ideal 
inviscid/frictionless flow. If it were not for the viscous effect, the velocity profile would be 
full/uniform. 

Momentum flux deficit of the actual mass flux: dM& x,deficit = (U - u)(pudA). [Note that the 
local free stream velocity is used as a reference/datum.] 


r\ 

— pOU , Equating Eqns. (i) and (ii), we have 
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Viscous/friction effect in a boundary layer causes fluid in the boundary layer to be retarded, 
losing its velocity/momentum, and resulting in momentum flux deficit relative to ideal 
in viscid/frictionless flow. If it were not for the viscous effect, the velocity profile would be 
full/uniform. 

Note that we account for the momentum flux deficit (U - u ) of the actual mass flux pud A , 

not of an ideal mass flux p UdA . [We want to know the lost of momentum of the actual 
mass.] 

Comments on Thicknesses: Mathematical Definitions and Physical Meanings 

In “deriving/arriving at” the “definitions”! of the displacement and momentum thicknesses 
given by Eqn. (1) and (2), respectively, we have (implicitly or explicitly) restricted 
ourselves to, say, incompressible flow, in order to gain some physical meanings of the 
thicknesses for such flow, namely, Definitions Physical Meanings Displacement thickness 

- related to - Mass flux deficit. 

Momentum thickness - related to - Momentum flux deficit. 

The question may naturally arise whether these thicknesses “can be used” for, say, 
compressible flow. In this regard, there are couples of points to consider: 

1. If we want to preserve the definitions: 

These definitions Eqns. (1) and (2), can simply be taken as definitions, hence, the 
“arriving/deriving” that was done is irrelevant and not necessary - one need not do it in 

order to so define “5* ” and “0 .” [However, one would probably and initially then not 

know what it was for or represented.] If so desired, both “5*” and “0” can be calculated for 
compressible flows. 

[At any instant in time t and at any cross section x, one can evaluate the integrals regardless 
of whether it is incompressible or compressible flow.] 

The questions are then rather: 

0 Whether they preserve the physical meanings (still have the same physical meanings) as 
in the case of incompressible flow, i.e., 5* corresponds to mass flux deficit and 0 
corresponds to momentum flux deficits? 

0 If they do not, what physics do they represent/mean in compressible flow? 

2. If we want to preserve the physics/physical meanings: 

On the other hand, if we want to preserve the physical meanings (if the two previously 
defined definition do not preserve the physical meanings), we have to redefine “5*” and 

“0” accordingly. New definitions - preferably with wider scope - must be sought to 
represent the mass and momentum flux deficits in compressible flow (and preferably still 
applicable also to incompressible flow). 

Question: 

Preserve its physical meaning in compressible flow? That is, does still represent mass flux 
deficit. If not: 

- What physics does it represent? And - How should we redefine 5* in order to preserve 
the physical meaning of mass flux deficit? What about 
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1 - 2-4 Energy thickness (SJ 

Energy thickness is defined as the distance, measured perpendicular to the 
boundary of the solid body, by which the boundary should be displaced to compensate 
for the reduction in K.E. of the flowing fluid on account of boundary layer formation. 
It is denoted by S e . 

Refer fig. 1.2 

Mass of flow per second through the elementary strip - pudy 

K.E. of this fluid inside the boundary layer 

1 2 1 

-—mu — 

2 2 


(pudy) 


u 


K.E. of the same mass of fluid before entering the boundary layer 
= ^ (pudy)U 2 

Loss of K.E. through elementary strip 


- ( pudy)U 2 ~ ~ (pudy)u 2 ~ ~ pu{u 2 — u 2 ^d 


:. Total loss of K.E. of fluid J —pu(u 2 — u 2 ]c/y 


Let d e = distance by which the plate is displaced to compensate for the reduction in 
K.E. Then loss of K.E. through d e of fluid flowing with velocity U ^ ( pUS e )U 2 


Equating eqns. (i) and (ii). We have 
\{pus e p 2 =\\pu[u 2 -u 2 )dy 

Z 0 2 
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Example LI . 


The velocity distribution in the boundary layer is given by: — = — , where u is 

U 8 

the velocity at a distance y from the plate and u = U at y =S, 3 being boundary 
layer thickness. Find : 

( i) The displacement thickness, ( ii) the momentum thickness, 

8* 


(iii) The energy thickness, and 


(iv) the value of 


9 


Sol. Velocity distribution: — = — 

US 


(Given) 


8 / \ 

U ' 


( i) The displacement thickness, 8 : 8* = j* |^1 — — dy 


S' = \l-l\ly v“=> 

J I sP jjs 


-iS 


y 


28 


8 = 


:2\ 


8- 


28 


=S--=- 

2 2 


(Ans.) 


(ii) The momentum thickness, 6: 


e= \^\i-F\ dy 


o 

s 




A 


y 


5f 2\ 

y y 


8 8< 


dy 
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f 2 


6 == 


y y 


.3 A 


25 35 2 


J 0 


” 2<5 3 < 5 2 ” 6 (A “' ) 


(iii) The energy thickness, S e : 




u 


2 A 


[/ 


dy 


*•=! 


8 f 

y_ 

5 


y 


y_ 

s : 


2\ 


dy = ) 


y y 


.3 A 


y 


25 45' 

, 8 . = S - 


nS 


JO 


ov^ S; 


5 2 5^_5 

25 45 3 ~ 4 


dy 


( iv) The value of — : 

0 


d* 

^ = ^7 = 3 - 0 (A™-) 
" % 


Example.] -2. 


The velocity distribution in the boundary layer is 

— - — , S being boundary layer thickness. 

JJ 2 5 2 S 2 


Calculate the following: 


5* 


(i) The ratio of displacement thickness to boundary layer thickness — . 

5 


(ii) The ratio of momentum thickness to boundary layer thickness 


0 


given by 
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Sol. Velocity distribution: - = — — 2— 

U 2 8 2 S 2 


(given) 


(i) 

8 


8 = 1 - 


u 

U 


dy = )\l- 


3 y ly 
+ 


2 A 


2 8 2 8 : 


dy 


y 


3 y + i r 


2 2 8 2 3 8 : 


Jo 


_ 3 8 2 18 

8 + 


4 8 


2 3 8 2 

8 ^ = 

8 


3 8' 

8--8 + 


V 


12 


4 6 

(Ans.) 


58 

12 


(h'J 


<9 


e = \-\i 

Uj 


u 


3 


dy 


M 


& f r\ 2 '\f 

3 y ly 


o v 


2 <5 2 <5 2 


i 2 i + *y 


2 A 


yv 


2 2 8 2 


dy 


H 


S f r\ r\ 2, q 3 't 2 q 3 -i 4 A 

3 t 9y 3 y 1 y 3 y 1 y . 

+ _2J_ + — ^ 


o v 


2 8 4 8 2 4 8 3 2 8 2 4 8 3 4 8 Z 


J 


s f 


H 


ov 


3 Z 
2 8 


f 9y 2 1 ^ 

+ -Ar + 


4 (J 2 2£‘ 


y 


3 y_ + 3 y_ 
4£ 3 4£ 3 


33 1/3 


4 / 




m 


S f r\ -I -I 2 Q 3 1 4 A 

3y 11 y + 3y r __ly_^ 


OV 


2 £ 4 2 S 3 4 <T 


J 
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0 = 


3 y 11 y s 

— x x — - 

2 28 A 3 S 2 


3 y 1 y 
+ -x^—~ — x — — — 

2 4 8 3 4 5 S 4 


-\8 


Jo 


_ 3 8 1 11 8 5 3 8 4 

6 — — x x — -H — x — - 

2 25 4 3 S 2 2 4 8 3 


1 8 

x 7 

4 58 4 


nS 


Jo 


d =—8 

120 


6>_ 19 
5 ~ 120 


Example. 1 .6. 

Explain what you understand by boundary layer thickness and displacement 
thickness. 

Determine the relationship between the two for a boundary layer which is (i) 
laminar throughout and (U) turbulent throughout. 

Assume (I) in the laminar boundary layer the flow obeys the law, shear stress 

r — p — , where p the viscosity, which leads to the velocity profile 

dy 


(fJ —u) = k(<5 — y) 2 , where U is the free stream velocity, u is the velocity at a 
distance y above the plate and k is a constant; 


u 


(2) The velocity distribution in the turbulent boundary layer is given — = 




Sol. (a) Boundary layer thickness: It is defined as the distance from the boundary of 
solid body measured in Y-direction to the point, where the velocity of fluid is 
approximately equal to 0.99 times the free stream (U) velocity of fluid. 

Displacement thickness: It is defined as the distance, measured perpendicular to 
the boundary of the solid body by which the boundary should be displaced to 
compensate for the reduction inflow rate on account of boundary layer formation. 
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(i) When the flow is laminar throughout: 
Velocity profile: (U —u) = k(S - y) 2 
where, 

U = free stream velocity, 
u - velocity at a distance y above plate, and 
k = constant. 

Let <5* = displacement thickness, and 

S - boundary layer thickness. 

sr \ 




iy 


Dividing Eqn. (A) by U, we get 


l- U - = ^(S-yf 
U U 


(A) 


yfdy = 

k (S-y) 3 

A 

8 

k(5 - y) 3 


U 3 

0 

3 U 


When y = 0, u = 0 ie. on the surface of the plate. 


kS 3 

3U 


Substituting these parameters in Eqn. (A), we get U = kS‘ 


Substituting the value of U from (C) in (B), we get 8 — 


kS 3 8 


3 k8‘ 


= - (Ans.) 


(B) 

...(C) 


(ii) When the flow is turbulent throughout: 

i 

Z N 


u 


r 


Velocity profile: — = 
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7 y 

y — x — 

8 


-is 


1 

5~ 7 


8 


Jo 


(Ans) 


Example. 1.8. 

The velocity distribution in the boundary layer over the face of a spillway was observed to 


be: 


u 

U 


s \ 0.22 

' yi 

kSj 


The free stream velocity U is 20 m/s and boundary layer thickness 5 cm at a certain 
section. The discharge is 5 m /s per meter length of spillway. Calculate displacement 
thickness, energy thickness and loss of energy up to section under consideration. 


u 


Sol. Given: Velocity distribution: — = 

U 


/ \ 0.22 

fy\ 


\° J 


J 

U = 20 m/s; 6 = 5 cm; Q = m/s per meter length of spillway, 
S\S e , E l : 

The displacement thickness is given by the equation: 

~\S 


■>wBH 

o v ^ J o 


1- 


/ X 0.22 

( y\ 


dy = 


y 


0) 


1 .,W 


0.22+1 


0.22 


1.22 


Jo 


= S 


1.22 1.22 


0.22 0.22 

6 = — — x5 = 0.9016 cm 


1.22 


(Ans.) 


The energy thickness is given by the equation: 
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s e = t- 
l U 


s '=\ 




2 

5 


f \ 

3 

s 


0.22 


0.66 

1- 

U 1 


=i 

VI 

u ) 


dy = \ 

y- 

— 

y 



Jj) 


J 

0 

u 

Jjj 


y J 

0 

|UJ 


UJ 



dy 


. 7 ) 


i x (y) 


0.22+1 


0.22 


1.22 


7) 


1 .+) 


0 . 66+1 




0.66 


1.66 


S e =- 


8 

<5 

1 

1 

= 5 

"1.66-1.22“ 

1.22 

1.66 

1.22 

1.66 

_ 1.22x1.66 _ 


= 1.086 cm 


(Ans.) 


The energy loss per m length of spillway is: 

E L =^(pxS t xU)xU 2 = \pS e U 2 

1 ( 

E T - — xlOOOx 


Ix(20) 3 xl0 -3 kNm/s = 43.44 kNm/s 
v 100 ) 


Energy loss in terms of m of head: 


43.44x1000 


pgQ 1000x5x9081 


= 0.8856m 


(Ans.) 


Example. 1-9. for steady Poiseuille flow in a pipe of radius R, obtain an expression for 
ratio of the displacement thickness (75*) to momentum thickness (6). 


Pipe 
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Sol. Radius of the pipe = R 

For steady Poiseuille flow in a circular pipe, the velocity distribution is given by, 


u = -J-M(R 2 -y 2 ). 
4ju dx 


Where y being measured from the center of the pipe. 
At y=0, u=U (maximum Velocity) U — 


•(0 




A/u dx 


m 


From the definition of displacement thickness, we have 
nR 2 -tt{r-5*) 2 \j = \2tv(U -u)dy 


R 


or 


2jfcS*U = \27ty(U -u)dy 


[Neglecting the term containing (S*) , since d* is very small] 
Rr 1 R v 2 R 

4 


or 


i R F \ \ R 2 

<$* = — f ( 1 - — ydy = — [ — ydy = 

j?J| tt Y r[r 2 


R\ 'A U J 


Dividing (i) by (ii), we have, — - 


u R 2 - y 2 1 y‘ 


= 1 


or 


U R 




V"' 

U ) 


y 


R‘ 


From the definition of momentum thickness, we have 


\tzR 2 — 7t(R — 6) 2 ]p U = J 2p n y.dy.u(u — u) 


On simplification, we get 


2ttR6U 2 = J 27tyu(U - u)dy 
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23 - 3 . Momentum Equation for Boundary Layer by Vonkdrman 

Vonkdrman suggested a method based on the momentum equation by the use of 
which the growth of a boundary layer along aflat plate, the wall shear stress and the 
drag force could be determined (when the velocity distribution in the boundary layer is 
known). Starting from the beginning of the plate, the method can be used for both 
laminar and turbulent boundary layers. 

Fig. 1-4. shows a fluid flowing over a thin plate (placed at zero incidence ) with a 
free stream velocity equal to U. Consider a small length dx of the plate at a distance x 
from the leading edge as shown in Fig. 13-4 (a); the enlarged view of the small length 
of the plate is 

Mass diagram: 
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X-Momentum diagram: 

In this diagram, we shall include the non-uniform pressure force. However, in 
the derivation below we shall consider only the case of uniform pressure. For the 
case of non-uniform pressure, we shall postpone it to later section. 



Fig. 1-4. Mass and momentum diagrams for the derivation of Vonkarman's momentum integral 
equation 

Consider unit width of plate perpendicular to the direction of flow. 

Let ABCD be a small element of a boundary layer ( the edge BC represents the outer 
edge of the boundary layer). 

Mass rate of fluid entering through AB 

5 

0 

Mass rate of fluid leaving through CD 

d 

+ — 

dx 



v 

jpudy 


dx 
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Mass rate of fluid entering the control volume through the surface AB = mass 
rate of fluid through BC - mass rate of fluid through CD 


= J pudy 


+ ■ 


d 

dx 


J pudy 


dx — J pudy = 


d 

dx 


J pudy 


dx 


The fluid is leaving through BC with a uniform velocity U. 

Momentum rate of fluid entering the control volume in X-direction through AB 


= \pu 2 dy 
o 

Momentum rate of fluid leaving the control volume in X-direction through CD 

d 


= J pu 2 dy 


+ ■ 


dx 


\pu 2 dy 


dx 


Momentum rate of fluid entering the control volume through BC in X-direction 


d_ 

dx 


8 


\ pudy 


dxxU = 


d_ 

dx 


8 


J puUdy 


dx 


Rate of change of momentum of control volume = momentum rate of fluid through 
AB - momentum rate of fluid through BC - momentum rate of fluid through CD 


= \pu 2 dy 


+ ■ 


d 

dx 


\pu 2 dy 


dx- J pu 2 dy 


d 


dx 


J puUdy 


dx 


d_ 

dx 

d 


8 


8 


J pu 2 dy-^ puUdy 


dx 

d 


J (pu 2 dy - puUdy ) 


dx 


dx 


dx 


8 


p\(u 2 -uu)dy 

o 


dx 


(p is cons tan tfor incompressible fluid ) 

d 


= P 


dx 


~s f 

K“ 

2 -uU ) 

_o 



ly 


dx 
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As per momentum principle the rate of change of momentum on the control volume 
ABCD must be equal to the total force on the control volume in the same direction. 
The only external force acting on the control volume is the shear force acting on the 
side AD in the direction D to A (fig. 1-4 b). The value of this force (drag force) is 
given by. 

A F d - z Q xdx 

Thus the total external force in the direction of rate of change of momentum: 

= -r 0 xdx 

Equating the equations (1-4) and (1-5), one has 







But 



— cly — momentum hickness 

U) 



pU 2 dx 
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Eqn. (1.7) is known as: 

Von Karman momentum equation for boundary layer flow, and is used to find out 
the frictional drag on smooth flat plate for both laminar and turbulent boundary 
layer. 

The following boundary conditions must be satisfied for any assumed velocity 
distribution: 


(i) T the surface of the plate: y — 0, u = 0, 


(ii) At the outer edge of the boundary layer: y - 

x dii 

y = S, — = 0 
dy 


— = finite value 
dy 

S, u = U, 


The shear stress for a given velocity profile in laminar, transition or turbulent zone 
is obtained from equation (1-6) or (1-7). Then drag force on a small distance dx of a 
plate is given by, 

AF d = shear stress x area = r c x(B x dx) [assuming width of the plate as unity], 

where B - width of the plate) 

Total drag on the plate of length L one side, 

L 

A F d = J A F d = | r 0 xBxdx 
o 


The ratio of the shear stress r a 



, is known as "Local coefficient of drag" (or 


coefficient of skin friction) and is denoted by C D 


i.e. 


r* - T o 

1 


pU‘ 


1 


- The ratio of the drag force to the quantity — pAU is called "Average coefficient of 
drag" and is denoted by C D 
i.e. C n = F ° 


D 


1 pAU 2 
2 


Where, p = mass density of fluid, 

A - area of surface/ plate, and 
U - free stream velocity 


NOTE: It is interesting to note here that, from Eqn. (7), the increase in drag force 
dD is directly proportional to the increase in the momentum thickness dO. That is, 
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we can then summarize the von Karman’s momentum integral equation for a 
boundary layer in uniform pressure as follows: 


von Karman’s Momentum Integral Equation: 


Note: la fluid, the dmmmlm form is often preferred owing partly to the concept of measuring things with its 
own scales which often gives insight into the problem as well as the relative order of magnitude of various 


Certainly, it can be a little cumbersome in some situations, e.g., calculate the magnitude of t v in Pa, etc. 

As we can see, the momentum thickness so defined from the momentum flux 
deficit relates directly to the drag force. 


1-4. Laminar Boundary layer 

Let us find out boundary layer thickness (S), shear stress (j 0 ), local coefficient of 

drag (C* D ) and coefficient of drag (C D ) for the following velocity distribution in the 
laminar boundary layer: 



d& 
2 dx 


dD = wpU 2 dd , 


Z)(* 1 , x 2 ) = wpU 2 [0(x 2 )-^(-X 1 )] 


C i] 

terms. For example, in a wide range of flat plate flows, Cj ~ 10" -10" . This means that, the wall shear 

l| 

stress r w is relatively much less than the dynamic pressure (l/2)pi7 in the flew, e.g., 0.001-1%. 
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1 . — = 2 
U 


f y ] 

f y ' 

UJ 

UJ 


2. 

U 3 {Sj 


3. — = 2 [ — 

u U 


Ifz 

2 U 


\ 3 


)-{y) 

3 

+ 

fU 

J UJ 


UJ 


Prandtl's velocity distribution. 


A U 

4. — = sin 
U 


f n 

V 2^y 


u 


Casel. Velocity distribution: — = 2 
(i) Boundary layer thickness: 

We know, 


|Yi (V 2 

UJ 


j 


t 0 _d 


pU dx 

Substituting the value of u/U, we get 

d 


W 1 -±)dy 

M U) 


o _ 


pU 2 dx 


d 

dx 


8( o 2 V 

2 y y ' 


l 


pu ■ 


0 V 

sf 


S S : 


J 


1- 


r o 2 v 

d d 2 


v 


y 


dy 


1 


ov 


2y_4y 2 , 2 y 
S 8‘ 


+ 


8 : 


y + 2y 


3 


<9 _ 


pu- 


d 

dx 

d 


C-.2 „_.3 


I 


ov 


2 Z _5/ 
d d 2 d 2 


d 2 S 3 

A \ 


S' 


dy 


J 


y 


pU 2 dx 


S A 

y 


+ 

28 3 S 2 4 8 3 5 8‘ 


J 

5 V 


dy 

iS 


yjo 


pU dx 


d ( 5 1 ^ 

' d--d + d--d 

3 5 J 


v 


d f 2 Vi 

U5 J 


dx 
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T 0 = P U Z X 


d 


( 


dx 


\ 


— S 

V 15 j 



dS_ 

dx 


Also, according to Newton's low of viscosity, 


*o=M 


r dif 

\ d yj,= o 


But 


u=U 


2 A 

?Z_ Z_ 

S 5 2 


du 


And 


dy 


= U 


(2 2 y 


\ 


V 


8 8 : 


J 


du^ 


\ d yj y =o 

Substituting this value in (ii), we get 

2 juU 


f 


= U 


V 


8 


-0 


A 


J 


, U being costant 

2JJ_ 

8 


= 


8 


Equating the values ofT 0 , given by Eqn (1.11) and ( 1. 12 ), we get 


2 dS T.uU 

— pU = — 

1 5 dx 5 

_ dS 15 juU 15 ju 

o = — = 

dx pU 2 pU 


S.dS = dx 

pU 

Integrating both sides, we get 


g2 ^ 

— = — — x + C,( where C = constant of integration) 
2 pU 
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At 


x = 0, £ = 0 :.C=0 


S 2 _15 px 
2 ~ pU 


S= M^ = 5.48 ,/K 


P U 


P U 


5 = 5.48 r XXX = 5.48. 1 X 


p Uxx 


Re 


Where Re y — 


pUx 

P 


or <5 = 5.48 




(ii) Shear stress, z G : 


From Eqn. (1.12), we have 


= ML But S = 5.48 

s 


X 



X 



2fjU$r x 

5.48x 


0.365 


pU 


x 



( Hi) Local co-efficient of drag, C D : 


(1.14) 


T o = 0.365 fRg~ x 

X 



* 

D 


pu 2 

2 


(Where C D — local co-efficient of drag) 

Equating the two values of t 0 , given by equation (1.14) and ( 1.9) we get 


[FLUID MECHANICS (1-B)] 


PAGE XvCH.2 


EXTERNAL FLOW - THEORY OF BOUNDARY LAYERS 


c* - 


0.365 

x 


pU 


hence C D = 


0.73 


( iv)Coefficient of drag, C D : 

p 

We know that, C D = - — — — 

-pAU 2 

2 


Where 


F D = | r 0 x B x <±x 


F D = f 0.365^^ 

oL x 


x B x dx 


=1 


0L 


0.365 


pU 
x \ 


pUx 


B 




\ 


x B x dx 


Re. = 


pUx 
B ) 


= 0.365 fp[/J^- x -j= xBxdx- 0365 pU J^^-x Bjx 2 xdx 


1 

1/ 

hxB 

2 

V B 

1/ 

L/2 J 


= 0365 pU 1^-xB^L 
B 


F d =0.73 pUB^ 


pUL 


P 


0.15pUB 


pUL 


• • C o ~ 


P 


-pAU 2 

2 


( Where A - area of plate = Lx B, L and B being width of the plate respectively ) 
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0.73 pUB \ /?UL - 

V m _ 1.46// pUL 

-pxhxBxU 2 pLU V M 
2 


c D = 12 ^ = 1.46 = 

VpL/7 V P UL 

Table 1-1 shows the values of S ( boundary layer thickness), C D ( local co-efficient of 

drag), C L) (average co-efficient of drag) in terms of Reynolds number (Re) for various 
velocity profiles/distributions. 

Table 1.1. Shows the values of 3, C D and C D in terms of Re 


S.No. 


Velocity Profiles 


s 

c* 

^ D 

Cd 

5.48x 

0.73 

^x 

1.46 

V Re L 

4.64x 

0.646 

V R ^ 

1.292 

V Re L 

5.84x 

V Re x 

0.686 

1.372 

V Re L 

4.795x 

0.654 

1.3 lx 

yl^ Q x 


1?1 

CD 

5x 

sl^x 

0.664 

1.328 

V Re £ 


3 


4 


— = 2 

U 


r y \ (y ) 2 

\3J fSy 


U _ 5 

U ~2 


fy) 1 fyV 


\° J 




— = 2 
V 


/ \ / \3 / \4 

I v \ ( v ( y 


y 

\3j 


y_ 

vSj 


+ 


j 


u 

— = sin 

U 


7i y^ 


v 2 8 j 

3 Blasius Results Re < 3.2^10 5 
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Example 1 . 10 . 

The boundary layer thickness at a distance of 1 mfrom the leading edge of a flat 
plate kept over zero angle of incidence to the flow direction is 1 mm. The velocity 
outside the boundary layer is 25 m/s. The boundary layer thickness at a distance of 
4 mis (i) 4 mm, ( ii) 2 mm, ( iii) 1 mm. 

Select the correct answer. Assume that the boundary layer is entirely laminar 

Sol. Free stream velocity, U = 25 m/s 

The boundary layer thickness atX] - 1 m, S ; - 1 mm 

The boundary layer thickness at a distance of 4 m, d 2 : 

Thickness of boundary layer is given by, 







or 


8 2 =2 mm (Ans.) 


Example 13-15. 
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Air is flowing over a smooth flat plate with a velocity of 12 m/s. The velocity 


u (y\ /vY 
profile is in the form: - = 2 - 


\ U J 


y_ 


The length of the plate is 1.1 m and width 0-9 m. If laminar boundary layer exists 
up to a value of Re =2 x 10 5 and kinematic viscosity of air is 0.15 stokes, find: 

(i) The maximum distance from the leading edge upto which laminar 
boundary layer exists, and 

( ii ) The maximum thickness of boundary layer. 

Solution: 

/ \ / \ 2 

u 


M 


A 

UJ 


UJ 


Velocity profile (distribution): — = 2 

Velocity of air, U—12m/s 
Length of plate, L-l.lm 
Width of the plate, B=0.9 

Reinhold's number upto which laminar boundary exists, Re = 2x10 s 

Kinematic viscosity of air, V=0.15 stokes = 0.15 x 10' 4 m 2 /s 

(i) the maximum distance from the leading edge upto which laminar boundary 

layer exists, x: 

n 12xx 

or 2x10 5 = — 
v 


D Ux 
Re = — 


0.15x10 


-4 


2x10 x0.15 x 10" 


x - 


12 


= 0.25m 


( ii) the maximum thickness of the boundary layer, 3: 

for the given velocity profile, the maximum thickness of the boundary layer is given 

5.48x 
by, 8 = -== 

8 = = 0.00306m or 3.06mm (Ans.) 


a/2x10" 
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13-5. Turbulent Boundary Layer: 



As compared to laminar boundary layers, the turbulent boundary layers are thicker. 
Further in a turbulent boundary layer the velocity distribution is much more uniform, 
than in a laminar boundary layer, due to intermingling of fluid particles between 
different layers of the fluid. The velocity distribution in a turbulent boundary layer 
follows a logarithmic law i.e. u ~ log y, which can also be represented by a power law 
of the type: 


u 

U 




1 7 5 

, where n = — (approx.) for Re<10 but >5x10 


i 


u 

U 


(y) 

UJ 


This is known as one-seventh power law. 

Eqn. (1-40), however, cannot be applied at the boundary itself because at 


y = o, 




ydy j 


1 


1 -6 


- — U8 1 y 1 =oc. This difficulty is circumvented by considering the 


velocity in the viscous laminar sublayer to be linear and tangential to the seventh-root 
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profile at the point, where the laminar sub layer merges with the turbulent part of the 
boundary layer. 

Blasius suggested the following relation for viscous shear stress: 


r 0 - 0.0226pU : 


,( 


P 


^4 


pVS, 

( for Re ranging from 5 x l(f to 10 7 ) 

Let us now find the values of 5, z () , C and C d for the velocity distribution given by 


eqn. (13-40) i.e. 


u 

U 


V 

Jy 


Boundary layer thickness, S: 


u 


Substituting the value of — in Von Karman integral eqn. (13-6), we have 


d 


pJJ 1 dx 


d 


s 

f- 



dy 

l u{ 

V u) 



pU 2 dx 


1 

f 

1 'N 


M 

l- 


dy 

J oUJ 

V 

v^J 

J 



_ d 

fi 

r y) 

7 

fy\ 

pU 2 

dx 








8 


9 " 

T 0 

_ d 

7 

y 1 

7 

Y 7 

pU 2 

dx 

8 

i 

9 

2 



5 1 


s 1 


dy 


d_ 

dx 


5_ 

1 9~ 


l_dS_ 
72 dx 


[In the expression above, the limits have been taken from 0 to S instead of S' to S 
since the laminar sub layer (S') is very thin] 


[FLUID MECHANICS (1-B)] 


PAGE vr CH.2 


EXTERNAL FLOW - THEORY OF BOUNDARY LAYERS 


7 


dS 


r 0 = — pU z — 


72 


dx 


Also, 


r 0 - 0.0225 pU 


t± 

pU8 


Now equating the eqns. ( 13-42 ) and ( 13-41 ), we have 

i 


7 


dS 




t 0 — — pU z — = 0.0225 pU 


Or 


72 dx 
7 dS 


p 


\ 


pus. 


72 dx 


= 0.0225 


f jT\ 4 
KPUj 


1 2 

x — — , (canceling pU on both sides ) 

8 4 


5 4 dS — 0.0225.x — x 
7 


72 ( p V 


\PUj 


dx 


8 4 d8 = 0.2314 


p V4 


dx 


\PU 

Integrating both sides, we have 


4 - ( p ^ 

-8 4 =0.2314 

5 \pUj 


x + C ( Where C = constant of integration) 


Let boundary layer be assumed to be turbulent over the entire length of plate. 
Hence, at x=0, d=0 :. C = 0 


:.-8 4 =0.2314 
5 


( jT\ 4 
KpUj 


x Qr 8 4 =(5/4x0.2314) 


P ^ 


KPUj 


x x 


£ = [5/4x0.2314]; 


[puj 


X X' 
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.. Ac 1 


8 = 0.371 


M 


pUx j 


x 5 xx 5 =0.0371 


1 ^ 


8 = 


0.37 lx 

i 

)s 


(Re,): 

(iii) shear stress, r 0 : 


v Re ,y 


1 

5 0.37 lx 

xi = 


(Re,) 


r 0 = 0.0225/?f/ : 


ypUSj 


Substituting the value of d from Eqn. (13.43), we get 

i 




t 0 =0.0225 pU 


\ 


p 


pU 


0.37 lx 


Re^s j 


0.0225 rr ^ 

r 0 = T f*J‘ 

0.3714 
1 P 


n < Re 7j] 

1 

4 =0.0288/?t/ 2 

(Re,)5 

pUx 


1 

X 

Q 

i 


Re A . pUx 


Or 


*0 = 


pU 1 0.0576 

- x - 


( Re Js 

(in) Local co-efficient of drag, C D : 


, pU 2 0.0576 

We know r„ x 


(Re,) 
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Also 


t o= c *d x \p u1 


Now equating the eqns. (13.44) and (13.9), we have 

0.0576 


* 1 tt2 pU 0.0576 

C D *-pU = -^-x r 


(ReJ 


or C D = , 


( Re A 


(iv) Drag force, F D : 

The total drag force (F D ) on one side of the plate of width B and length L is 
given by, 

L 

F d = Jr 0 xBx<7x 
o 


FD= |p|l xM 576 xBx(fa = j^l x ft0576. xBx<fc 


0 


(ReJ 


L ~JJ 2 

2 


0 


f pUx 

v M J 


pU 2 0.0576 

x 


pU 


f 7 . pU 2 0.0576 _ 

JC 5 = — x 5 

5 5 

-xr 

J o i 

° f yoUV 

4 


\ M 


J 


-\L 


pU 0.0576 

x 


5 4 

-xflx— ( l ) 5 

1 A v ’ 


pU 0.072 D T 

x x B x L 

o i 

5 



5 

f 

M J 


l F J 


^ pU 1 0.072 „ r 

Or F n = x x B x L 

u o 1 

(Rej5 

(v)Co-efficient of drag, C D: 
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1 ? 

We know, F D = C D x — pAV 2 
Now equating eqns. (13.10) and (13.46), we have 


C n x — px B x LxU 2 
D 2 2 


0.072 


xfixL 


Or 


C D = 


0.072 


(Re.U 

This is valid for 5 x 10 5 < Re L < 10 7 

For Reynolds number between 10 7 and 10 9 the following relationship suggested by 
Prandtl and Schlichting holds good. 

0.455 


C D — 


\ 2.58 


(log io Re J 

13-6. Total Drag due to Laminar and Tujjmlent Layers 

When the leading edge is not very rough, the turbulent boundary layer does not 
begin at the leading edge, it is usually preceded by the laminar boundary layer. The 
point of transition from laminar to turbulent layer depends upon the intensity of 
turbulence. The distance x c (Fig. 13-7) of the transition from the leading edge can be 
obtained from critical Reynolds number which normally ranges from 3 x l(f to 3 x 10 6 . 
V 


.Laminar boundary 
layer 


Transition 


Turbulent boundary . 
layer 



[FLUID MECHANICS (1-B)] 


PAGE VV CH.2 


EXTERNAL FLOW - THEORY OF BOUNDARY LAYERS 


13.7. Drag due to laminar and turbulent boundary layers. 

Drag force (F D — F) for the turbulent boundary layer can be estimated from the 
following relation: 

^turb (F turb. ) total turb . ) x c 

Where ( F uture ) tota i = the drag which would occur if a turbulent boundary extends 
along the entire length of the plate, and 

(Ftrub)Xc = the drag due to fictitious turbulent boundary layer from the leading 
edge to a distance X c . 

Let us assume that the plate is long enough so that Reynolds number is greater 
than 10 , then the turbulent drag is given by, 

„ 0.455 pU 2 , D 0.072 pU 2 

F,„, — — x xLxB x xx xB 


turb ■ n D , \ 2.58 

(log 10 Re L ) 


(Re = ) 


Where: 

L = length of the plate, 

B = width of the plate, and 
U - free stream velocity. 

The laminar boundary layer prevails within the length x c and its contribution to 
drag force is given by 


la min ar 


1.328 p\J 2 „ 1.328x „ pU ' 

x — xi xB = — , x B x 


2 

F total F laminar F turbulent 


_ 1.328* pU 

Frnai = xBx + 


J7 - 
total 


1.328*. 


VRe ( 


0.455L 


(log jo Re L ) 


2.58 


pU 2 0.072* pU 

x B x — -VxBx- — 


(Re c ) 


1/5 


+ 


0.455L 


0.072x„ 


x/R® 7 (log 10 Re L ) 2 - 58 (Re c ) 1/5 


BpU 
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Re {pUx / //) x 

Also, — - = -7 — f = — Or 

Re L (pUL/ ju) L 


X c = 


Re c L 
Re r 


Substituting the value ofx c in eqn. (iii), we have 

0.455 0.072 Re 


r - 

total 


1.328^/Re 


0.8 


Re, 


+ 


(log 10 Re L ) 


2.58 


(R ej 


1/5 


LBpU 


Assuming that transition occurs at Re c = 5x10 s 


F tolal =C D d-pAU 2 =C D d^fT 


(where Cd = average co-efficient of drag.), Equating the above two equations 

0.455 1700 


C D = 


(logjoReJ 258 Re t 


Example 13-24 

A submarine can be assumed to have cylindrical shape with rounded nose. Assuming 
its length to be 50 m and diameter 5.0 m, determine the total power required to 
overcome boundary friction if it cruises at 8 m/s velocity in sea water at 20° C (p = 
1 030 kg/m), v=lxl0' 6 m 2 /s. 


Sol. Length of submarine, L — 50 m, Diameter of submarine, D = 50 m 

Velocity of submarine, U = 8 m/s, Density of sea water, p= 1030 kg/m 3 
Kinematic viscosity of sea water, v= 1 x 10' 6 m/s 

Total power required to overcome boundary friction, P: 


Reynolds number, Re L = 


UL 

v 


8x50 

lxlO ' 6 


= 4xl0 8 


The length over which boundary layer will be laminar is given by, 


Ux c m5 

— = 5x10 or 


x = 


5x10 5 xi/ 

U 
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x ■ 


5xl0 5 xlxl0' 6 

8 


= 0.0625m 


This being very small, contribution to total drag from laminar boundary layer is 
negligible; hence Cd is given by, 

0.455 0.455 


C D ~ 


(log 10 Re L ) 2 58 [log 10 (4xl0 8 )] 2 ' 58 


= 0.001765 


Area, 


A = ttDL = n x 5x 50 = 785 .4m" 


.-. Drag Force, F D =C D x-pAU 2 = 0.001765x^x1 030x785.4x8 2 = 456902^ 

2 2 

Hence total power required to overcome boundary friction, 


P=W kW = 456902X8 


1000 


1000 


365.52 kW (Ans) 


13 . 7 Separation 

What causey separation ? 

The increasing downstream pressure slows down the wall flow and can make it go 
backward-flow separation, dpfdx > O adverse pressure gradient, flow 

separation may occur. c ^ ! < ^ , favorable gradient, flow is very resistant to 

separation. Previous analysis of BL was valid before separation. 
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dp 


> O 


U 



( d ) Critical adverse 
gradient: 



gradient: 


dU 

djc 


< O 


dp 

dx 


> G 


No separation, 
PI in the flow 


Zero slope 
at the wall: 

Separation 


Backflow 
at the wall: 

Separated 
flow region 



Note : 1. Due to backflow close to the wall, a strong thickening of the BL takes place 
and BL mass is transported away into the outer flow 

2. At the point of separation, the streamlines leave the wall at a certain angle. 

Boundary Layer Separation and its Control 


[FLUID MECHANICS (1-B)] 


PAGE i \ CH.2 


EXTERNAL FLOW - THEORY OF BOUNDARY LAYERS 



Fig* 1 + 8 Separation of the boundary layer 
and vortex formation at a circular cylinder 
(schematic)* S = separation point 


Notes: 

1. D to E, pressure drop, pressure is transformed into kinetic energy. 

2. From E to F, kinetic energy is transformed into pressure. 

3. A fluid particle directly at the wall in the boundary layer is also acted upon by the 
same pressure distribution as in the outer flow (inviscid). 

4. Due to the strong friction forces in the BE, a BL particle loses so much of its 
kinetic energy that is cannot manage to get over the “pressure gradient” from E to 
F. 

5. The following figure shows the time sequence of this process: 

a. reversed motion begun at the trailing edge 

b. boundary layer has been thickened, and start of the reversed motion has moved 
forward considerably. 

c. and d. a large vortex formed from the backflow and then soon separates from the 
body. 
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In a flowing fluid when a solid body is immersed, a thin layer of fluid called the 
boundary layer is formed adjacent to the solid body. The forces acting on the fluid in the 
boundary layer are: 

1- Inertia forces, 2- Viscous forces and 3- Pressure force. 


— When the pressure gradient in the direction of flow is negative 


^<0 

dx 


i.e. when the 


pressure decreases in the direction of flow, the flow is accelerated. In this case, the 
pressure force and inertia force add together and jointly tend to reduce the effect of 
viscous forces in the boundary layer. This results in a decrease in the thickness of 
boundary layer in the direction of flow, as a consequence of which there are low losses 
and high efficiencies in accelerating flows. 


— When the pressure increases in the direction of flow 


^>0 

dx 


the pressure forces 


acts opposite to the direction of flow and further increases the retarding effect 
of the viscous forces. Subsequently the thickness of the boundary layer increases 
rapidly in the direction of flow. If these forces act over a long stretch, the 
boundary layer gets separated from the surface and moves into the main stream. 
This phenomenon is called separation. The point of the body at which the 
boundary layer is on the verge of separation from the surface is called point of 
separation. 


Consider a flow of fluid over a curved surface as shown in Fig. 1.8. 


— As the fluid flows round the surface (the area of flow decreases ) it is accelerated 
over the left hand section until at point B the velocity just outside the boundary is 
maximum and the pressure is minimum (as shown by the graph below the surface). 
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Fig.13.8. Separation of boundary layer. 

Edge of the y boundary layer Separations stream line A and B the pressure 


gradient is negative, As long as 


^<0 

dx 


, the entire boundary layer moves forward. 


Beyond B (i.e. along the region BCDE), the area of flow increases and hence velocity of 
flow decreases; due to decrease of velocity the pressure increases ( in the direction 

dp 
dx 


of flow) and hence the pressure gradient 


is positive i.e. 


^>0 

dx 


. The value of the 


velocity gradient 


du 

dy 


at the boundary is zero at the point C, this point is known as 


a separation point ( the boundary layer starts separating from the surface because 
further retardation of flow near the surface is physically impossible). Large turbulent 
eddies are formed downstream of the point of separation. The disturbed region in which 
the eddies are formed is called turbulent wake. 

The flow separation depends upon factors such as: 

1- The curvature of the surface; 
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2- The Reynolds number of flow; 

3- The roughness of the surface. 

The velocity gradient, for a given velocity profile, exhibits the following characteristics 
for the flow to remain attached, get detached or be on the verge of separation. 


1- 


du^ 

v3>Vv- o 


is positive Attached flow (the flow will not separate) 


2 - 


^ du^ 


aJ v=0 


is zero The flow is on the verge of separation 


3- 


1 ' du^ 


A y v =o 


is negative separated flow. 


- Boundary layer separation is unstable, inefficient process and entails large losses 
due to appreciable eddying region. 

-Separation occurs in the following cases: 

1. Diffusers , 

2. Open channel transitions, 

3. Pumps, 

4. Fans, 

5. Aerofoil, 

6. Turbine blades etc. 


Examples of BL Separations (two-dimensional): 

Features: The entire boundary layer flow breaks away at the point of zero wall 
shear stress and, having no way to diverge left or right, has to go up and over the 
resulting separation bubble or wake. 
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Fig. 2.10- Stagnation point flow, after H. Pcrttinger (1939). (a) free stagnation- 
point flow without separation, (b) retarded stagnation— point flow, with separation 

(a) . Plane stagnation-point flow: no separation on the streamlines of symmetry (no 
wall friction present ), and no separation at the wall ( favorable pressure gradient) 

(b) .Flat wall with right angle to the wall: flow separate, why? 

2. Diffuser flow: 


1 . Plane wall(s) 


Tliin wall 



Fig. 2.9. Flo* in t widening channel (diffuser! (a separation at both diffuser 
wadis, (b) auction of the boundary lam at the upper diifuseT wall, (c) suction at 
both diffuser walls (alter L Prandtl; 0. Tiecjens (1931)) 
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(b) 


3. Turbulent Boundary Layer 

Influence of a strong pressure gradient on a turbulent flow 

(a) a strong negative pressure gradient may re-laminarize a flow 

(b) a strong positive pressure gradient causes a strong boundary layer top thicken. 
Examples ofBL Separations ( three-dimensional ) 

Features: unlike 2D separations, 3D separations allow many more options. 

There are four different special points in separation: 

(1) . A nodal Point, where an infinite number of surfaces streamline merged 
tangentially to the separation line. 

(2) . A saddle point, where only two surface streamlines intersect and all others 
divert to either side 

(3) . A focus, or spiral node, which forms near a saddle point and around which an 
infinite number of surface streamlines swirl. 

(4) . A three-dimensional singular point, not on the wall, generally serving as the 
center for a horseshoe vortex. 


Profile NACA 0015 



Separation regions in corner flo* between airfoifs. [After Gersien (1959). ) 
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1. 3D separations on a round-nosed body, at angle ofattack 



Tfireenjiirncnsionsil separation on a round -nosed body at angle of attack, first 
described by Legendre Point A ns a nodal auachment point, point $ h a saddle pointy and 

point F U a focus of separation. 


Methods of preventing the separation of boundary Layer: 

Following are some ofjhe methods generally adopted^ to retard or arrest the flow 
separation: 

1. Streamlining the body shape. 

2. Tripping the boundary layer from laminar to turbulent by provision of surface 
roughness. 

3. Sucking the retarded flow. 

4. Injecting high velocity fluid in the boundary layer. 

5. Providing slots near the leading edge. 

6. Guidance of flow in a confined passage. 

7. Providing a rotating cylinder near the leading edge. 

8. Energizing the flow by introducing optimum amount of swirl in the incoming 
flow. (Note: Refer Example 13.29 also) 

Example.13.29 

Explain what is meant by separation of boundary layer. Describe with 
sketches the methods to control separation. 

Sol. 

When a solid body is immersed in a flowing fluid, a thin layer of fluid called 
boundary layer is formed, adjacent to the solid body. In this thin layer of fluid, the 
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velocity varies from zero to free stream velocity in the direction normal to the solid 
body. Along the length of the solid body, the thickness of the boundary layer increases. 
The fluid layer adjacent to the solid surface has to do work against surface friction at 
the expense of kinetic energy. This loss of kinetic energy is recovered from the 
immediate fluid layer in contact with the layer adjacent to the solid surface through 
momentum exchange process. Thus the velocity of the layer goes on decreasing. 
Along the length of solid body, at a certain point a stage may come when the boundary 
layer may not be able to keep sticking to the solid body, if it cannot provide kinetic 
energy to overcome the resistance offered by the solid body. In other words, the 
boundary layer will be separated from the surface. This phenomenon is called the 
boundary layer separation. The point on the body at which the boundary layer is on 
the verge of separation from the surface is called point of separation. 

Methods to control separation: 

1. Motion of solid boundary: 

By rotating a circular cylinder lying in a stream of fluid, so that the upper side of 
cylinder where the fluid as well as the cylinder moves in the same direction, the 
boundary layer does not form. However on the lower side of cylinder where the fluid 
motion is opposite to that of cylinder separation would occur (Fig. 13.9). 
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Fig. 13. 10. Injection fluid into boundary layer 



Fig. 3. 11 slotting wing 



Fig. 13. 12 suction of fluid from boundary layer 


2. Acceleration of fluid in the boundary layer: 

This method of controlling separation consists of supplying additional energy to 
particles of fluid which are being retarded in the boundary layer. This may be 
achieved either by injecting the fluid into the region of boundary layer from the 
interior of the body with the help of some available device as shown in Fig. 13.10 or 
by diverting a portion of fluid of the main stream from the region of high pressure to 
the retarded region of boundary layer through a slot provided in the body (fig. 13.11) 

3. Suction of fluid from the boundary layer: 

In this method, the slow moving fluid in the boundary layer is removed by suction 
through slots or through a porous surface as shown in the Fig. 13.12. 

4. Streamlining of body shapes: 

By the use of suitably shaped bodies the point of transition of the boundary layer 
from laminar to turbulent can be moved downstream which results in the reduction of 
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the skin friction drag. Furthermore by streamlining of body shapes, the separation may 
be eliminated. 


1.1.5. Comments on "Boundary layer, pressure gradient, and separation" 

- Fluid in the BL in the region near wall is continuously subjected to (net) shear 
stress (acting on fluid in the -x direction) from the start of the leading edge. (See 
Fig. 2a.) 

Inviscidlrrotational external flow 




Fig. 2. Boundary layer as a layer of viscous/shear/rotational flow. Two views of shear: fluid 
element and CV. External flow is often inviscid and irrotational (if in fact the fluid is static and the 
body is moving). 


1.1.6. Wake: Region of high shear/friction behind a body. Vortex shedding. 

• Fluid that was subjected to shear/viscous force in the BL is convected 
downstream from the body, forming another viscous/shear layer called wake. 
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• Bulk of rotational fluids from the upper and lower BL ’s often leave the upper 
and lower surfaces alternately, resulting in a phenomenon called vortex 

shedding. 

• Vortex shedding: shedding of BL ’s vortices. 

• Similar to BL, one can think of mass and momentum deficits in the wake. 

1.1.7. Aerodynamics: Aerodynamic forces and moments (lift, drag, etc.): 

• They are simply resultants of stress distribution (pressure + friction) due to 
fluid on solid surface. So long as there is no large-scale separation: 

• Lift is due mainly to pressure and can be estimated by inviscid effect/flow. 

• Drag is due mainly to wall shear stress (skin friction drag) and is the result of 
viscous effect. If there is large-scale separation: 

• Loss of lift. 

• Pressure distribution is such that it causes relatively larger amount of drag 
(pressure drag) than wall shear stress (skin friction drag). 

• Total drag D = Skin friction drag (Df, T w ) + Pressure drag (D p ). 

• Streamline Body and Bluff Body (and The In-Between): [ Phenomenon wise, 
whether a body is streamline or bluff depends upon how it appears to the free 
stream, not to us. Hence, a well-designed, streamlined airfoil is a bluff body if 
it moves in a free stream upright or at very high angle of attach] 

• Streamline body: No large-scale separation: 

• Drag of streamline body is due mainly to wall shear stress (skin friction drag). 

• Bluff body: Large-scale separation: 

• Drag of bluff body is due mainly to pressure distribution (pressure drag). 
However, this is not to mean that wall shear stress is not important. Recall 
that wall shear stress causes the large-scale separation in the first place. 

• Certainly, there are the In-Betweens. 


1.8 AXIOM 

The foregoing may now be generalized into the following simple axiom: 
"Acceleration of real fluids tends to be an efficient process, deceleration an 
inefficient one ". 

Accelerated motion occurs along the surface of front part of submerged object. 
(l.U is accompanied by a favorable pressure gradient which serves to stabilize the 
boundary layer and thus minimize energy dissipation. Decelerated motion is 
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accompanied by an adverse pressure gradient which tends to promote separation, 
instability, eddy formation and large energy dissipation. 

1.9 Secondary Flow 

Another consequence of wall friction is the creation of a flow within a flow - a 
secondary flow superposed on the main primary flow. 

Fig.21 shows a cross section through a river bend. At the mean line fluid particles 
are in equilibrium under the centrifugal and pressure force. Near walls the pressure 
forces dominate the centrifugal forces. And considering the application of 
Bernoulli's equation, one expect that the flow velocity at the outer side is lower than 
that at the inner side with a pressure distribution opposite to the velocity 
distribution. 


Therefore is originated the flow shown, which called secondary flow. The 


scouring the outer side, thus producing the well-known meandering characteristics 
of natural stream. 




secondary flow serves to deposit material at the inside of the bend and to assist in 
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Another example of secondary flow is shown in Fig. 22, where shown the horseshoe- 
shaped vortex that is produced by a projection from a boundary surface. The 
velocity profile along the wall leads to stagnation pressure at a being larger than 
that at B. This pressure difference maintains a downward secondary How from A to 
B, thus inducing a vortex-type of motion, the core of the vortex being swept 
downstream around the sides of the projection. This principle is used on some 
aircraft wings, the devices ( called vortex generators) being used to draw higher 
energy fluid down to the wing surface to forestall large-scale separation. 
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SHEET (1) 

"BOUNDARY LAYER THEORY'" 


— = 2 

M 


(y) 

U 

UJ 


W 


, 8 being 


boundary layer thickness. Calculate the following: 

a) Displacement thickness, 

b) Momentum thickness, and 


c) Energy thickness. 

2. If velocity distribution in laminar boundary layer over aflat plate is assumed to be given 
by second order polynomial u = a + by + cy , determine its form using the necessary 
boundary conditions. 


3. The velocity distribution in the boundary layer is given by — = 


mV 

~ 8 ) 


Calculate the 


following: 

a) Displacement thickness, 

b) Momentum thickness, 


c) Shape factor 

d) Energy thickness and, 

e) Energy loss due to boundary layer if at a particular section, the boundary layer 
thickness is 25 mm and the free stream velocity is 15 m/s. If the discharge through the 
boundary layer region is 6 m 3 /s per meter width, express this energy loss in terms of 
meters of head. Take p = 1.2 kg/m 3 . 
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4. In the boundary layer over the face of a high spillway, the velocity distribution was 
observed to have the following form: 


u 

U 


r \°- 22 


\° J 


The free stream velocity U is 20 m/s and boundary layer thickness 5 cm at a certain 
section. The discharge is 5 m /s per meter length of spillway. Calculate displacement 
thickness, energy thickness and loss of energy up to section under consideration. 

5. The free stream velocity U at a certain section was observed to be 30 m/s and a boundary 
layer thickness of 60 mm was estimated from the velocity distribution measured at the section. 
The discharge passing over the spillway was 6 m /s per meter length of spillway. Calculate: 

a) The displacement thickness, 

b) The energy thickness, and 

c) The loss of energy up to the section under consideration. 

6. A smooth plate 2 m wide and 2.5 m long is towed in oil (sp. gr. = 0.8) at a velocity of 1.5 
m/s along its length. Find the thickness of boundary layer and shear stress at the trailing 

edge of the plate. v oil =\Q m Is. 


7. A plate 450 mm x 150 mm has been placed longitudinally in a stream of 
crude oil (specific gravity 0.925 and kinematic viscosity of 0.9 stoke) which flows with 
velocity of 6 m/s. Calculate: 

a) The friction drag on the plate, 

b) Thickness of the boundary layer at the trailing edge, and 

c) Shear stress at the trailing edge. 


8. Air is flowing over aflat plate 5 m long and 2.5 m wide with a velocity of 4 m/s at 15°C. If 

3 5 2 

p= 1.208 kg/m and v = 1.47 x 10' m /s, calculate: 


a) Length of plate over which the boundary layer is laminar, and thickness of the 
boundary layer (laminar), 


b) Shear stress at the location where boundary layer ceases to be laminar, and 
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c) Total drag force on both sides on that portion of plate where boundary layer is 

laminar. 


9. Atmospheric air at 20° C is flowing parallel to a flat plate at a velocity of 2.8 m/s. 
Assuming cubic velocity profile and using exact Blasius solution estimate the boundary 
layer thickness and the local co efficient of drag (or skin friction) at x = 1.2 mfrom the 
leading edge of the plate, Also find the deviation of the approximate solution from the 
exact solution. Take the kinematic viscosity of air at 20°C=15.4 x 10' 6 m 2 /s. 


10. Airflows over a plate 0.5 m long and 0.6 wide with a velocity of 4 m/s. The velocity profile 


is in the form — = sin 
U 


f n y ^ 
\2~d J 


If p = 1.24 kg/m 3 and v = 0.15 x 10 ~ 4 m 2 /s, calculate: (i) 


Boundary layer thickness at the end of the plate, (ii) Shear stress at 250 mm from the leading 
edge, and (iii) Drag force on one side of the plate. 


11. Find the ratio of friction drag on the front half and rear half of the flat plate kept at zero 
incidence in a stream of uniform velocity, if the boundary layer is laminar over the whole 

plate. 


12. Air at standard conditions is flowing over aflat plate which is 1 m long and 0.3 m wide. 
The flow is uniform at the leading edge of the plate. The velocity profile in the boundary 

as the free stream velocity is U = 30 m/s. Assume 

that the flow conditions are independent of Z Using control volume abed, shown by 
dashed line, calculate the mass flow rate across the surface ab. [Density of air may be 


— = 2 

M 


(y) 

U 

W 


U J 


taken as 1.23 kg/m 3 , refer to Fig. 1.1] 



Fig. 1.1 
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DIMENSIONAL AND MODEL ANALYSIS 


Introduction 

Dimensional analysis is a mathematical technique which makes 
use of the study of the dimensions for solving several 
engineering problems. 

Each physical phenomenon can be expressed by an 
equation giving relationship between different quantities; such 
quantities are dimensional and non-dimensional. Dimensional 
analysis helps in determining a systematic arrangement of the 
variables in the physical relationship, combining dimensional 
variables to form non-dimensional parameters. It is based on the 
principle of dimensional homogeneity and uses the dimensions of 
relevant variables affecting the phenomenon. 

The method of dimensional analysis is used in every field of 
engineering, especially in such fields as fluid dynamics and 
thermodynamics where problems with many variables are 
handled. This method derives from the condition that each term 
summed in an equation depicting a physical relationship must 
have same dimension. By constructing non-dimensional quantities 
expressing the relationship among the variables, it is possible to 
summarize the experimental results and to determine their 
functional relationship. 


Uses of dimensional analysis: 

The uses of dimensional analysis may be summarized as follows: 

1. To test the dimensional homogeneity of any equation of fluid 

motion. 

2. To derive rational formulae for a flow phenomenon. 

3. To derive equations expressed in terms of non-dimensional 
parameters to show the relative significance of each 

parameter. 
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4. To plan model tests and present experimental results in a 
systematic manner; thus making it possible to analyze the 
complex fluid flow phenomenon. 

Advantages of dimensional analysis: 

Dimensional analysis entails the following 
advantages: 

1. It expresses the functional relationship between the variables 
in .dimensionless terms. 

2. in hydraulic model studies it reduces the number of variables 
involved in a physical phenomenon, generally by three. 

3. By the proper selection of variables, the dimensionless 
parameters can be used to make certain logical deductions 
about the problem. 

4. Design curves, by the use of dimensional analysis, can be 
developed from experimental data or direct solution of the 

problem. 

5. It enables getting up a theoretical equation in a simplified 
dimensional form. 

6. Dimensional analysis provides partial solutions to the 
problems that are too complex to be dealt with 

mathematically. 

7 . The conversion of units of quantities from one system to 
another is facilitated. 

2. Dimensions 

The various physical quantities used in fluid phenomenon can 
be expressed in terms of fundamental quantities or primary 
quantities. The fundamental quantities are mass, length, time and 
temperature, designated by the letters M, L, T, 9 respectively. 
Temperature is especially useful in compressible flow. The 
quantities which are expressed in terms of the fundamental or 
primary quantities are called derived or secondary quantities, ( e.g., 
velocity, area, acceleration etc.). The expression for a derived 
quantity in terms of the primary quantities is called the 
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dimension of the physical quantity. 

A quantity may either be expressed dimensionally in M-L-T 
or F-L-T system ( some engineers prefer to use force instead of 
mass as fundamental quantity because the force is easy to 
measure). The table- 1 gives the dimensions of various quantities 
used in both the systems. 

Fluid Mechanics and Heat Transfer and their Dimensions 


Quivtitits used u fluid MeduolcG uid Hnl 
Transfer ind (hetr ItimeiHians- 


S A/* 




Quantity 

Af -kr 

w.r. 

Xijften 

1. 

(■} K^ndutnul 

M 

FL -, T l 

2, 

Len|ifa. k 

L 

L 

3. 

Time. T 

T 

T 

4. 

Ci> (jrfmwlrk OdmUiIcS 
A re A 

L a 

L s 

5 

VoflUBlfifV 

L? 

L 3 

$ 

MeflWH iiKflja. 

LC 1 

Id 

7. 

jt) EintfCJlJr Quantities 
Linear velocity, - n, ¥ h U 

LT-' 

LT" 1 

ft. 

Arif □ Lit velociiy r oouiiui^ spud. JV 

T" 1 

T‘ ( 

ft 

Acceleration. 0 

LT -1 

LT -2 

10. 

Angular accdermnvi, a 

T' 1 

T-* 

II 

Discharge. fi 


L'T' 1 

l± 

! Grpi-ity, ^ 

LT -1 

LT -2 

13. 

Kinematic vixtosiLy, v 

L^T 1 

L*?' 1 

14. 

Stream i it notion ,'4^ cLmliLieJm.. F 

L?T' 1 

Lh 4 * 1 

15- 

Vonkity. ft 

T"“ 

r 

T -‘ 

l* 

(rf"J Dynorole QmnllEles 
Forte. F 

MLT~ Z 

f 

n. 

Densiriy; p 

MX-' 3 

FL" 4 T= 

ifl. 

fipthiLl'-t vf'c cjlr k L, il' 

ML' J T‘" J 

FL" 3 

19. 

Dynamic vi^oiily. |1 

ML' 1 T ~ 1 

FL _Z T 

1ft. 

Ptcssute, shear stress, r 

ML‘ I T' S , 

FL' 3 

2k 

Modulus of tlu'kky. £, AC 

ML" 1 T _z 

FL" 1 

22_ 

McifTKTilUm 

MLT~‘ 

FT 

23. 

Angular morTHMLluin ur Q'iaeru±iii <sf napnncpicuirri 

ML^T" 1 

F1.T 

24. 

Wert, cncrEj,. £ 

ML 1 T -2 

FL 

25. 

Tongue, T 


FL 

26. 

Pome. P 

ML Z T J 

FLT' l 

27. 

irj Tbeirnodyn antic QuanElEks 
T" cm pe ra Lu rc 

■ ft 

0 

2G. 

Thermal conducitvicy 

MLT’ST 1 

ft" 1 a -1 

29. 

E-ti du 1 py pci- uniL rruLi 

L’T 2 

l*t' z 


onsiam 

l “ r-V 1 

L*T" V 1 

M. 

Entropy 

ML*T" 3 B" 1 

FLe‘ 

32- 

Imcrnal per vm\ nlMB 

L^r- 1 

'L’T' 1 j 


tkal transfer 

ML*T"* ! 

RJT" 1 | 
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Example. 3-1. Determine the dimensions of the following quantities: 
( i ) Discharge, ( ii ) Kinematic viscosity, 

( iii ) Force and ( iv) Specific weight. 

Sol. (i) Diwharft*’ = are* x velocity 

= D x ^ = ^ = L 3 T~ ‘ (A ns.) 

T T 

(rl) JCintfJiioliO viscosity <v) = ^ 


where u is by ; 


it* 


T shear Sires* _ fpTce/ grea 


ilu 

dy 


k'i 


, M *T* 

m a&s x acceleration _ 1 


ML 


M 

LT 


area x 1/T 
~ML ’T' 1 P 

It- i 


l ! x- t 2 r 2 x 




1 volume L 5 


•. t — I — = /Tt * 1 (Am.) 

Kinematic viscosity l>) “ p j^// 3 

(Jti) Fore* = m*s* * acceleraboo 

„ AJ x „ MLT~ 2 (Ans.) 

time 2 . T 

weight ' 1 _ siV^' 2 (Ans.) 

(rv) Specific **tghl * vollmlC - volume Lr 


3-3. Dimensional Homogeneity 

A physical equation is the relationship between two or more 
physical quantities. Any correct equation expressing a physical 
relationship between quantities must be dimensionally 

homogeneous (according to Fourier's principle of dimensional 
homogeneity ) and, numerically equivalent. Dimensional 

homogeneity states that every term in an equation when reduced to 
fundamental dimensions must contain identical powers of each 
dimension; A dimensionally homogeneous equation is applicable 
to all systems of units. In a dimensionally homogeneous equation, 
only quantities having the same dimensions can be added, 
subtracted or equated. Fet us consider the equation: P— y . h 
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Dimension!* of L.H.S- - ML T 

. f n u c - AfL ri T' ! xL = ML 7”* 
DicncTisions of R.ri-^ -ml 

_ . . r i il t - j^nifrncirtni Of R.cC.Ji 

tit L.M.5 


Equation P= y.h is dimensionally homogeneous; so it can be 
used in any system of units. 

Applications of Dimensional Homogeneity: 

The principle of homogeneity proves useful in the following 

ways: 

1. It facilitates to determine the dimensions of a physical 
quantity. 

2. It helps to check whether an equation of any physical 
phenomenon is dimensionally homogeneous or not. 

3. It facilitates conversion of units from one system to another. 

4. It provides a step towards dimensional analysis which is 
fruitfully employed to plan experiments and to present the 
results meaningfully. 

3-4. Methods of Dimensional Analysis 

With the help of dimensional analysis the- equation of a 
physical phenomenon can be developed in terms of dimensionless 
groups or parameters and thus reducing the number of variables. 
The methods of dimensional analysis are based on the Fourier's 
principle of homogeneity. The methods of dimensional analysis 
are: 

1. Rayleigh 's method 

2. Buckingham 's it- method 

3. Bridgman 's method 

4. Matrix-tensor method 

5. By visual inspection of the variables involved 

6. Rearrangement of differential equations. 
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Here only first two methods will be 
dealt with. 

3-4-1. Rayleigh's Method 

This method gives a special form of relationship among the 
dimensionless group, and has the inherent drawback that it does 
not provide any information regarding the number of 
dimensionless groups to be obtained as a result of dimensional 
analysis. Due to this reason this method has become obsolete and 
is not favoured for use. 

Rayleigh's method is used for determining the expression for a 
variable which depends upon maximum three or four variables 
only. In case the number of independent variables becomes more 
than four, then it is very difficult to find the expression for the 
dependent variable. 

In this method a functional relationship of some variables is 
expressed in the form of an exponential equation which must be 
dimensionally homogeneous. Thus ifX is a variable which depends 
on X„ X 2 , X 3 ...Xn, the functional equation can be written as: 

X=f(X b X 2 , X 3 , X n ) 3-1 

In the above equation X is a dependent variable, while X„X 2 , 
X 3 , ... X n are independent variables. A dependent variable is the 
one about which information is required while independent 
variables are those which govern the variation of dependent 
variable. Eqn. (3-1) can also be written as: 

x = (x 1 ,x 2 ,x 3 , x n ) 3-2 

Where C is a constant and a, b, c, ...., n are the arbitrary powers. 
The values of a, b, c, ... n are obtained by comparing the powers of 
the fundamental dimensions on both sides. Thus the expression is 
obtained for dependent variable 

Example 3.2: Find an expression for the drag force on smooth 
sphere of diameter D, moving with a uniform velocity V in a fluid 
density p and dynamic viscosity p. 
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Sol. The drag force F is a function of 
(»') Diameter D, (ii) Velocity V, 

(Hi) Fluid density p, and (fv) Dynamic viscosity p. 

Mathematically, F=f(D, Kp.p) or F=C(tf. V i .p c .\i d ) ,„(1) 

where C is a non-dimensional constant. 

Using M-L-T system the corresponding equation for dimensions is: 

MLT ~ 2 = [CL* . (17*' ')* . (ML 3 ) ( . (ML' Y 

For dimensional homogeneity the exponents of each dimension on both sides of the 
equation must be identical. Thus 

For M : 1 = c + d . (f) 

For L : 1 = a + b - 3c - d ...(/i) 

ForT: -2-b-d 

There mfour unknowns (a, b, c, d) but equations are three in number. Therefore, it is not 
possible to find the values of a, b, c and d. However, three of them can be expressed in terms 
of fourth variable which is most important. There the role of viscosity is vital one and hence a, 
b, c are expressed in terms of d (i.e. power to viscosity) 


c - i - d ...from (0 

b = 2-d ...from («T) 

Putting these values in (i), we get 
a=l-b + 3c + d=l-2 + d + 3(l-d) + d 

= )- 2 + d+3-2>d + d = 2-d 
Substituting these values of exponents in eqn. (1), we get 
F = C [D*~ d . V 2 ' 




r d - 

= c[ & V 2 p {D~ d . V* . p-** . = c|_ p D 2 j 


» p D 1 V 3 <1> 




2. Buckingham 's n-Method/Theorem 

When a large number of physical variables are involved 
Rayleigh's method of dimensional analysis becomes increasingly 
laborious and cumbersome. Buckingham's method is an 
improvement over the Rayleigh's method. Buckingham designated 
the dimensionless group by the Greek capital letter 71 (Pi). It is 
therefore often called Buckingham 7t -method. The advantage of this 
method over Rayleigh's method is that it let us know, in advance of 
the analysis, as to how many dimensionless groups are to be 
expected. 
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The Buckingham's re-theorem states as follows: 

"If there are n variables ( dependent and independent 
variables) in a dimensionally homogeneous equation and if these 
variables contain m fundamental dimensions (such, as M, L, T, 
etc.), then the variables are arranged into (n-m) dimensionless 
terms. These dimensionless terms are called 7t -terms. 

Mathematically, if any variable X h depends on independent 
variables, X 2 , X 3 , X 4 , ... X n the functional equation may be written as 

XU(X 2 , X3, X4...X,,) ...(3-3) 

Eqn. (3-3) can also be written as 

f=(X h X 2 ,X 3 ...Xn)=0 ...(3-4) 

It is a dimensionally homogeneous equation and contains n 
variables. If there are m fundamental dimensions then according 
to Buckingham's 7t-theorem, it [eqn. (3-4)] can be written in terms 
of number ofn-terms (dimensionless groups) in which number ofn 
-terms is equal to (n-m). Hence eqn. (3-4) becomes as 

fl tt?, TTj, Ttn-m) — 0 ...(3-5) 

Each dimensionless 71 -term is formed by combining m 
variables out of the total n variables with one of the remaining 
(n-m) variables i.e. each n -terms contain (m+ 1 ) variables. 
These m variables which appear repeatedly in each of 7t-terms 
are consequently called repeating variables and are chosen from 
among the variables such that they together involve all the 
fundamental dimensions and they themselves do not form a 
dimensionless parameter. Let in the above case X 2 , X 3 , and X 4 
are the repeating variables if the fundamental dimensions m (M, 
L, T) =3. Then each term is written as: 

n l =X$>.X%' .Xj .X x 

n 2 = (. X ° 2 * . X% . X % . X 5 ) 

• ► 

• * 3 — • ■ *n) n 
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where ci[- bj- Ci ; a 2 - b 2 - C 2 etc. are the constants, which 
are determined, by considering dimensional homogeneity. These 
values are substituted in eqn. (3-6) and values of 7ii, 7 T 2 , H 3 , ... 7t n _m 
are obtained. These values of it's are substituted in eqn. (3-5). The, 
final general equation for the phenomenon may then be obtained by 
expressing anyone of the 7t-terms as a function of the other as 


Selection of repeating variables: 

The following points should be kept in view while selecting m 
repeating variables: 

1. m repeating variables must contain jointly all the fundamental 
dimensions involved in the phenomenon. Usually the fundamental 
dimensions are M, L and T. However, if only two dimensions are 
involved, there will be 2 repeating variables and they must contain 
together the two dimensions involved. 

2. The repeating variables must not form the non-dimensional 
parameters among themselves. 

3. As far as possible, the dependent variable should not be selected 
as repeating variable. 

4. No two repeating variables should have the same dimensions. 

5. The repeating variables should be chosen in such a way that 
one variable contains geometric property ( e.g. length, I; diameter, 
d; height, H etc.), other variable contains flow property (e.g. 
velocity, V; acceleration, a etc.) and third variable contains fluid 
property (e.g. mass density p; weight density y, dynamic viscosity 
p, etc.). The choice of repeating variables, in most of fluid 
mechanics problems, may be: 


Tlj - ^ (% 2 > *£ 3 . ^4* ■ 

tt-2 ~ ^ 7 * 4 * ■** 71 n-m ) 


(3.7) 


( i ) / t Kp 
im) LV, \i 


00 d , V, p 
(iv) d , V, p. 
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The procedure for solving problems by Buckingham's n- 
theorem is outlined in the Example 3-4 below: 

Example 3-4. The resistance R experienced by a partially submerged 
body depends upon the velocity V, length of the body l, viscosity of 
the fluid p, density of the fluid p and gravitational acceleration g . 
Obtain a dimensionless expression for R. 

Sol, Step. 1. The resistance R is a function of: 

( i ) Velocity V, ( ii ) Length l 

( in) Viscosity p, and ( iv) Density p, 

(v) Gravitational acceleration g. 

Mathematically, R = f(V, l, p,p,g) or fi (R,V, l,p,p,g) = 0 

Total number of variables, n = 6, m is obtained by writing 
dimensions of each variable as: 

R = MLT 2 , V = LT p = ML 'T p = ML 3 , g = LT 2 . 

Thus the fundamental dimensions in the problem are M, L, T and 
hence m-3. 

Number of dimensionless Pi-n- terms -n - m=6-3=3. 

Thus three n-terms say iti, it 2 , 713 are formed. 

The eqn. may be written as f l (Tt l ,n 1 ,Tt ?) ) = 0 

Step 2: Selection of repeating variables: Out of six variables 
R,V ,l,p,p,g three variables (as m-3) are to be selected as 

repeating variables. R is a dependent variable and should not be 
selected as a repeating variable. Out of the remaining five variables 
one variable should have geometric property, second should have 
flow property and third one should have fluid property; these 
requirements are met by selecting l, V and p as repeating variables. 
The repeating variables themselves should not form a 
dimensionless term and must contain jointly all fundamental 
dimensions equal to m i.e. 3 here. Dimensions of l, V and p are 

L, LT , ML' and hence the three fundamental dimensions exist 
in l, V and p and also no dimensionless group is formed by them. 
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Step 3: Each n -term (=m + 1 variable) is written as given in eqn. 
(3-6), i.e. 

tc , = /“ i . \y ri> \ . p c i _ ft 
TCy — l‘*- . V&l p c 2 . p. . 

TC 3 = 1°* . V*3 . ^ 


Step 4: Each n -term is solved by the principle of dimensional 
homogeneity, as follows: 

ftj-term : 

71, = l a ' . V h > . p C '.R 

M°L°T° = L a ' . (LT ')*' . (ML'Y 1 • ( MLT~ 2 ) 

Equating the exponents ofM, L and T respectively, we get 

For M : 0 = c, + 1 

For L : 0 = a l + b x — 3c, + 1 

For T : 0 = - - 2 

c, = - l ; Z>, = - 2 

and a, = -^,+3c,-l=2-3-l = - 2 


Substituting the values of a h b h and Ci in ic h we get 


7C 2 — term: 


7t, = l~ 2 . V~ 2 . p~ 1 . R = 


n 2 = 1°^ . V b 2 . p c 


R 


/ 2 V 2 p 


M°L°T° = L a * . (J_T~ ') b 2 , (MIS' 3 ) Cz • ') 


Equating the exponents ofM, L and T respectively, we get 


For M : 
For L. : 
For T : 

and 


O = r 2 + i 

0 = a 2 + i> 2 - 3 c 2 — I 
0 = - - 1 
c 2 — — 1 , b 2 = — \ 

~~ "I” 3c y + 1 ™ 1 ■ — 3+1 - — 1 
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JC 2 — l ! - v 1 . p 1 . |i = 

7E3— term; 

= J"» - V ** . . g 

— ZJ"* . 1 ) 6i * i MLT 3 > c ’ - <Z-J 2 ) 

Equating the exponents of A / , £, and f respectively, we get 

For M : O = c 3 

For 1- : O = t 2 3 + ir 3 — 3c 3 + 1 

For T : O == — — 2 

= O ; t *3 = — 2 

and ^ fc 3 + 3 e 3 — 1 = 2 + 0—1= l 

Substituting the values of fr 3 . andc 3 in 7T 3 * we get 



Step 5: Substitute the values of n h n 2 , 713 . The functional 
relationship becomes 


b 


R 


R 

V z p 


/Vp 


* a V 2 p 


-*( 


/Vp 

pVI 


]- 

) 


_£s_ 

V 




The above step has been made on the postulate that reciprocal 
of pi-term and its square root is non-dimensional. 


R = Tv 2 p <t> 



The resistance R is thus a function of Reynolds number 



and Froude 's number 



EXAMPLE 3-5: Friction in a Pipe 

■ Consider flow shown in Fig.; V is the average speed across 
the pipe cross section. The flow is fully developed, which 
means that the velocity profile also remains uniform down 
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the pipe. Because of frictional forces between the fluid and 
the pipe wall, there exists a shear stress tw on the inside 
pipe wall. The shear stress is also constant down the pipe in 
the region. We assume some constant average roughness 
height, s, along the inside wall of the pipe. In fact, the only 
parameter that is not constant down the length of pipe is the 
pressure, which must decrease (linearly) down the pipe in 
order to “push” the fluid through the pipe to overcome 
friction. Develop a nondimens ional relationship between 
shear stress t w and the other parameters in the problem. 



Solution 


List of relevant parameters : 

t w = f(V, e, p, p. D) n = 6 

primary dimensions 

t w V e 

p p D 

{m'L-'r 2 } {L't -1 } {L 1 } 

{m'L -3 } {m'L _, r'} {L 1 } 

Reduction: 

7 = 3 —^k=n-j= 3 

Repeating parameters: 

V. D. and p 

n , = r/'DV' -» n '= p ;> 

Modified IT j . n i , modified 

8t k . 

^2 = Darcy friction factor = / 


n 2 = p r 2 £>V 2 
n 3 = er’DV 3 -» 


n 2 

n 3 


pVD 

= = Reynolds number = Re 

A 

e 

= — = Roughness ratio 


=> f = 


pV 2 



3-4-3. Limitations of Dimensional Analysis 
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Following are the limitations of dimensional analysis: 

1. Dimensional analysis does not give any clue regarding the 
selection of variables. If the variables are wrongly taken, the 
resulting functional relationship is erroneous. It provides the 
information about the grouping of variables. In order to decide 
whether selected variables are pertinent or superfluous 
experiments have to be performed. 

2. The complete information is not provided by dimensional 
analysis; it only indicates that there is some relationship 
between parameters. It does not give the values of co-efficient 
in the functional relationship. The values of co-efficient and 
hence the nature of functions can be obtained only from 
experiments or from mathematical analysis. 


MODELANALYSIS 

Introduction 

In order to know about the performance of the hydraulic 
structures (e.g. dams, spillways etc.) or hydraulic machines (e.g. 
turbines, pumps etc.) before actually constructing or 
manufacturing them, their models are made and tested to get the 
required information. The model is the small scale replica of the 
actual structure or machine. The actual structure or machine is 
called Prototype. The models are not always smaller than the 
prototype, in some cases a model may be even larger or of the same 
size as prototype depending upon the need and purpose 


Advantages of model testing 

The following are the advantages of model analysis: 

1. The model tests are quite economical and convenient ( because the 
design, construction and operation of a model may be changed 
several times if necessary, without increasing much expenditure, 
till most suitable design is obtained). 
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2 . With the use of models the performance of hydraulic 
structures/hydraulic machines can be predicted in advance. 

3. While designing a particular portion of the structure if clear cut 
analytical and reliable method is not available then in such cases it 
becomes absolutely necessary to know about the safety and 
reliability of such parts which is possible by means of model testing. 

4. Model testing can be used to detect and rectify the defects of an 
existing structure which is not functioning properly. 

Applications of the model testing 

Following are the important fields where applications of the model 
testing is of great use: 

1. Civil engineering structures such as dams, spillways, weirs, canals 
etc. 

2. Scour in rivers, irrigation channels. 

3. Turbines, pumps and compressors. 

4. Design of ships and submarines. 

5. Aeroplane, rockets and missiles. 

6. Tall buildings (to predict the wind loads on buildings, the stability 
characteristics of the buildings and airflow patterns in their vicinity). 

Next, in order to determine the characteristics of a full-scale 
device through model tests, besides geometrical similarity, 
similarity of dynamical conditions between the two is also 
necessary. When the above dimensional analysis is employed, if 
the appropriate non-dimensional quantities such as Reynolds 
number and Froude number are the same for both devices, the 
results of the model device tests are applicable to the full-scale 
device. 

When the characteristics of a water wheel, pump, boat or 
aircraft are obtained by means of a model, unless the flow 
conditions are similar in addition to the shape, the characteristics 
of the prototype cannot be assumed from the model test result. In 
order to make the flow conditions similar, the respective ratios of 
the corresponding forces acting on the prototype and the model 
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should be equal. 

7. Forces Influencing Hydraulic Phenomena (e.g.), the forces 
acting on the flow element: 

The forces acting on the flow element are due to gravity Fq, 
pressure Fp, viscosity F v , surface tension F r ( when the prototype 
model is on the boundary of water and air), inertia F h and 
elasticity F E . The forces can be expressed as shown below. 

The forces which may affect/influence the flow characteristics of a 
problem are: 

1. Inertia force (F,) 

It always exists in the fluid flow problem (and hence it is 
customary to find out the force ratios with respect to inertia force). It 
is equal to the product of mass and acceleration of the flowing fluid 
and acts in the direction opposite to the direction of acceleration. 

2. Viscous force (F T ) 

It is present in fluid flow problems where viscosity is to play an 
important role. It is equal to the product of shear stress ( T) due to 
viscosity and surface area of the flow. 

3. Gravity force (F G ) 

It is present in case of open surface flow. It is equal to the 
product of mass and acceleration due to gravity, 

4. Pressure force (F P ) 

This type of force is present in case of pipe-flow. It is equal to the 
product of pressure intensity and cross-sectional area of the 
flowing fluid. 

5. Surface tension force (F s ) 

It is equal to the product of surface tension and length of surface of 
the flowing fluid. 

6. Elastic force (F E ) 

It is equal to the product of elastic stress and area of 
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the flowing fluid. 


gravity force 
pressure force 

viscous force 

surface tension force 

inertial force 

elasticity force 


Fa — mg = pL?g 
Fp = (A p)A = (A p)L 2 

Fr = ^y = = » vL 

F, = TL 

F f = mat = pi? = pL 4 T ~ 2 = ptfL? 
F e = KL 1 


Since it is not feasible to have the ratios of all such 
corresponding forces simultaneously equal, it will suffice to 
identify those forces that are closely related to the respective flows 
and to have them equal. In this way, the relationship which gives 
the conditions under which the flow is similar to the actual flow in 
the course of a model test is called the law of similarity. In the 
following section, the more common force ratios which ensure the 
flow similarity under appropriate conditions are developed. 

6. Similarity 

To find solutions to numerous complicated problems in hydraulic 
engineering and fluid mechanics model studies are usually 
conducted. In order that results obtained in the model studies 
represent the behavior of prototype, the following three similarities 
must be ensured between the model and the prototype. 


1. Geometric similarity; 

2. Kinematic similarity, and 

3. Dynamic similarity. 

1. Geometric similarity: 
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For geometric similarity to exist between the model and the 
prototype the ratios of 

corresponding lengths in the model and in the prototype must be same 
and the included angles between two corresponding sides must be the 
same. Models which are not geometrically similar are known as 
geometrically distorted models. 

Let L m = length of model, 

H m = height of model, 

D m - diameter of model, 

A m ~ area of model, 

V m = volume of model, 

and L p , B p , H p> D p , A p and V p = corresponding values of the prototype. 
Then, for geometric similarity, we must have the relation; 

A- = = 

L p -B f -HpD p L ' 

where L f is called the scale ratio or the scale factor. 

Similarly A r = area ratio = -^-11 

k t> 

y 

aod V r ~ volume ratio = = tf 

'p 

2. Kinematic similarity: 


Kinematic similarity is the similarity of motion. If at the corresponding 
points in the model and in the prototype, the velocity or acceleration 
ratios are same and velocity or acceleration vectors point in the same 
direction, the two flows are said to be kinematically similar. 


Similarly 


(Y,) m 

iPdrn 

( a l )p 


(yfhn 

(V 2 ) p 


= V n velocity ratio 


(<h)r 


a n acceleration ratio 


The directions of the velocities in the model and prototype should be 
same.The geometric similarity is a pre-requisite for kinematic 
similarity. 

3. Dynamic similarity: 
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Dynamic similarity is the similarity of forces. The flows in the 
model and in prototype are dynamically similar if at all the 
corresponding points, identical types of forces are parallel and bear 
the same ratio. In dynamic similarity, the force polygons of the two 
flows can be superimposed by change in force scale., 


Let (F,) m = inertia force at a point in the model, 
[F v ) m - viscous force at the point in the model, 
(^)m = gravity force at the point in the model, 


an d (Ff, (F v ) p , ( F g ) p = corresponding values of forces at the corresponding points in 
prototype, 

Then for dynamic similarity, we have 
(Fjim 

Wp = (FJ P ~ (FJ b 


■ = F r , force ratio 


i-p 


The directions of the corresponding forces at the corresponding 
points in the model and prototype should also be same. 

8. Dimensionless Numbers and their Significance 

The dimensionless numbers (also called non-dimensional 
parameters) are obtained by dividing the inertia force (which always 
exists when any mass in motion ) by viscous force or gravity force or 
pressure force or surface tension force or elastic force. The 
important dimensionless numbers: 


1. Reynolds Number Re: It is defined as the ratio of the inertia force 
to the viscous force. 

2. Froude number (Fr): It is defined as the square root of the 
ratio of the inertia force and the gravity force. 

3. Euler's Number (Eu): It is defined as the square root of the 
ratio of the inertia force to the pressure force. 

4. Weber Number (We): It is defined as the square root of the 
ratio of the inertia force to the surface tension force. 

5. Mach Number (M): It is defined as the square root of the ratio 

of the inertia force to the elastic force. 


[FLUID MECHANICS (1-B)] 


Page 19-ch3 


DIMENSIONAL ANALYSYS 


Table : Dimensionless groups/numbers 


SI. 

Dimensionless 



Aspects 

No. 

numbers 

Symbol 

Group of 
variables 

Significance 

Field of application 

i. 

Reynolds 

number 

Re 

pVL 

P 

Inertia force 
Viscous force 

Laminar viscous flow in 
confined passages (where 
viscous effects are significant) 

2. 

Froude 's 
number 

Fr 

V 

TQ 

Inertia force 
Gravity force 

Free surface flows (where 
gravity effects are important) 

3. 

Euler's number 

Eu 

V 

Vp 

Inertia force 
Pressure force 

Conduit flow (where pressure 
variations are significant) 

4. 

Weber's 

number 

We 

V 

Vo/pL 

Inertia force 
Surface tension 

Small surface waves, capillary 
and sheet flow (where surface 
tension is important) 

5. 

Mach's number 

M 

yfH’ 

' K/p 

Inertia force 
Elastic Force 

High speed flow (where 
compressibility effects are 
significant). 


Similarity Laws 

To ensure dynamic similarity between the model and prototype 
it is necessary that the ratio of the corresponding forces acting at 
the corresponding points in the model and prototype be made 
equal. It implies that dimensionless numbers should be same for 
the model as well as the prototype; this condition is difficult to be 
satisfied for all the dimensionless numbers. Hence models are 
designed on the basis of the force which is dominating in the flow 
situation. The laws on which the models are designed for 
dynamic similarity are called model or similarity laws; these are: 

1. Reynolds model law, 

2. Froude model law, 

3. Weber model law, and 

4. Mach model law. 

5. Euler model law, 

1. Reynolds model law, 

In flow situations where in addition to inertia, viscous force is the 
only other predominant force. Where the compressibility of the 
fluid may be neglected and in the absence of a free surface, e.g. 
where fluid is flowing in a pipe, an airship is flying in the air (Fig. 
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3.2) or a submarine is navigating under water, only the viscous 
force and inertia force are of importance. The similarity of flow in 
the model and its prototype can be established if Reynolds number 
is same for both the systems. 

This is known as Reynolds law and according to this law 

Re Model = Re Prototype 



inertia force F t pv 2 L? L.vp JLv 

viscous force F v pvL, p — v ~ 

which defines the Reynolds number Re, 

Re = Lv/v 

Consequently, when the Reynolds numbers of the prototype and 
the model are equal the flow conditions are similar. 

Following are some of the phenomena for which Reynolds 
model law can be a sufficient criterion for similarity of flow in the 
model and the prototype: ( i) Motion of air planes, 

(i) Flow of incompressible fluid in closed pipes, 

(ii) Motion of submarines completely under water, and 

(iii) Flow around structures and other bodies immersed 
completely under moving fluids 

For similar dynamical flow conditions, the ratio of 
corresponding forces acting at corresponding points in the model 
and prototype should be equal. The ratios of forces are 
dimensionless numbers. It means for dynamic similarity between 
model and prototype, dimensionless numbers should be same for 
model and prototype. 
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Since the submarine has to overcome the viscous resistance, there 
has to be dynamic similarity between the model and the prototype; 
which implies equality of Reynolds number. 

2. Froude model law, 

When the resistance due to the waves produced by motion of a 
boat ( gravity wave) is studied, the ratio of inertia force to gravity 
force is important: 

inertia, force P ¥ pv^lJ 2 v z 

gravity force pMJ*& G^ 

In general, in order to change v2 above to v as in the case for Re, 
the square root ofu2/gL is used. This square root is defined as the 
Froude number Fr, 

Fr — — - — 

When the gravitational force can be considered to be the only 
predominant force which controls the motion in addition to the inertia 
force, the similarity of the flow in any two such systems can be 
established if the Froude' s number for both the systems is the same. 
This is known as Froude Model Law. 

Some of the phenomena for which the Froude model law can be 
sufficient criterion for dynamic similarity to be established in the model 
and the prototype are: (1) Free surface flows such as flow over 
spillways, sluices etc.; (2) Flow of jet from an orifice or nozzle; (3) 
Where waves are likely to be formed on the surface; (4) Where fluids 
of different mass densities flow over one another. 

If a test is performed by making the Fr of the actual boat (Fig. 
3.3) and of the model ship equal, the result is applicable to the 
actual boat so far as the wave resistance alone is concerned. This 
relationship is called Froude ’s law of similarity. For the total 
resistance, the frictional resistance must be taken into account in 
addition to the wave resistance. Also included in the 
circumstances where gravity inertia forces are important are 
flows in an open ditch, the force of water acting on a bridge pier, 
and flow running out of a water gate. 
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Fr = — — 



3. Weber model law, and 

When a moving liquid has its face in contact with another fluid or 
a solid, the inertia and surface tension forces are important: 


In this case, also, the square root is selected to be defined as 
the Weber number, We. 


In a fluid system where surface tension effects predominate in 
addition to inertia force, the dynamic similarity is obtained by equating 
the Weber number for the model and its prototype, which is known as 
Weber model law. 

According to this law: 


Weber model law is applied in the following flow 
situations: 

Flow over weirs involving very low heads; Very thin sheet of liquid, 
flowing over a surface; Capillary waves in channels; Capillary rise in 


inertia force F t pv 2 L 2 



surface tension F T TL 


T 


We = VyfpLpT 


We = u/pL/T 
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narrow passages; Capillary movement of water in soil. 

4. Mach model law. 

When a fluid flows at high velocity, or when a solid moves at high 
velocity in a fluid at rest, the compressibility of the fluid can 
dominate so that the ratio of the inertia force to the elasticity force 
is then important. 

inertia force F t pv 2 L 2 _ v 2 $ 

elastic force F E KL K)p a 2 


Again, in this case, the square root is selected to be defined as 
the Mach number M, 

M = v/a 

M < 1, M = 1 and M > 1 are respectively called subsonic flow, 
sonic flow and supersonic flow. When M = 1 and M < 1 and M > 
1 zones are coexistent, the flow is called transonic flow. 



Boeing 747: full length, 70.5 m; full width, 59.6 m; passenger capacity, 498 
persons; turbofan engine and cruising speed of 891 km/h (M = 0.82) 


When in any fluid system only the forces resulting from 
elastic compression are significant in addition to inertial force, then 
the dynamic similarity between the model and its prototype may be 
achieved by equating the Mach numbers, which is known as Mach 
model law. 
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M — v/a 


m iltZ 1 ™' e “ dms lhc spced ” f *»* 

(<v) Umfcr-waier testing of WpeJoes, 

wloaif of M m/i, Calcukie tht miini m thtVLr, ' ,**“ M ll ' mg> ' ™ tf « “ 

. . „ ! ^ I°r cakHlmg ,fe jfc „ 

8 'w by, C f =Qm/{Rtf\ " J * 

wm (w m water) * m tint" ^ ** ™ ,fr) W ^ Wo * f ^ ^ specific night of 
Sdi. Linear scale ratio. £,^20 


Submerged area of the model, ft 1 

Lct^ih of the model, = $ m 
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Total drug of model. * 20 N 
Velocity of model. V m * 1-5 mA 

Let A r 1^, R r V p be the corresponding values for the prototype. 
Kinematic viscosity of sea water, 

v_«ir,= CK>l stokes = 001 cmVs = 0-01 x 10' 4 mVs 

Total drag on the prototype 

(l) Analysts of model : 


Reynolds number. 


* 01 - 


V-ia 

V_ 


1-5x8 


Also 


001 x 10 
00735 
C/= ^ 
c 0073 S 

(I2x 


77=12x10* 
1 * 


= 000282 


Frictional or skin resistance of the model is given by : 

r 


F}-±C /P AV l 


-i x0O0282xf ^2 |x5 x 1-5 1 

- 15 86A * 

f wh « c k v.ave rcsistancccxpencnccd by the model and 1 
L (*/•)- = frictional or skin resistance experienced by the model. 

20 = (*„),. *15-86 
(A.) w = 20- 15-86 = 4 14 N 


(ft) Analysis of prototype: 

Since resistance to wave formation exists, condition will be dynamically similar if the 
Froudt S numbers are equal. 


it. 




V r~V m * ^ 15 * ^ = 6 - 70R t>/* 
Reynolds number, Rt • 

i v i v f 


S^?=.(h73x,0» 
0-01 x 10 4 




0-0735 


(V tp m gj 


= 20 


^ = 8 x 20= 160m 
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FriccL-anal or %kia mkruuicc of (he pfOGOEyps is givtn by. 

iF dr c " ^ c »Pp A f l £ *5*™ 14 a* HIS* 2000 k S- 70 S 2 


■ 51637 N 


[ 2000 m* 


WflT^ rtwlwicc ft* for dir prototype mn be evaluated fry misfying the faUowiJig 
ctmiSjiion for dynamic simJliu^Ly: 


Bui 

and 


ff- |_f *» 

■ v ’l l^**JL 


j£ , 


P,“P* 


v ^ 


= k 

A, 


■ -From cht Gqutvjitern ofFnemde’i number^ 


<*.!,= '*(^1 


5< (^Oi. = < iOV K J- 14 = 3 J ! 20 W 


Hence coral drug an the pmlvylypr. 

A, = * <fy # * 3 3 120 4 5 I 657 - 84777 N (Ant) 


5. Euler model law, 

In a fluid system where pressure forces alone are the controlling forces in addition to the 
inertia force . the dynamic similarity is obtained by equating the Euler number for both the 
model and its prototype. This is known as Euler model law. According to this law: 

(^ w )model - (^ u )pnxcxypc 

If, V m = velocity of fluid in model, 

p„ = pressure of fluid in model, 
p m = density of fluid in model, 

and V p , p p , p p are the corresponding values in prototype, then by substituting these values 
in eqn. ( 34), we get 


fa7(>m fa? V P 

when p m = p /) , (i.e. same fluid in model and prototype) the above equation becomes 


fa fa 

This law is applied in the following flow problems : 

(0 Enclosed fluid system where the turbulence is fully developed so that viscous forces 
are negligible and also the forces of gravity and surface tensions are entirely absent; 
(it) Where the phenomenon of cavitation occurs. 
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Example , k an aeroplane model of size ~ of us prototype the pressure drop is 

75 kN/ m'. The model is tested in water. Find the corresponding pressure drop in the prototype. 
Take: Density of air = 1-24 kg/ Density of water - 1000 kg/m* 

Viscosity of air = 0-00018 poise; Viscosity of water - 0-0! poise. 


”• Li "' ai ** * «• <top t» model, (ip). = 7 5 kN/n,’- 

Demi* of w«,r. p.= I000kg/m>; V.scos,,, of wuer, =M, p„ ise . 

Densily of air. P,= l 24k W; Viocooi.y of ^-WWOUpoiso ' 

Pressure drop in the prototype, (AP) p • 

similarity between r*’ *^ fore for 

considered. P ypt ’ v cr s numbcr Reynolds number .should be 

Making Reynolds number equal, we get 


or 


or 


or 


or 


cn-m 


or 


or 


PmKLm 




Pp 

h 




K X L. 


Substituting the values, we have 

V m I 24 


x40x— 2£L_ 

V. 1000 X W X 0-000 18 


: 2-755 


Now making Euler's number equal, we get 

Wsliih] 


vl d 

(; 


p~ 

Ap)„ 

i V( Wpp V <H 

V^^ x a/^. ? 7 „ A /i^r 
V <H ^ V 2 ‘ 755x v T ^= 78 - 24 
M 
(ApV 

(4p>, 

Hence pressure drop in the prototype = 1-225 N/m J (Ans.) 


Pp 

>2-, 


- = (78 24) 2 = 6121-5 

- to)* _ 7-5x1000 _ , 
r 61215 — 6121-5 ' N/m = ' 225 N/m* 
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INTERNAL FLOWS 
FLOW IN PIPES AND DUCTES 


For the flow of real fluid in ducts, channels and pipes, many of 
the concepts developed before can be applied with only minor, but 
important, changes needed to account for viscous effects, boundary 
layers, separation and secondary flows. 

The problem of fluid flow in pipelines-the prediction of flowrate 
through pipes of given characteristics, the calculation of energy 
conversions therein, and so forth-is encountered in many areas of 
engineering practice. The subject of pipe flow embraces only those 
problems in which pipes flow completely full; pipes that flow partially 
full, such as sewer lines and culverts. 

The solution of practical pipe flow problems; results from application 
of the work-energy principle, the equation of continuity, and the 
principles and equations of fluid resistance. Resistance to flow in pipes 
is offered not only by long reaches of pipe but also by pipe fittings, such 
as bends and valves, which dissipate energy by producing relatively 
large-scale turbulence. 

4. 1- Steady Flow FUNDAMENTAL EQUA TIONS 

4.1.1- SHEARS Stress and Head Loss 

The effect of friction forces on the boundary surface inside a pipe, can be calculated 
by applying the linear impulse momentum principle, derived before in the form, 

^ R + Fj - ^ (mC 2 - thC to a c [ rcu i ar cylindrical pipe, in which the surface 

abed abed 

shear stress To represent the effect of friction forces on the solid periphery (boundary) of 
passage opposing the direction of fluid motion. Fig. 3.1 
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Let us write the impulse momentum relation along the direction of this 
stream tube to the control volume bounded by sections 1, 2 and the stream tube 
boundary, as result: 


pA — (p + dp) A — r 0 P dl — 



A dl 


dz 

dl 


= (V + dV?A(p + dp) - V 2 Ap 


Where: S = jiD, (perimeter of the stream tube) 

m=pAV=(p+dp)A (V +dV) = const 

Introducing hydraulic radius; 

For an established incompressible flow in a tube of constant cross section, t 0 
is not a function of l, y = const and d (1/ y) =0, and we can write the last 
relation in the form: 



+ z ) 


t 0 dl 
-y^h 

ai 

~Y r h 


Integrating this equation between 1 & 2, one has: 





Tjh - lj) 
yRh 


The left hand side of the above equation represents the drop of the energy 
line A (EL) between points 1&2. The energy line is also called the total head 
line, and the drop of the energy line is the head loss hL 


Pi v i 

7 + 2g„ + Z| 




Z 2_ 
2 g n 


+ Z 2 = A(EL) 


Tjh - lx) 
y R h 


And we can write: 


L<?fh ~ (if 


The forgoing analysis may be similarly applied to any streamtube of the flow, 
e.g. streamtube with radius r and concentric with the axis of cylindrical pipe, 
fig. 3. 2. The frictional stress T will be that exerted on the outer most fluid 
layer of the streamtube by the adjacent (more slowly moving) fluid. By 
analogy, we can write T for T 0 , r/2 for R h , hi for hu. 2 and l for f-h, from 
which: 
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This shows that, in established pipe flow, the shear stress in the fluid 
varies linear with distance from the pipe centerline. As note both of the last 
two relations are valid for both laminar and turbulent flow in pipe as they 
have been developed without regard to the flow regime. 



Example. 1 

Water flows in a 0.9 m by 0.6 m rectangular conduit. The head lost in 60 m of this conduit 
is determined (experimentally) to be 10 m. Calculate the resistance stress exerted between 
fluid and conduit walls. 


Solution 


Since 


, rfll 2 - f) 

hr ~ 

^ yR h 

R h = A/P 

and the data (given or obtained from Appendix 2) are 

y = 9.8 kN/m 3 h L ^ = 10 m / 2 - /, = 60 m 
A = 0.9 X 0.6 = 0.54 m 2 P = 2(0.9 + 0.6) = 3 m 

we have 


10 X 9.8 X 10 3 X 0.54/3 
60 


= 0.29 kPa * 


If the results cited in the preceding problem are obtained for water flowing in a cylindrical 
pipe 0.6 m in diameter, what shear stress is to be expected {a) between fluid and pipe wall, 
and (b) in the fluid at a point 200 mm from the wall? 
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The additional relevant equation is 

r 

and 



d = 0.6 m r = 0.1 m 

See the previous problem for Eq. 7.36, y, h L = h L , l = l 2 — / ( , R h , etc. 
Clearly, 


Rn 


d/ 4; then 

/ 10 X 9.8 X IQ 3 " 

V 60 


0.15 = 0.25 kPa ® 


From the linear variation of r with r (Eq. 7.37), 

t = ^§(0.25 X 10 3 ) = 83.3 Pa = 0.083 3 kPa ® 


4.1 .2- General Energy Equation for Steady Incompressible 
Flow: 

Considering the work energy equation; for the control volume which 
coincides with the walls of streamtube or, full-sized conduit, Fig, 3.3. 

dQ + dW _ dE 
dt dt dt 



Starting from the equation considering the unit mass of substance, which has the 
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form: 

m = PiA l V 1 — p 2 A 2 V 2 in this steady flow. Therefore, 


dQ . 

1 dQ 

"N-m" 


joule 


'ft -lb" 

H = mq 

o 

£5 

II 

3- 1 

*1 

II 

il kg _ 


. kg _ 


_s!ug_ 


Now dividing the equation by g and considering incompressible fluid; p/=p 2 =p, 
noting that: 

~~~f~ = Qy(Ep - E r ) = rhg n (E p - E t ) 


1 ^^^shaft 

rh dt 


SnYp 


SnY T 



Where: Et, is the energy extracted by the turbine per unit weight of fluid 
flowing, [J/N] E p , is the energy added by the pump per unit weight of fluid 
flowing, [J/N =m] 


dW now 

dt 




Pi^tY 7 . 


Inflow 
rh dt 


P± 

Pi 


Pz 

P2 



ft lb 
slug 


The total energy of the system is given by 


= /// i • dm = III (kinetic + potential + internal energies) * dm 


System 


System 


= III + 8nZ + ie ) ' p d V 

System 

U 2 -U ], change of specific internal energy of fluid. Then the work energy 
equation takes the form: 


dE 

dt 


II 


i • (pv * dA) -b 


Jf i • (pv • dA) 

C. S. in 


Which; is known as "first law of thermodynamics for a control volume". For 
which, when we neglect the effect of pump and turbine, it takes the form: 

= Jf y2 + 8n z + ie j(P y • dA I + If [k v2 + SnZ+ iejipv • dA) 

C.S.ou, ' ' C.S-in ' ' 

The left hand side of this equation represent the head loss hff-2 (as 
explained in section 8.10), the flow in both cases apply precisely the same conditions 
( incompressible flow in tubes). And we can write 
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dE 

dt 


\ y2 + SnZ + ie 


(P2^2^2) 


\ V 2 +g„z + ie 


However, m = p 2 \ ' 2 A 2 = p 1 V 1 A 1 ; thus, 




This equation offers proof that head loss is not a loss of total energy but rather a 
conversion of energy into heat, part of which leaves the fluid, the remainder 
serving to increase its internal energy. 


J_ dE 
rh dt 



+ g n z + ie 



g n Z 


ie 


Here head loss is a permissible and useful concept because heat energy leaving 
the flow and energy converted into internal energy are seldom recoverable and are in 
effect lost from the useful total energy (pressure, velocity and potential 
energies). 

Ph + (g n E P ~ g„E T ) + - ~j = (5 y2 + + /e ) “ y2 + SnZ + 'ej 


P 1 
Ti 


Yl 

2g n 


+ z 1 I + 


E P = 


P 2 


. V? 

+ 7, h Z 2 

72 2 8 n 


4- E'r + 


Sn 


(ie 2 


ie. 


Ph) 


This is riot generally true, since the useful Total of energies will include the 
internal energy. Introducing die last equation generates a reasonably (but not 
completely ) general energy equation for steady-incompressible flow, in the 
form: 


(r + ft + *■) + “ ( 


P2 

y 


yl 

7s „ 


z 2 ) 


This equation is derived considering the mean streamline concept therefore; 
correction factors should be applied to allow for nonuniform distribution 
of fluid" parameters at each section; (velocity, pressure, density and 
temperatures). 


Example3: 

A flowrate of 1.42 m 3 /s of water occurs in a streamtube, say a pipe of variable cross sec- 
tion, containing a pump, but no turbine. The pump delivers 300 kW to the flowing fluid. 
Measurements at Lwo points, i.e., 1 and 2 along the streamtube, show that A L — 0.4 m 2 , 
A 2 = 0.2 m 2 , Z\ — 9 m, z 2 = 24 m, p l = 138 kPa, p 2 = 69 kPa. Calculate the head lost 
between sections 1 and 2. 


Solution: 
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Taking y = 9.81 kN/'m 3 and using Eqs. 4.4 and 5.8 produces 

Q = AV, 1.42 = 0.4V, - 0.2 V 2 V, = 3.55 m/s V 2 = 7.1 m/s 

P = QyE P and so E P = (300 X 10 3 )/(9 810 X 1.42) 

= 21.54 J/N = 21.54 N-m/N = 21.54 m 


Apply Eq. 7.46 with E T = 0 to find (recall g n = 9.81 m/s 2 ) 


pi + Yi 

Tl 2 g n 

138 000 . (3.55) 2 


+ z, 


+ E P — 


= (Hz + Yi 

V72 2 g n 


9 810 


+ 


2g„ 


+ 9 + 21.54 = 


69 000 
9 810 


Z* 
(7.1) 2 

2g„ 


+ E r + h Lt2 


+ 24 + h Li _ 2 


h L = 11.7 m 

" 1—2 


4. 1.3- Velocity Distribution and its Significance 

Due to viscous effect (shear stress effect), we have non-uniform velocity 
distribution, characterized by reduced velocities near boundary surface, Fig. 3.4. 
In this case we can't consider the mean streamline concept, in which the 
values on the cross sectional area is considered constant and equal to the 
values on the mean line. We must count for the non-uniform velocity 
distribution by taking the integration over the cross section 



Total kinetic energy (J/’s or ft lb/s) 



dA 


Momentum flux (>J or lb) 



u 2 dA 


For the volume flow rate: 
Q = fd<2 = f VdA 

A A 

For total kinetics energy: 



A 


J* pvt/A 

A 
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AAA A 


For momentum flux: 

MF = J* pdQy = J* pvcMv = J* p v 2 dA = pj* v 2 dA = yj* ^ dA 

However we can use the mean velocity V, and total flow rate Q for 
expressing total kinetic energy and momentum flux by introducing correction 
factor a, p (dimensionless) defined as: 

v 2 ( p 2 \ 

Total kinetic energy = aQy - — = Qy I a - — 

2 8n \ 2 8nJ 

Momentum flux - f3 QpV 


For momentum flux and total kinetics energy 


II v 3 dA 

- 1 A 

and 

J J v 2 clA 

n ^ a 

V 2 ff 

JJ "" 

A 

13 V ff 

J I v dA 

A 


Where a=fl=l for uniform flow, and a >ff>l for non uniform velocity profiles. 

Notes : For ideal incompressible fluid we prove that -.the quantity ( z+p\y ) is 
constant over flow cross section normal to the streamlines when they are 
straight and parallel, which is often called the hydrostatic pressure distribution as 
it is the same relation for a fluid at rest. For the established flow of a real fluid in 
prismatic passage (where the shear stresses perpendicular to the direction of 
motion may be neglected, which is the case of almost engineering problems except for 
very low Reynolds number and low velocity), the hydrostatic pressure distribution is 
not affected by viscosity, but only by the streamline curvature. 

Therefore the quantity any flow cross section, and the Energy line concept is no 
more valid, and must be reexamined. 

Taking A and C as typical streamlines, Fig. 4. 5, it is noted that each streamline is 
associated with a different energy line, i.e. the flow along different streamlines 
possesses different, amount of total energy. 

However, for flow-in passage, as pipe, duct, or open channel, properties of 
individual streamlines are seldom of interested the whole flow is characterized by a 

single effective energy line a distance above the hydraulic gradient line. 

2g 

Another consequence ( effect ) of non-uniform velocity distribution can be shown 
by considering a flow through a short constriction, Fig. 3.6. From continuity V 2 >V1 
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Hydrostatic relation is valid 


P + z 


= const > 


f P ' 

— + z 


Ji 


\Y 


= const 


J 2 


Frictional effect of constriction is relatively small, and as the acceleration is an 
efficient process, the velocity profile at sec is flatter than that at sec 1 , then 
accordingly d 2 (& 1 (correction factor) 






— (EL) C 

— (EL), 
HGL 



v Horiz. datum 

And now let us calculate exact and conventional head loss, exact hi 12 and hL 12 
( head loss is given by drop in energy line) 

Exac, + 7 + + V + 


And hu-2 is given by ignoring a, 


h., = 





) 



comparing both equations gives 
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h, 


Exact h r 


+ («2 — !) 




Oi — 1) 


Yl 

2g n 


as V 2 >Vj and Oi 2 ( Ot 1 , some compensating features exist, and therefore Exact hu 2 ~hus 


• Only one value of change of (+ z), between section 1&2 

The work-energy equation for incompressible fluid motion in pipes is: 

Z 1 + P± + a 1 ^- = Z 2+ ^ + a 2 ^ + h L (3) 

y 2g y 2g 

However, in most problems of pipe flow, a term may be 
omitted for several reasons: 


(1) Most engineering pipe flow problems involve turbulent flow in 
which a is only slightly more than unity. 

(2) In laminar flow where a is large, velocity heads are usually 
negligible when compared to the other terms. 

(3) The velocity heads in most pipe flows are usually so small 
compared to other terms that inclusion of a has little effect on the final 
result. 


(4) Engineering answers are not usually required to an accuracy 
which would justify the inclusion of a in the equation. 

Example 2 


Assuming Fig. 7.24 to represent a parabolic velocity profile in a passage bounded by two 
infinite planes of spacing 2 R and maximum velocity v c , calculate q, a, and /3. 

Solution: 


The necessary equations are Eqs. 4.8 and 7.51, namely, 

q = VA = V(2R) = J J v d A 

A 

ff v’dA ff 


dA 


V 2 


//- 


and 


dA 


B = 

P V 


II 


v dA 


Taking r as the distance from centerline of the passage to any local velocity, v, and element 
of area dA, dA - dr. The equation of the parabola is v = u c (l - r 2 /R 2 ). 
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Since cj also equals 2/? V, V = 2i; c / 3 . 

~ g) S <» 

(^) 2 I < 2 *^> 

^ _ lA —t 1 z zz 

(c¥) I 

The meaning of these figures is that the exact velocity head is more than 54% greater than 
V 2 /2 g n , and the exact momentum flux 20% greater than QpV. Differences of this mag- 
nitude warn the engineer that a and should be considered when applying the energy and 
momentum equations to one-dimensional laminar flow problems — unless their effects can 
be shown to have negligible consequence in the results desired. 


35 


1 .54 


- 5 * 1 20 


4.1.4- NA VIER-STOKES EQUATIONS 


We would like to derive the differential form of the momentum conservation 
law for incompressible flows . We shall first choose an infinitesimally small 
control volume AV = Ax-Ay-Az around the point (x, y, and z) and then 
shrink this volume to a point by allowing AV Recall from the Newton’s 


Z F b +S F s =ma 

second law of motion: 


Where, the left hand side represents the sum of the external forces and the 

right hand side represents the inertia forces. Since a = 2X represent the total 

Dt 


acceleration at a point, 


Let us review some facts about stresses. The stress tensor that you were 
taught in the strength of materials and represented by the matrix below has 9 
components. 


Note that “cj” represents a normal stress and “x” represents a shear stress. 
Since stress is defined as a force over an area, each stress component has 

F 

two subscripts. For example, g xx = — represents a force in the x-direction 

A x 

that is also acting over an area that is in the x-direction. 

Now, an area may be considered positive or negative, depending on whether 
the outward normal to the area points in the positive or negative coordinate 
direction. For example, in the diagram below, we have chosen a positive and 
a negative y-area. 
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Thus, the stress due to the force F y shown on the positive y-face will 
represent a positive stress. Similarly, a positive stress on the negative y-face 
must be due to a force pointing in the negative y-direction, P y . We may be 
able to represent shear stresses the same way. Question: Can you identify the 
shear stress due to the force Q x shown in the above diagram? What will its 
subscript be and will it yield a positive or negative shear stress? 

In the next set of figures we assume that the stress tensor is known at the 
center of a small control volume AV = Ax • Ay • Az which is centered at (x, y, 
z). To eliminate the stress values on the surfaces of the control volume we 
apply the Taylor series expansions as before. We assume the stresses are all 
positive (i.e., forces act on positive faces in positive direction, or on negative 
faces in negative direction.). One such typical arrangement is shown on two 
y-faces below. 
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Let us examine the term 


da yy Ay^ 

5- — — 

yy 8y 2 

V : J 


AxAz .The origin of 


c?yy(x,y,z) is at the center of the volume. However, since we would like to 

obtain the value of the stress on the “-y ’’-face, we must expand 
Ay 

Qyy (x, y- — ,z) using Taylor series about the cfyy (x,y,z) and its y- 
derivatives at (x,y,z): 


Ay 


da 


cJx,y-—,z,t) = (Tjx,y,z,t) 

2 oy 


yy 


Ay 1 d 2 a 


Ay, 


2 2 oy 2 


But this provides only the shear value, the corresponding force in the “-y 
' SCyy Ay' ' 

yy a, 2 


direction is given by 


AxAz as a first order approximation. 


Similarily the other stress terms are written. 
Derivation ofNavier's Equations 


Now let us focus on only the forces that act in the y -direction ( instead of 
those that act on y -faces). Once this result is derived we may be able to write 
similar expressions for the x- and z-directions. For this purpose, use the 
following sketch: 
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AF = 

Sy 


yy 


dcr a A 

XV Ay 

dy 2 


f 


Ax A z 


dcr a A 

,yv Ay 


yy 


r 


xy 


Zx A A 
xy Ax 

Zx 2 


f 


AyAz — 


xy 


V 


dr 


T + - 

dz 


zy A z 

2 


AxAy - 


dr 


^ & 


z j Az 

2 


dy 2 dx 2 


; dy 

2 

dx 

xy 

A 

Ax 

Zx 

A 

' 2 ^ 
J 

■ Ax Ay 

J 


dz 

2- zy 

A. t 


AxAz 


Ay Az 


dz 2 


ZF S y +ZF B y =ma y 

/Vote that; the above shown are the only possible surface forces in the y- 
direction on the control volume. Therefore: 


f 


dx 


xy 


da 


+ 


dx ^ 
yy _j_ zy 


v Zx Zy Zz y 


AxAyAz 


Body Force along y => AFgy = pBy • AxAyAz 

Now we are ready to write the y-component of the momentum equation: 

Since the mass of the small control volume chosen is pAV = pAxAyAz , and 


Dv dv dv dv dv , 
a v = — = — + u — + v — + w — , thus, 
y Dt dt dx dy dz 


d't 


xy 


Zcr , dr. 


\ 


dx 


+ ■ 


yy 


zy 


Zy dz 


AxAyAz + pB AxAyAz 


. . . ( dv dv dv dv 

pAxAyAz \-u b v \-w — 

^ dt dx dy dz 

Cancelling AV from both sides we get: 


r C)X Xy d<Jyy C'T Zy A 


dx 


- + - 


- + - 


dy dz 


+ pB v =p 


dv dv dv dv 

bU bv b w — 

dt dx dy dz 
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The first three terms on the left hand side represent surface forces per unit 
volume and the fourth one is the body force per unit volume. Similiarly, the 
right hand side terms collectively represent the inertia force per unit volume. 


3X xy 

Let us examine the term — once again. This term is the surface force on 

dx 

the control volume that arises due to the net shear force between two “x” 
surfaces but acting along the y-direction. Likewise the second and third terms 
are due to the normal force and shear force between two y -surfaces and two 
z-surfaces respectively. 


Note there is a pattern which emerges from the summation of all the y -forces 
on the control volume. The stress terms' second subscripts are all y, 
representing the direction of the force, but the first subscripts that represent 
which areas the forces act on, change from x to z for each successive term. 
Similarly the derivatives change also from x to z for each term. The vector 
component is represented by y, and the velocity component is accordingly, v. 
Likewise the x-and the z-equations may also be written rather than derived. 
Finally all three equations are represented by: 


^ da X x 
v <9x 



dr. 


\ 


dz 


PB 5 


( d\x 


du du 





hu h v h w — 

dt dx dy dz J 


v dx dy 


dx \ 

zy 

dz j 


r 



dv dv dy dy ^ 

1- u 1- v h w — 

dt dx dy dz J 




dx. 


dx 


dx 


yz 


dy 


dz 




+ pB z =p 


v dt 


dw dw d\y ^ 
u h V h w 


dx 


dy 


dz 


The above three equations are collectively called the Navier's equations, 
named after their originator. Note that these equations have 4 independent 
variables (x, y, z, and t) but 12 dependent variables (u, v, w, and the stress 
components ). We shall assume that the body force component (which is 
usually due to gravity in mechanical engineering problems ) are known. 
Therefore, the Navier's equations are not solvable since there are 9 more 
unknowns in these equations than the number of equations. 

To make these equations solvable, Stokes proposed a set of constitutive 
relations given below. These relations, together with the continuity equation 
derived before make the Navier ’s equations solvable. The Stoke ’s relations in 
Cartesian Coordinates are: 
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^xy ^yx M' 


yx 


Dv Du 
Dx Dy 


T yz = T zy = P 


zy 


d 5w 05v ^ 
1 

dy <3z y 


^ <3u <3w^ 

X«=T„=^— + — j 


2 - du 

^xx =-P-tMV- V + 2p — 

3 ox 

2 - dv 

=-P--H V V + 2 ,a — 

2 __ — _ 5w 

O-zz = -P - ^ ^ V ’ V + 2p — — 

3 Oz 

In the above relations, /a is the dynamic viscosity and p is the thermodynamic 
pressure. Also note that since the stress tensor is symmetric (i.e., x xy = x yx , etc), 
the first three relations are written in a more compact form rather than using 
6 lines. 


Since our focus in this class is incompressible flows, we may be able to set 
V • V = 0 in the above relations by the use of the continuity equation for 
incompressible flows. If we substitute the above equation into the Navier’s 
equations and repeatedly use the continuity equation for incompressible 
flows (this work is left out for you to try as an exercise ), we obtain the 
simplified form of the Navier ’s equations as: 


Du 

~dt 

dv 

dt 


u 


du 

dx 

dv 


u 


f dw 


X 


dx 

dw 


d t 


u 


dx 


+ v- 


d u 
dy 

dv 

dy 

dw 

1 

dy 


w 


w 


du 

dx 

dv 

dx 


T3 dP 

PB X - ^ + P 

dx 

0 dp 
PBy - — + P 


d 2 u d 2 u d 2 u 


w 


dw 

dx 


dy 

OZ 


dx 2 

d 2 v 


5)2 


dy 2 

d 2 v 


dx 2 

d 2 v 


J 

2.A 


P 


v dx 2 

^ d 2 w 
dx 2 


dy 2 dx 
d 2 w d 2 w 


J 


dy z 


dx 


In the above equations, popularly called the Navier-Stokes equations, we 
have only 4 dependent variables (u, v, w, and p), once again assuming B x , B y , 
and B z to be known. Thus adding the continuity equation to this set: 


- du dv dw 

V = 1 1 

dx dy dz 


0 
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We now have a set of 4 equations in 4 dependent variables, or unknowns. 
Therefore, the fluid dynamic problems for incompressible flows are now 
clearly seen to be solvable. 

Or, simply: 


Note that the first equation (A) is a scalar equation (Since 

V- V = — + — + is a scalar). However, the second equation, (B), hides 

dx dy dz 

3 equations ( one each for the x, y, and z components ) into a single equation. 
Thus we still have 4 equations in 4 unknowns. 

We have already discussed the physical significance of the equation (A). 
Equation (B) is the vector form of the Newton ’s second law of motion (except 
each term has been divided by AV during the derivation ). Thus the left hand 
side represents inertia force per unit volume, while each of the three terms in 
the right hand side represents a type of external force per unit volume (body 
force, pressure force, and viscous force respectivily). 

4.1.5- DARCY- WEISBA CH EQUATION: 

Early experiments (circa 1850) on the flow of water in long, 
straight, cylindrical pipes (Fig.) indicated that head loss varied 
(approximately) directly with velocity head and pipe length, and 
inversely with pipe diameter. Using a dimensionless coefficient of pro- 
portionality f, called the friction factor, which known by 


V ■ V = O 





( B ) 





Fig. 4.6 
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f 


t v 2 
d 2g 


The friction factor f depended primarily on pipe roughness in addition 
to velocity and pipe diameter; more recently it was observed that the 
friction factor f depended on the viscosity of the fluid flowing. From the 
relation between shear stress and head loss derived before; 




Where Ri,= hydraulic radius or hydraulic depth: 

^ A tc c / d R 

h R ~4 Ter / 2 2 

Equations 4.2 and 3.36 may now be combined to give a basic relation 

between frictional stress, T, and friction factor, this is: 

x = 


In this fundamental equation relating wall shear to friction factor, 
density, and mean velocity, it is apparent that, with dimensionless, 
Vt 0 /p must have the dimensions of velocity; this is known as the 

friction velocity, V*, which is related to the friction factor and the mean 
velocity by: 


V* 



(4.4) 


However, the physical meaning of the friction velocity is not 
revealed by this algebraic definition; since it is a velocity which 
embodies only wall shear and fluid density, it is defined by the same 
equation whatever the flow regime ( laminar or turbulent ) or whatever 
the boundary texture (rough or smooth). For this reason it is a useful 
generalization that finds wide application in further developments. 

Example 

Water flows in a 150 mm diameter pipeline at a mean velocity of 4.5 m/s. The head lost 
in 30 m of this pipe is measured experimentally and found to be 5.33 m. Calculate the 
friction velocity in the pipe. 


SOLUTION: 
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We will use Eq. 9.4 to calculate the friction velocity but first we must find the Darcy- 
Weisbach friction factor /. We resort to Eq. 9.2 because all the terms in that equation are 
known except /. 


K = f 


LYl 

d2g n 


(9.2) 


Solving Eq. 9.2 for f gives 


/ = 



2 X 9.81 0.150 m 
(4.5 m/s) 2 30 m 


5.33 m = 0.026 


Now, bringing in Eq. 9.4, 


ft = 4.5 m/s = 0.26 m/s « 


V * = V J s = 45 m/s J T 


(9.4) 


4.2- LAMINAR FLOW 

Analysis of laminar flow in a pipeline (Fig. 4.2) may be begun with 
the following established facts: 

(1) symmetrical distribution of shear stress and velocity. 

(2) maximum velocity at the center of the pipe and no velocity at the 
wall (the no-slip condition). 

(3) linear shear stress distribution in the fluid given by the following 
Eqn (from application of the impulse-momentum principle. 



(4) The shear stress in the laminar fluid flow is proportional to the 
velocity gradient as following: 

dv 

= ^ 
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In Fig. 4.2 it is clear that r / R - y and dr = - dy; thus, equating the 
expressions for r yields 



yh 

At the wall r—R, y—0 and T = T 0 , then i Q = — — R, from which, 

2 / 


^ = lhL, then , * = J_ X= -L 

R 21 dr (j, jlx 

Integrating this relation gives: 
x„ r 2 

V = — 


r yhj^ 

\2l j 


1 T o 

r -r 

JLX R 


pR 


x r 

0 +c 


I iR 2 


x„ R z 


L IV 

From the no slip condition at r=R, V-0, thenO — I -C, then 


\ xR 2 

t R 2 

C = — , from which the velocity profile is found to be parabolic in 

2 pR 

the form: 


2 julR 


(R 2 - r 2 ) 


t„R z 


f 


( 9 . 5 ) 

.2 A 


For r=0 and v — V r — — , from which v = V. 

c 2pR 


R z 


J 


2 T 

Noting that v* = — , so , the previous equation can be written as: 


p 




v = — (R 2 ~ r 2 ) 
2 vR 

or - = (R 2 - r 2 ) 

v* 2 vR 

(9.6) 


or with r 2 — (R — y ) 2 , 

2L _ = ?*( yL\ 

v * V \ 2 R) 


(9-6) 
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V 

— = V 
V* 


V 


y- 


y 


2 \ 


2R 


Thus, the velocity profile can be expressed in terms of distance from the 
wall and is seen to be a characteristic "velocity" of the laminar profile. 

It is interesting to note that near the pipe wall, the term y / 2 R 
becomes negligible compared to y and the velocity profile near the 
wall is essentially linear. In dimensionless form, Eqn. becomes: 


— = — where y(((R 
v* v 


We will approach development of a head loss equation by 
finding the flowrate Q by integrating the velocity profile given by Eqn. 


R 


R 


Q = f v(2ndr ) = — — f (r 2 — r 2 W = — - — , substituting the relation 

J J 4 \ Ll 


for the wall shear stress x = then: 

° 21 


Q = 


TrR^yh L 

8/xZ 


'7T d*yh L 
128 /jlI 


( 9 . 8 ) 


2 

This is Hagen-Poiseuille law. Since Q=nR V, the last equation takes the 
form: 

v = = -y fflflfyL, 

S/mI 


From which head loss h L equal to: 


h L 


32jxlV 
y d 2 


( 9 . 9 ) 


Equation 4.8 shows that in laminar flow the flowrate, Q v which will 
occur in a circular pipe, varies directly with the head loss and with the 
fourth power of the diameter but inversely with the length of pipe and 
viscosity of the fluid flowing. These facts of laminar flow were 
established experimentally, independently, and almost simultaneously 
by Hagen (1839) and Poiseuille (1840). 
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Equation 4.9 shows that, in a laminar flow, head loss varies with 
the first power of the velocity. Equating the Darcy -Weisbach 
equation 4.2 for head loss to Eqn. 4.9 yields an expression for the 
friction factor; it is 


64{x 64 

Vdp~ ' R 


(9.10) 


Thus, in laminar flow the friction factor depends only on the 
Reynolds number. 

Example 


A fluid flows from a large pressurized tank through a 100 mm long, 4 mm diameter tube. 
In a 600 sec time period, 1 300 cm 3 of fluid are collected in a measuring cup. If the head 
loss in the tube is 1 m, calculate the kinematic viscosity v. Check to verify that the flow 
is laminar. 


SOLUTION: 

First, we will calculate the flowrate by the fundamental relationship Q - volume/time. 
The volume is 1 300 cm 3 /(10 6 cm 3 /m 3 ) = 13 X 10 -4 m 3 . The time is 600 s. 


Volume 

Time 


13 X 10~ 4 
600 


= 2.17 X 


10 6 n C/s 


We will use a modified version of Eq. 9.8 to find v. 

= tt d A yh L = rr = 77 d 4 g n h L = 77 d 4 g n h L 

U 128 /j./ 128 /r/ 128 (/i/p)/ 128 v/ 


(9.8) 


Now, solving the above equation for v , 

_ 77 d 4 g n h L _ 3.14 X (4/1 000) 4 X 9.81 X 1 

128(2/ ~~ 128 X (2.17 X 10~ 6 ) X (100/1 000) 


To calculate the Reynolds number, we need the average velocity V. 


Q _ 2.17 X 10~ 6 

A ” (tt/4) X (4/1 000) 2 


0.17 m/s 


Using Eq. 8.11 to obtain Reynolds number, we get 


Vd _ 0.17 X (4/1 000) 
v ~ 0.000 28 


( 8 - 11 ) 


Since R = 2.4 is well below 2 100, we have a strongly laminar flow and the analysis is 
correct. • 
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4.3- TURBULENT FLOW 
4.3.1- SMOOTH PIPES 

Pipe flow with friction effects has been seen to be a viscosity-inertia 
phenomenon and, thus, to be characterized by a Reynolds number. For 
smooth pipes, the discussion about flow past solid boundaries is 
relevant and strongly suggests the existence of a viscous sublayer near 
the pipe walls. 

The total shear stress can be written as: 


Where, p dv/dy is the viscous stress and — pv x v is the turbulent 
(Reynolds) stress. 

It has already been shown that the total stresses are a linear 
function of radius r in a pipe. Experiments show that over most of the 
flow, the turbulent stress dominates, but the maximum stress equals the 
viscous stress at the wall, where the turbulent shear stress is zero. 

The analytical treatment is begun by assuming that the viscous 
stress in Eqn. 9.11 is negligible over most of the flow, employing the 


stress, and equating this to the linear total shear stress relation 


where y is measured from the pipe wall. Now in Section 7.2 it was 
argued that near the wall (y«< R), l = Ky. However, here on the 
centerline (y = R), x = 0, so either 1= 0 or dv/dy = 0 there (or both 
equal zero). Thus, some insight is needed to proceed further; it comes 
from experiment. Nikuradse's systematic and comprehensive 
measurements of velocity profiles, in the form: 



( 9 . 11 ) 



for the turbulent shear 
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ZJ c 17 

— 2.5 in — 


f< 


C^ll pipes) 


[Viscous sheaf stress 



Fig.4. Typical variation of shear stress in a turbulent pipe flow 


dv 1 


y 


From which — oc — , and from equation l oc 1 


dy y 


R 


Then it is follows that, in pipe flow: 


Equation 9.12 now becomes 


i = «y i 



-rjp^y 2 - = (dv/ dy ) 1 


or 


dv/dy — 


Ky 


Ky 


(9.14) 


(9.15) 


Attempts to represent the Reynolds stress have only been partly 
successful; this area represents one of the classic unsolved problems of 
fluid flow. As derived in Section 7.2, Prandtl mixing length theory is a 
plausible and practical, although not theoretically rigorous, approach. 
It is used here to derive the velocity distribution and other quantities 
because the results obtained are close to experimental data and the 
process gives insight to the physics of the flow. 
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Which illustrates again the pervasive presence of the friction velocity 
v*. Integrating the last equation, one has: 

v = — In y -+- C (9.16) 

K 

There are two ways to approach evaluation of this unknown constant 
C. 

- First, ifv = v c at y = R (on the centerline ), then and 

v — v 1 y 

C = v c — (vJK) InR The result is: — = In — 

v* K R 

Second, one can try to find C in terms of the no-slip condition (v = 0) 
air = 0. But at y = 0, v = - oo, according t o this Eqn. This is at least 
unrealistic (and not unexpected because the neglected viscous stresses 
dominate at the boundary)! It must be that the " turbulent " profile is 
replaced by a viscous-dominated profile near the wall and that there is 
some appropriate distance from the wall which marks the onset of this 
process- distance y' (nominal sublayer thickness 8 V , ( 8 V = y' ). 

Now very near the wall in a laminar flow, the velocity profile is linear, 
and if we match a linear laminar flow velocity profile with the turbulent 
profile at some distance y' from the boundary, and if we include 
additional experimental data from Nikuradse work, we come to the 
general equation of the velocity profile for turbulent flow in smooth 
pipes in the form: 

v luv 

— = 2.5 In 4- 5.5 (9.17) 

u* v 


Or — — 5. 75 log — ^ + 5.5 

V* V 


v — v 1 y 

Equating the equation — = In — , and the equation for laminar 




K R 


v v* y 

shear stress near the wall , we have: — = 

v* v 
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in smooth pipes (5 X 105 < R < 3 X 106) and in pipes with a 
uniform sand grain roughness showed that all velocity profiles ( not just 
in smooth pipes) could be characterized by the single equation 

Therefore, it must be that dv/dy l/y and, from Eqn. 9.12, (y/R). Near the 
wall = Ky it follows that, in a pipe flow, (9.14) Equation 9.12 now 
becomes dy/dy 

We can find the intersection of both where velocity is the same at 
which y=y'=S v , nominal sublayer thickness which gives: 

^ = 11.6 or 5 V = 11.6 — 
u v* 


These results are illustrated, together with some experimental data in 
the following figure. 



The flow rate in turbulent pipe flow is given: 

Q = J vClirr dr) = 2-7 ru„ J j^5.75 log + 8.5 J r dr 

yielding 

V = Q/rrR 2 = L>* |^5.75 log j + 4.75 J 


V v*R 

— = 5.75 log h 7. 75 for smooth pipe and 

v* u 

From which, 

— = 1 + 4.07 
V 


also for smooth pipes 


[FLUID MECHANICS (1-B)] 


Page 26 


INTERNAL FLOW 


Which is adjusted to conform with experiment. 

Fig. 4.5 gives a comparison (for the same flow rate and mean 
velocity) of typical velocity profiles for laminar and turbulent pipe 
flow. 

Next it is possible to developed a relation between f and Re, as results 




(Smooth pipes) 


( 9 . 21 ) 


Where Re = Vd/v is the pipe Reynolds number. 


Also we can reach the following relation for laminar sublayer: 


8 V _ 11. 6u 

d v*d 


11.6v 



3Z8 

Re 47 


, which gives Re 



3Z8 
(5, /d) 


Then substituting in the right hand side of relation between f and Re we 


have: 


1 

47 


2.0log[32.8/(bjdj\-0.8 


Which gives, that for turbulent flow over smooth wall the friction factor 
is a function only of the ratio of the sublayer thickness to the pipe 
diameter. 
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PROBLEM 9.4 


Water at 20°C flows in a 75-mm diameter smooth pipeline. According to a wall shear 
meter (Section 14.11), t q — 3.68 N/m 2 . Calculate the thickness of the viscous sublayer, 
the friction factor, the mean velocity and flowrate, the centerline velocity, the shear stress 
and velocity 25 mm from the pipe centerline, and the head lost in 1 000 m of this pipeline. 

Solution: 

Because the friction velocity appears in so many of the relevant equations, we will calculate 
that quantity first. By definition, 


Water at 20°C flows in a 75-mm diameter smooth pipeline. According to a wall shear 
meter (Section 14.11), r Q — 3.68 N/m 2 . Calculate the thickness of the viscous sublayer, 
the friction factor, the mean velocity and flowrate, the centerline velocity, the shear stress 
and velocity 25 mm from the pipe centerline, and the head lost in 1 000 m of this pipeline. 

solution 

Because the friction velocity appears in so many of the relevant equations, we will calculate 
that quantity first. By definition, 


Equation 9.19 will be used to compute the average velocity although Eq. 9.4 would do as 
well. 



Now we employ Eq. 9.18 to compute sublayer thickness. 

8 V = 11.6(v/u # ) = 11.6(1 X 10“ 6 m 2 /s)/(0.061 m/s) (9.18) 

= 1.9 X 10 -4 m = 0.19 mm • 

PROBLEM 9.4 



Now we employ Eq. 9.18 to compute sublayer thickness. 

8 V = 11.6(v/u*) = 11.6(1 X 10“ 6 m 2 /s)/(0.061 m/s) 

= 1.9 X 10 -4 m = 0.19 mm • 

We turn to Eq. 9.23 to calculate the friction factor. 


(9.18) 


'—p = 2.0 log[32.8/(S /ri)] - 0.8 
Vf 

= 2.0 log[32.8/(0.19 mm/75 mm)] - 0.8 = 7.42 
f = 0.018 ® 


(9.23) 


V 

— = 5.75 log 
v* 


v*R 


v 


+ 1.75 


(9.19) 



+ 1.75 = 21.1 
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V = 1.29 m/s • 

Q = VA = 1.29 m/s X (tt/4) X (0.075 m) 2 = 0.005 7 m 3 /s • 
To compute the centerline velocity, we utilize Eq. 9.20. 


— = 1 + 4.07 V/78 = 1 + 4.07 V.01 8/8 = 1.193 (9.20) 


v c = 1.54 m/s ® 

We know that shear stress varies linearly with radius. The point in question is 25 mm 
from the pipe centerline or 2/3 of the way to the wall. As a consequence, the shear stress 
at this point is 2/3 of the value at the wall or 2.45 N/m 2 • 

Equation 9.17 can be used to find the velocity at this point (remember, y is measured 
outward from the wall). 


v v+y 

— = 5.75 log + 5.5 = 5.75 log 
v * v 


0.061 m/s X .012 5 m 
1 X 10" 6 m 2 /s 


+ 5.5 


22.1 (9.17) 


v = 1.35 m/s ® 

Finally, we calculate the head loss in 1 000 m of pipe using the Darcy- Wei sbach formula. 


hL , f il j _ Q018 1 000 m (1.29 m/sf 
L d2g„ 0.075 m 2 X 9.81 


20.4 m © 


(9.2) 


Before the development of the foregoing generalizations by Prandtl, von Karman, and 
Nikuradse, a pioneering effort was made by Blasius to relate velocity profile, wall shear, 
and friction factor for turbulent flow in smooth pipes. Although Blasius’ s work has been 
superseded by these generalizations, it is still of some importance in engineering (in spite 
of its limited scope and empiricism) because of a mathematical simplicity which allows 
easy visualization and leads directly to useful (but approximate) results. 

First, Blasius 3 showed that the curve representing the friction factor (for 3 000 < R 
< 100 000) could be closely approximated by the equation 


0.316 

r0.25 


(Blasius) 


(9.24) 


When this is substituted in the Darcy- Weisbach Eq. 9.2, it is noted that h L « y 1 - 75 for the 
turbulent flow in smooth pipes with R < 10 5 . 

Substituting Eq. 9.24 into Eq. 9.3 produces 

T " = aRvl)^ = 0 033 2 ^ AR ~ , 74 V 7 / V /4 ( 9 . 25 ) 


Blasius then assumed that the turbulent velocity profile could be approximated (see 
Fig. 9.6) by a power relationship 



For this equation the mean velocity V may be related to the center velocity v c by applying 
Eq. 4.8: 


V'ttR 2 



2ttCR - y)(~dy) 
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which gives 

V _ 2 

v c (m + 1 )(m -f 2) 


(9.26) 


from which, by substitution into the profile equation, we derive 

2 


V = 


(m + 1 )(m + 2) 


f R 

v I - 


Substituting this expression for V into Eq. 9.25, 

1 7 /4 


t 0 = 0.033 2 


(m + 1 )(m + 2) 


- 1/4 + (Jm/4)y7/4y - (7m/4)^3/4 


However (Blasius reasoned), wall shear r 0 could depend only on the form of the velocity 
profile and the physical properties of the fluid but could not be affected by pipe size R\ 
thus the exponent of R must be zero, and from this the exponent m must be equal to 1/7. 
Validation of this depended on experimental measurements of turbulent velocity profiles 



Fig. 9.6 

which were found to agree quite well with the hypothesis above; thus the so-called seventh- 
root law for turbulent velocity distribution has been widely accepted. It is written 



(9.27) 


A useful corollary of this law is an equation for wall shear r 0 in terms of v c and R. 

Using m = 1/7 in Eq. 9.26 yields V/v c = 49/60, and, substituting 49u c /60 into Eq. 9.25 
for V and rearranging, 


t 0 = 0.046 4 


JL 

■PcPR; 


1/4 


pv] 


(9.28) 


which will be used later in the approximate analysis of turbulent boundary layers. 
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PROBLEM 9.5 

For the conditions of Illustrative Problem 9.4, calculate the friction factor, wall shear stress, 
centerline velocity, and the velocity 25 mm from the pipe centerline using the seventh- 
root law. 

Solution: 

The friction factor can be calculated from Eq. 9.24. However, we will first calculate the 
Reynolds number R. 



v 


1.29 m/s X 0,075 m 
1 X 10 -6 m 2 /s 


= 96 750 


Checking to see if we are in the range of applicability of the Blasius approximation, we 
note that 96 750 is greater than 3 000 and less than 100 000; so, the Blasius approach 
should yield good results. 

Now, calculating the friction factor. 


0.316 

R 0.25 


0.316 

(96 750) 0,25 


= 0.018 ® 


(9.24) 


Before calculating wall shear stress, we must find the centerline velocitv We co to 
Eq. 9.26. S 





1.29 m/s - 1.58 m/s ® 


(9.26) 


which were found to agree quite well with the hypothesis above; thus the so-called seventh- 
root law for turbulent velocity distribution has been widely accepted. It is written 



(9.27) 


A useful corollary of this law is an equation for wall shear t 0 in terms of v c and R. 
Using m = 1/7 in Eq. 9.26 yields V/v c = 49/60, and, substituting 49u c /60 into Eq. 9.25 
for V and rearranging. 


t„ - 0.046 4 


— 

v cP Rj 


1/4 


pu2 

2 


(9.28) 


which will be used later in the approximate analysis of turbulent boundary layers. 
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4.3.2- ROUGH PIPES 


Pipe friction in rough pipes at high Reynolds numbers will be 
governed primarily by the size and pattern of the roughness, since 
disruption of the viscous sublayer will render any laminar flow action 
in that region ineffective. However, experiments show, as noted above, 

v — v 1 y 

that the logarithmic velocity profile given by— = — — /ft , is 

applicable for both smooth and rough pipes. 

A development similar to that used for smooth pipe flow, except 
with the average roughness height e used instead of the viscosity, and 
tempered by the experimental work of Nikuradse on pipes with sand 
grain roughness, yields the following equation for the velocity profile. 


u y 

= 5.75 log - H- 8.5 


(Rough pipes) 


(9.29) 


As before the flowrate Q 

is found by 

integration: 

R 

f 

R r 

rr 

( r -A 


Q- \v( 2nrdr ) = 2nv * 

5.75 log 

+ 8.5 

J 

0 

{L 

l e J 

- 


V = -Qj = v* 

nR 2 


5.75 log - + 4.75 

e 


rdr , yielding to 


V R 

, thus — — 5.75 log — I -4.75 one 
v* e 


v if 

can used this equation to verify that = 1 + 4.07 J— , which is 
clearly, then valid for all turbulent flows cases. 

Using v* = V \Jf/8 , in the previous equation this gives: 


1 Re 

-i= = 2.03 log — + 1.68 , 

V/ e 


this equation is adjusted by help of 


Nikuradse experimental results, and the final equation is: 


-f= = 2.0 log — + 1.75 
A If e 


Or 


— = 2.0log — + 1.14 

V/ e 


For turbulent flow in rough pipes the friction factor is a function 
only of the relative roughness e/R or e/d and is not a function of 
Reynolds number (compare Eqn. 4.31 to Eqn. 4.21). Friction effects 
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are produced in fully rough flow by roughness alone, without 
dependence on viscous action. 

PROBLEM "REPORT" 

The mean velocity in a 3 00 -mm pipeline is 3 m/s. The relative 
roughness of the pipe is 0.002 and the kinematic viscosity of the water 
is 9 X 10~ m /s. Determine the friction factor, the centerline velocity, 
the velocity 50 mm from the pipe wall, and the head lost in 300 m of 
this pipe under the assumption that the pipe is rough. 


4.4- FRICTION FACTOR 

4.4.1- CIRCULER PIPE 

The evaluation of pipe friction factors is extensively based on 
experimental work so it is informative to approach the representation 
of the friction factor from the standpoint of dimensional analysis and 
similitude. From the analyses presented in previous sections, it is 

known that the wall shear stress T 0 depends on the mean velocity V, the 
pipe diameter d9 the mean height of the wall roughness projections e, 
the fluid density p, and the fluid viscosity /x. Thus 

<f> ( To , V", d, e, p, pty — O 


Using the Il-theorem (Section 8.2), with V, d, and p as the repeating 
variables, gives 



n = 


Vdp 

/x 


n 3 


€3 

d 


Setting j as a function of U2 and n3, we obtain 

- - •*** iff- 1 ) - ^ («■ §) - (*■ 5 ) 

Comparison with equation 4.3 shows that 

( 9 . 38 ) 

Thus, the dimensional analysis approach confirms the experimental 
results and previous analyses that the friction factor depends only on 
the Reynolds number of the flow and the pipe relative roughness e/d. 
The physical significance of Eqn. 4.38 may be stated briefly as the 
friction factors of a number of pipes will be the same if their Reynolds 


f = <f>‘ 


!' |^K, 
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numbers, roughness patterns, and relative roughness are the same. 
When this is interpreted by the laws of similitude, its basic meaning is 
the friction factors of the pipes are the same if their flow pictures in 
every detail are geometrically and dynamically similar. 

Eqn 4.38 suggests a convenient means of presenting the experimental 
data on the friction factor. This was first used by Blasius in 1913 and 
by Stanton6 in 1914 and consists of a logarithmic plot of friction factor 
versus Reynolds number, with relative roughness as the parameter. 
Figure 4.9 shows Nikuradse's data plotted in the Blasius -Stanton 
format, revealing the relationship between, f , R, and e/d. Recall that 
Nikuradse's experiments employed fixing sand grains on the inside of 
the pipe to establish an easily measurable index of roughness. The 
purpose was to demonstrate conclusively the inter-relationship among, 
R, and e/d. in that regard, the effort was successful and the following 
important fundamentals were clearly established. 



Fig. 9.9 Biasius-Stanton diagram with Nikuradse's experimental data. 

1. The physical difference between the laminar and turbulent flow 
regimes is indicated by the change in the relationship offto R near the 
critical Reynolds number of 2100. 

2. the laminar regime is characterized by a single curve, given by the 
equation = 64/Rfor all surface roughness. This confirms that head loss 
in laminar flow is independent of surface roughness and that hL V. 

3. in turbulent flow a curve of versus R exists for every relative 
roughness, e/d9 and the horizontal aspect of the curves confirms that 
for rough pipes the roughness is more important than the Reynolds 
number in determining the magnitude of the friction factor. 
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4. At high Reynolds numbers the friction factors of rough pipes become 
constant, dependent wholly on the roughness of the pipe, and thus 
independent of the Reynolds number. 

5. Although the lowest curve was obtained from tests on hydraulically 
smooth pipes, many of Nikuradse's rough pipe test results coincide with 
it for 5000 < R < 50000. Here the roughness is submerged in the 
viscous sublayer (Sections 4.3 and 4.5) and can have no effect on 
friction factor and head loss, which depend on viscosity effects alone. 
From Darcy-Weisbach equation h V 17 5 for turbulent flow in smooth 
pipes where the friction factor is given by the Blasius expression ( Eqn. 
4.24). 

6. The series of curves for the rough pipes diverges from the smooth 
pipe curve as the Reynolds number increases. In other words, pipes 
that are smooth at low values of R become rough at high values of R. 
This is explained by the thickness of the viscous sublayer decreasing as 
the Reynolds number increases, thus exposing smaller roughness 
protuberances to the turbulent region and causing the pipe to exhibit 
the properties of a rough pipe. 


MOODY DIAGRAM 

Earlier, Colebrook showed that Nikuradse's results were not 
representative of commercial pipes. Roughness patterns and variations 
in the height of roughness projections in commercial pipes resulted in 
friction factors which were considerably different than Nikuradse 's in 
the transition zone between smooth and wholly rough turbulent flow. 
Benedict gives a particularly lucid discussion of this apparent 
discrepancy. He notes that the work of Colebrook and others on 
commercial pipe rather conclusively demonstrates that the friction 
factor varies as described by Eqn. 4.34 rather than as shown in Fig. 
4.9. Later, Moody presented the equations of Colebrook in graphical 
form using the Blasius -Stanton format. The result of this work is shown 
M Fig. 4.10. This plot, known as the Moody diagram, along with the e- 
values for commercial pipes given in Fig. 4.11, can be used directly in 
the solution of engineering problems. A computer program, DARCY, 
which calculates a friction factor from the Moody diagram, is given in 
Appendix 7. However, a word of caution is necessary. The roughness of 
commercial pipe materials varies widely with the manufacturer, with 
years in service, and with liquid conveyed. Corrosion of pipe wall 
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material and deposition of scale, slime, and such can drastically 
increase the roughness of the pipe and the resulting friction factor. One 
of the practical considerations of pipe flow calculations, which pose a 
difficult problem for the engineer, is the prediction of friction factors 
which will actually be realized after a pipeline is constructed. 
Extensive practical experience is a valuable asset in his situation. 

Example 4.10 

Water at 100°F flows in a 3 inch pipe at a Reynolds number of 80 000. This is a commercial 
pipe with an equivalent sand grain roughness of 0.006 in. What head loss is to be expected 
in 1 000 ft of this pipe? 

Note that this is the same problem as Illustrative Problem 9.8 except we are using 
commercial pipe. 

SOLUTION 


We will use the Moody diagram to find the friction factor but first we must calculate the 
relative roughness ej d. 


e 

d 


0.006 in. 
3.00 in. 


= 0.002 


With R — 80 000 and e/d = 0.002, we go to the Moody diagram of Fig. 9.10 and find 


f = 0.025 5 

From Illustrative Problem 9.8, V = 2.36 ft/s. We can now compute the head loss using 
the Darcy-Weisbach equation. 


h L 


l V 2 

= f ~ — = 0.025 5 
d 2 Sn 


1 000 2.36 2 

3/12 2 X 32.2 


8.8 ft © 


(9.2) 


Note that the commercial pipe with an equivalent sand grain roughness equal to that of a 
pipe lined with real sand grains 0.006 inches in diameter has a 30% higher head loss at 
this Reynolds number. However, both of these pipes would have the same head loss under 
smooth pipe and wholly rough conditions. 


EXAMPLE 


Crude oil at 20° C flows in a riveted steel pipe 1.00 m in diameter at a 
mean velocity of 2.0 m/s. What range of head loss is to be expected in 1 
km of pipeline? 
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Fipft disrnttfcsr - mm 



SOLUTION 


The question posed above suggests there may be some uncertainty in the evaluation of the 
head loss. We will begin our investigation by calculating the Reynolds number. From 
Appendix 2, crude oil at 20°C has a density p = 855.6 kg/m 3 and an absolute viscosity 
p. = 71.8 X 1(T 4 Pa-s. 


Vdp _ 2.0 m/s X 1,00 m X 855.6 kg/m 3 
IX ~ 71.8 X ltr 4 Pa-s 


2.4 X 10 5 


Now, to find the relative roughness, we go to Fig. 9.11 for riveted steel pipe 1 m in diameter 
and read 


= 0.000 9 to 0.009 
a 
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So now we know what the uncertainty is about. Apparently there are several ways to 
construct riveted steel pipe with the number, type and spacing of rivets open to question. 
We will take the two extreme values and find the head loss for each. From the Moody 
diagram of Fig. 9.10, 

- = 0.000 9 / = 0.020 5 (By program DARCY, App. 7, / = 0.020 4) 

d 


= 0.009 f = 0.036 5 (By program DARCY, App. 7, f = 0.036 5) 


The Darcy-Weisbach equation 9.2 yields 
l V 2 

h L = f = 0.020 5 


K = f-y. 


d 2 g n 
l V 2 


= 0.036 5 


1 000 2r 

1 2 X 9.81 

1 000 2 2 


= 4.2 m • 


d2g n 1 2X9.81 

The range of head loss is 4.2 m to 7.4 m per km. 


= 7.4 m • 


4.4.1-PIPE FRICTION IN NONCIRCULAR PIPES THE HYDRAULIC 
RADIUS: 

Although the majority of pipes used in engineering practice are of 
circular cross section, occasions arise when calculations must be 
carried out on head loss in rectangular ducts and other conduits of 
noncircular form. The foregoing equations for circular pipes may be 
adapted to these special problems through use of the hydraulic radius 
concept. 

The hydraulic radius, RHI is defined as the area, A, of the flow cross 
section divided by its "wetted perimeter," P. In a circular pipe of 
diameter d, 


R 


h 


d 

4 


or 


d = 4 R h 


( 9 . 39 ) 


This value may be substituted into the Darcy-Weisbach equation for 
head loss and into the expression for the Reynolds number with the 
following results: 


. _ fJ_Yl 

4 R h 2 gn 


( 9 . 40 ) 


R 


V(4 RjPp 


( 9 . 41 ) 


from which the head loss in many conduits of noncircular cross section 
may be calculated with the aid of the Moody diagram of Fig. 4.10. 
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The calculation of lost head in noncircular conduits thus involves the 
calculation of the hydraulic radius of the flow cross section and the use 
of the friction factor obtained for an equivalent circular pipe having a 
diameter d equal to 4R H . In view of the complexities of viscous 
sublayer, turbulence, roughness, shear stress, and so forth, it seems 
surprising at first that a circular pipe equivalent to a noncircular 
conduit may be obtained so easily, and it would, therefore, be expected 
that the method might be subject to certain limitations. 

The method gives satisfactory results when the problem is one of 
turbulent flow but, if it is used for laminar flow, large errors are 
introduced. 

The foregoing facts may be justified analytically by examining further 
the structure of Eqn. 4.40 in which hL varies with l/Rh. From the 
definition of the hydraulic radius, its reciprocal is the wetted perimeter 
per unit of flow cross section and is, therefore, an index of the extent of 
the boundary surface in contact with the flowing fluid. The hydraulic 
radius may be safely used in the equation above when resistance to 
flow and head loss are primarily dependent on the extent of the 
boundary surface, as for turbulent flow in which pipe friction 
phenomena are confined to a thin region adjacent to the boundary 
surface and thus vary with the size of this surface. However, severe 
deviation from a circular flow cross section will prevent the hydraulic 
radius from accounting for changes in head loss, even for cases of 
turbulent flow. Tests on turbulent flow through annular passages, 10 for 
example, show a large increase in friction factor with increase of the 
ratio of core diameter to pipe diameter. 

In laminar flow, friction phenomena result from the action of viscosity 
throughout the whole body of the flow, are independent of roughness, 
and are not primarily associated with the region close to the boundary 
walls. In view of these facts the hydraulic radius technique cannot be 
expected to give reliable conversions from circular to noncircular 
passages in laminar flow. This expectation is borne out both by 
experiment and by analytical solutions of laminar flow in noncircular 
passages. 

EXAMPLE. 11 

Calculate the loss of head and the pressure drop when air at an 
absolute pressure of 101.3 kPa and 1 5° C flows through 600 m of 450 
mm by 300 mm smooth rectangular duct with a mean velocity of 3 m/s. 
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SOLUTION 


We will treat this problem in the conventional fashion substituting 4 R h for d in the for- 
mulas. Since the duct is smooth, we will only need the Reynolds number to find the friction 
factor. From Appendix 2 , we find values of p = 1.225 kg/m 3 and /jl = 1.789 X 10~ 5 
Pa*s for air at the specified pressure and temperature. The hydraulic radius for the duct is 


Rh = 


Area 


0.45 m X 0.30 m 


Wetted perimeter 2 X 0.45 m + 2 X 0.30 m 
Now, we can compute the Reynolds number from Eq. 9.41. 

R = V(4 R h )p = 3 m/s X (4 X 0.090 m) X 1.225 kg/m 3 
M 1.789 X 10” 5 

From the smooth pipe line on the Moody diagram of Fig. 9.10, 


= 0.090 m 


= 73 950 (9.41) 


f = 0.019 


Equation 9.40 can be used to calculate the head loss in 600 m of the pipe. 

600 3 2 


l V 2 

hc = f 4 ~R h 2g~„ = 0 019 4 X 0.090 2 X 9.81 


14.5 m of air 


(9.40) 


The pressure drop is found from the equation 

A p = yh L = pgh L = 1.225 X 9.81 X 14.5 = 174 Pa 


4.4.3- PIPE FRICTION-EMPIRICAL FORMULAS 

The Darcy-Weisbach equation 4.2 has provided a rational basis for the 
analysis and computation of head loss. However, historically, a number 
of empirical formulas have been and still are being used for pipe 
friction calculations in engineering practice. A treatment of pipe 
friction would be incomplete without mentioning at least two of these 
formulas and how they relate to the Darcy-Weisbach equation. 

By far the most widely used of these empirical formulations is the 
Hazen-Williamsl 1 formula. It was developed to permit calculating the 
capacity of pipes to convey water. 

! (XJ.S. Customary units) V — 1 .31 8C ftlv /?2’ 63 S 0-54 i (9.42a) 

I (SI units) V - O.X49 C m „^ 6 ‘.S' 1 ” 1 (9.42b) 

Where, Rh is the hydraulic radius, S is the head loss per unit length 
h L /U and is a roughness coefficient associated with the pipe material. 
Table 1 gives some typical values of the Hazen-Williams coefficient. 
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TABLE I Hazen- Williams Coeffk 


aiad Manning /?-vak 


Extremely smooth pipes — PVC 

Copper, aluminum tubing 

Asbestos cement 

New cast iron 

Welded steel 

Concrete 

Ductile iron (cement lined) 
Vitrified clay pipe 
Riveted steel 
Old cast iron 


Qiw 

n 

150-160 

0.009 

150 

0.010 

140 

0.011 

130 

0.013 

130-140 " 

0.012 

120-140 

0.011-0.014 

140 

0.011 

— 

0.011-0.013 

110 

0.013-0.017 

100 

0.015-0.035 


a These are typical values but, because of variabilities in fabrication, the user should 
consult the pipe manufacturer for recommended values of roughness coefficients. 


Unlike the Darcy -Weisbach equation which can be applied to both 
laminar and turbulent flow over a wide range of fluids and 
temperatures, the Hazen-Williams formula is restricted to turbulent 
flow of water at normal temperature in a limited size range of relatively 
smooth pipes. However, the formula is so constructed that it permits 
direct calculation of the velocity of flow (or discharge) for a known 
allowable head loss. This feature is attractive in the design of water 
pipelines. Conversely, use of the Darcy -Weisbach equation in the same 
circumstances often requires a trial solution as one estimate the 
Reynolds number to find the friction factor and makes successive 
adjustments as one closes in on the answer. 

In view of the formula configuration and the rather indefinite 
descriptions of the roughness coefficients, it is difficult to judge the 
range of validity of the Hazen-Williams formula without wide 
experience in its application. It is not clear whether are measure of 
absolute or relative roughness, whether there is any effect of Reynolds 
number in the formula, whether it applies only in smooth or rough pipe 
situations, etc. These questions can be answered fairly conclusively if 
the formula is rewritten in the form of the Darcy -Weisbach equation. 
Working with the U.S. Customary version, replacing Rh with d/4, S 
with h U and solving for h L gives 




[ 


194 

CL 85 (W )° 15 cT 


-1 


L Yl 


If (V is multiplied by y>) 0, 15, a Reynolds number appears and the 
equation may be written as 
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. I” 194 "1 l V 2 

^ L^ 015 CA^ 5 ^o.°, 5 r o.i 5 J d 

Comparing this equation with the Darcy-Weisbach equation 4.2 and 
taking a nominal v -value of 0.00001 ft /s for water, the 4 'equivalent /- 
value for the Hazen-Williams formula is 

= 1 090 

•t £^a.85 ^O.OlSj^O.15 C9.43) 

The most effective means of visualizing the meaning of this value " is to 
plot the equation on the Moody diagram for several values. This is 
done on Fig. 4.12. The following conclusions can now be drawn in 
regard to the range of applicability of the Hazen-Williams formula 

The Hazen-Williams formula is a "transition" formula working best as 
a smooth pipe formula or in the "early" transition zone between smooth 
and rough pipe flow. Note the = 150 smooth-pipe value applies fairly 
well over a range of Reynolds numbers from 1 00000 on but the plot 
suggests a value higher than 150, say 160, should be used for Reynolds 
numbers above 1000000. Also note that the formula grossly 
underestimates friction factors in smooth pipes for Reynolds numbers 
below 10000 (typical of trickle irrigation systems) and does not apply 
at all for laminar flow. Because the lines for the various values on Fig. 
4.12 all have a constant slope, the Hazen-Williams formula should 
never be used for wholly rough pipes. More specifically, the Hazen- 
Williams formula should be used only 

in zones where the pipe is relatively smooth and in the ' ’early" part of 
its transition to rough flow. If there is a question as to applicability, it 
is always possible to estimate the Reynolds number of the flow and, 
with the proposed value, check Fig. 4.12 to verify the applicability of 
the formula. 

2. Addressing the issue of relative roughness, we see it is clear from the 
above equation for that there is a small relative roughness effect built 
into the equation although the 0.015 exponent dampens the impact 
considerably. However, in the near-smooth'- zone where the Hazen- 
Williams formula applies, the effect of relative roughness is small, 
thereby minimizing the error. 

3. There is definitely a strong Reynolds number effect built into the 
Hazen-Williams formula which renders it accurate in the "near- 
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smooth” transition range for Reynolds numbers greater than 100000 or 
with the appropriate value The second formula selected for 
examination is the Manning equation, noted for its application to open 
channel flow but also used for pipe flow. The Manning formulas for 
pipe flow are given as 


(9.44a) 

(9.44b) 


where Rh and S are as described for the Hazen-Williams formula and n 
is a roughness coefficient. Some typical values are given in Table 1. 

We will investigate the range of applicability of the Manning formula 
in the same manner as was done with the Hazen-Williams formula. 
Again using the U.S. Customary version of the formula, and arranging 
in the Darcy-Weisbach configuration, we obtain. 


(U.S. Customary units) 

V = /t2/3 S l/2 

(SI units) 

V = i R?/ 3 S ,/2 




185n 2 

d ^ 3 


(9.45) 


It seems clear from this equation that there is no Reynolds number 
effect but there is a fairly strong relative roughness effect. This 
equation is plotted in Fig. 4.12 for two different n-values and two 
different diameters to show the effect of relative roughness. The 
following conclusions can be drawn regarding the applicability of the 
Manning formula. 

1. There is no Reynolds number effect so the formula must be used only 
in the wholly rough flow zone where its horizontal slope can accurately 
match Darcy-Weisbach values provided the proper n-value is selected. 

2. The relative roughness effect is correct in the sense that, for a given 
roughness, a larger pipe will have a smaller friction factor. 

3. In a general sense, because the formula is valid only for rough pipes, 
the rougher the pipe, the more likely the Manning formula will apply. 

Any other formula can be compared with the Darcy-Weisbach equation 
in the manner used above. For empirically-based formulas, this 
approach is strongly recommended unless the user is confident the 
formula is applicable to the case at hand. 
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4.5- LOCAL LOSSES IN PIPELINES 
Separation: 

Separation has great impact on design and performance of 
certain familiar internal flow systems. For example, separation of the 
flow from the blades of a pump impeller 'or turbine rotor (runner) can 
cause. Careful design can prevent separation in these cases, but for 
flow through a valve or metering orifice or' -sudden expansion, 
separation is inevitable. Fig shows ideal and real fluids for the three 
cases. 




Ideal fluid 








Real fluid 


The problem of separation is a major problem in a diffuser design (a classic 
internal flow). There is an optimum diffuser angle a - 6-8°, Fig. 32. 



For greater angle, adverse pressure gradient increases which cause 
flow separation accompanied by great losses, while for small a 
great friction surface exist and shear stress is high accompanied 
also with great losses. 

The "Axiom" mentioned before is also valid for internal flows, 
that is "Acceleration of real fluid tends to be efficient process, 
deceleration an' inefficient one". Accelerating motion through a 
convergent nozzle is accompanied by a favorable pressure gradient 
which stabilized the boundary layer and thus minimizes energy 
dissipation. Decelerating motion is accompanied by an adverse 
pressure gradient which "promotes separation, instability, eddy 
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formation and large energy dissipation. This axiom may be extended 
to explain the difference between the maximum efficiencies obtained 
in comparable complex machines such as hydraulic turbines and 
centrifugal pumps. In turbines the flow passage are predominately 
convergent and the flow is accelerated; in the pumps the situation is 
opposite. Hydraulic turbine efficiencies have been obtained up to 94%, 
but maximum centrifugal pump efficiencies are around 87%. The striking 
difference between these figures can be attributed primarily to the 
inherent efficiency and inefficiency of the acceleration and deceleration 
flow process, respectively. 

4.6- SECONDARY FLOW 

Let us consider the flow through a pipe bend, Fig. 4.33. 



The fluid particles velocities at the center are the same, due to 
which arise centrifugal force which is in equilibrium with pressure force. 
The reduction of velocity at the outer part A, of the bend reduces the 
centrifugal force acting on particles moving near the wall. 

As result of this, non-equilibrium between pressure and 
centrifugal forces, the twin eddy motion shown are developed, 
that is the secondary flow, when added to the main flow will cause 
a double spiral motion. The combination of double spiral flow and 
separation leads to increased head losses in bends, a large part of 
which may be caused by increased wall shear stress in and 
downstream from the bend because of the redistribution of the flow 
streamlines by the secondary flow. 

Into the category of local losses in pipelines fall those losses 
incurred by change of cross section, bends, elbows, valves, and 
fittings of all types. Although in long pipelines these are distinctly 
minor losses and can often be neglected without serious error, in 
shorter pipelines an accurate knowledge of their effects must be known 
for correct engineering calculations. 
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The general aspects of local losses in pipelines may be obtained from a 
study of the flow phenomena about an abrupt obstruction placed in a 
pipeline (Fig. 4.13), which creates flow conditions typical of those 
which dissipate energy and cause local losses. Local losses usually 
result from rather abrupt changes (in magnitude or direction ) of 
velocity; in general, increase of velocity ( acceleration ) is associated 
with small head loss but decrease of velocity ( deceleration ) causes 
large head loss because of the production of large-scale turbulence. In 
Fig. 4.13, useful energy is extracted to create eddies as the fluid 
decelerates between sections 2 and 3, and this energy is dissipated in 
heat as eddies decay between sections 3 and 4. Local losses in pipe 
flow are, therefore, accomplished in the pipe downstream from the 
source of the eddies, and the pipe friction processes in this length of 
pipe are complicated by the superposition of large-scale turbulence on 
the normal turbulence pattern. 
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Separate action of the normal turbulence and large-scale turbulence, 
although in reality a complex combination of the two processes exists. 
Assuming the processes independent allows calculation of the losses 
due to established pipe friction, h L2 and h u and also permits the loss 
hL, due to the obstruction alone, to be assumed concentrated at section 
2. This is a great convenience for engineering calculations since the 
total lost head in a pipeline may be obtained by simple addition of 
established pipe friction and local losses without detailed consideration 
of the above-mentioned complications. 


Early experiments with water (at high Reynolds number ) indicated that 
local losses vary approximately with the square of velocity and led to 
the proposal of the basic equation 



( 9 . 46 ) 


in which K L , the loss coefficient, is, for a given flow geometry, 
practically constant at high Reynolds number; the loss coefficient tends 
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to increase with increasing roughness and decreasing Reynolds 
number, 12 but these variations are usually of minor importance in 
turbulent flow. The magnitude of the loss coefficient is determined 
primarily by the flow geometry, that is, by the shape of the obstruction 
or pipe fitting. 


When an abrupt enlargement of section (Fig. 4.14) occurs in a pipeline, 
a rapid de-celebration takes place, accompanied by characteristic 
large-scale turbulence, which may persist in the larger pipe for a 
distance of 50 diameters or more downstream before the normal 
turbulence pattern of established flow is restored. Simultaneous 
application of the continuity, Bernoulli, and momentum principles to 
this problem has shown that (with certain simplifying assumptions), 




( 9 . 47 ) 


2Note that this is the same trend followed by the friction factor, in 
which K l = 1. 



Experimental determinations of K L confirm this value within a few 
percent, making it quite adequate for engineering use. A special case of 
an abrupt enlargement exists when a (relatively small) pipe discharges 
into a (relatively large) tank or reservoir. Here the velocity 
downstream from the enlargement may be taken to be zero, and when 
the lost head (called the exit loss) is calculated from Eqn. 4.47 it is 
found to be the velocity head in the pipe. 

The loss of head due to gradual enlargement is, of course, dependent 
on the shape of the enlargement. Tests have been carried out by 
Gibsonl3 on the losses in conical enlargements, and the results are 
expressed by Eqn. 4. 47, in which K L is primarily dependent on the cone 
angle but is also a function of the area ratio, as shown in Fig. 4.15. 
Because of the large surface of the conical enlargement which contacts 
the fluid, the coefficient loss coefficients for conical enlargements. 
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embodies the effects of wall friction as well as those of large-scale 
turbulence. In an enlargement of small central angle, K L will result 
almost wholly from surface friction; but, as the angle increases and the 
enlargement becomes more abrupt, not only is the surface area 
reduced but also separation occurs, producing large eddies, and here 
the energy dissipated in the eddies determines the magnitude of K L . 
From the plot it may be observed that: 

(1) There is an optimum cone angle of about 7° where the combination 
of the effects of surface friction and eddying turbulence is a minimum. 

(2) It is better to use a sudden enlargement than one of cone angle 
around 60°, since KL is smaller for the former. 



Gradual enlargements in passages ( termed diffusers) of various forms 
are widely used in engineering practice for pressure recovery, that is, 
pressure rise in the direction of flow. No attempt is made here to review 
the extensive literaturel4 on diffusers, but it should be realized that 
Gibson's results cannot give a reliable solution to this problem. His 
tests were made, as are all such tests for local losses, with long straight 
lengths of pipe upstream and downstream from the enlargement ( see 
Fig. 4.13). The designer, however, is frequently more interested in the 
pressure rise through the diffuser and substitution of a short nozzle for 
the upstream pipe length. The tests of Gibson and others have shown 
that the pressure will continue to rise for a few pipe diameters 
downstream from section 2, owing primarily to readjustment of velocity 
distribution from a rather pointed one caused by deceleration through 
the diffuser to the flatter one of turbulent flow. 15 From this it may be 
concluded that pressure rise through the diffuser, computed from 
application of the Bernoulli equation with data from Fig. 4.15, will be 
larger than that which will actually be realized. Substitution of a nozzle 
for the upstream pipe length will alter the inlet velocity distribution 
from the standard one of turbulent flow to a practically uniform one 
with a thin boundary layer. The effect is to reduce the losses by 
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stabilizing the flow and delaying separation; not only do smaller loss 
coefficients result, but the cone angle for minimum losses is larger, 
thus allowing a shorter diffuser for the same area ratio. 

EXAMPLE. 13 


A 300 mm horizontal water line enlarges to a 600 mm line through a 
20° conical enlargement. When 0.30 m /s flow through this line, the 
pressure in the smaller pipe is 140 kPa. Calculate the pressure in the 
larger pipe, neglecting pipe friction. 

SOLUTION 


In determining the head loss in a conical enlargement, Eq. 9.47 applies with the K L -value 
obtained from Fig. 9.15. First, we will calculate the velocities in each pipe. 


Q _ 0.30 m 3 /s 

A 300 ~~ (V4X0.300 m) 2 


4.24 m/s 


Q 0.30 m 3 /s 

A 600 (77-/4X0.600 m) 2 


1.06 m/s 


From Fig. 9.15, K L = 0.43. To compute the pressure in the large pipe, we turn to the 
work-energy equation 7.46 without pumps or turbines. 


^300 


7*300 

y 



Z 600 


7*600 

y 


^600 


+ h L 


(7.53) 


where h L = K L ^22 Ymff 

2 Sn 

Taking the datum as the pipe centerline eliminates z from the calculations leaving 

140 X 10 3 Pa (4.24 m/s) 2 _ p 600 (1.06 m/s) 2 (4.24 - 1.06) 2 

9 800 N/m 3 2 X 9.81 y 2 X 9.81 + ' 2 X 9.81 


P 600 

y 


14.6 m 


7*6oo = 14.6 X 9 800 = 143 000 Pa = 143 kPa • 

This is the pressure to be expected a metre or so downstream from the end of the enlarge- 
ment. 


Flow through an abrupt contraction is shown in Fig. 4.16 and is 
featured by the formation of a vena contracta and subsequent 
deceleration and expansion of the live stream of flowing fluid. 

Experimental measurements of K L are somewhat conflicting in 
magnitude although they exhibit a well-established trend from 0.5 for 
A 2 /Ai = 0 to 0 for A 2 /Ax = 1. In view of this it is entirely adequate for 
engineering practice to use a synthesis of analytical approaches and 
generally accepted experimental information. 16 The result is given in 
Table 2 where Cc = Ac/A2. 
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A square-edged pipe entrance (Fig. 4.17a) from a large body of fluid is 
the limiting case of the abrupt contraction, with A2/AI = 0 (A is 
virtually infinite here). The head loss is expressed by 

Eqn. 4.46 in which KL is close to 0.5 for highly turbulent flow, as 
mentioned above. 



TABLE 2 Contraction, & and Loss, A A, Coefficients for Abrupt 
Contractions 


TABU 

!2 Contrs 

iction, C c , 

, and Loss, K l: Coeffic 

ients for Al 

>rupt Coni 

Tactions 




A 2 Mi 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 1.0 


0.617 

0.624 

0.632 

0.643 

0.659 

0.681 

0.712 

0.755 

0.813 

0.892 1.00 

h 

0.50 

0.46 

0.41 

0.36 

0.30 

0.24 

0.18 

0.12 

0.06 

0.02 0 


The entrance of Fig. 4.11b is known as a re-entrant one. If the pipe 
wall is very thin and if the plane of the opening is more than one pipe 
diameter upstream from the reservoir wall, the loss coefficient will be 
close to 0.8, this high value resulting mainly from the small vena 
contracta and consequent large deceleration loss. For thick-walled 
pipes the vena contracta can be expected to be larger and the loss 
coefficient less than 0.8; Harris 17 has shown that pipes of wall 
thicknesses greater than 0.05J (if square-edged) will give a loss 
coefficient equal to that of the square-edged entrance. 

If the edges of a. pipe entrance are rounded to produce a streamlined 
bell-mouth (Fig. 4.18) the loss coefficient can be materially reduced. 
Hamiltonl8 has shown that any radius of rounding greater than 0.\4d 
will prevent the formation of a vena contracta and thus eliminate the 
head loss due to flow deceleration. The nominal value ofK L for such an 
entrance is about 0.1, but its exact magnitude will depend on the 
detailed geometry of the entrance and structure of the boundary layer. 
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The head loss caused by short well-streamlined gradual 
contractions (Fig. 4.19) is so small that it may usually be neglected in 
engineering problems. However, an appreciable fall of the hydraulic 
grade line over such contractions is to be expected; in the pipe 
downstream from the contraction the hydraulic grade line will be found 
to slope more steeply 



K l = 0.8 (b) Re-entrant 4.17 Pipe entrances. 


than that for established flow because of the change of velocity 
distribution and boundary layer growth, which cause an increase of a 
These effects are known to be small and may be ignored in many 
problems. A nominal value of K L (for use in Eqn. 4.46 with V — V 2 )for 
short well- streamlined contractions is 0.04; by careful design this 
figure may be lowered to 0.02, but for long contractions values much 
larger than 0.04 are to be expected because of extensive wall friction. 

Losses of head in smooth pipe bends are caused by the combined 
effects of separation, wall friction, and the twin-eddy secondary flow 
described in Section 7.14; for bends of large radius of curvature, the 
last two effects will predominate, whereas, for small radius of 
curvature, separation and the secondary flow will be the more 
significant. The loss of head in a bend is expressed by Eqn. 4.46 in 
which the head loss is the drop of the extended energy lines (see Fig. 
4.13) between the entrance and exit of the bend. Reliable and extensive 
information on this subject will be found in the work of ltd, 19 a small 
but typical portion of which is shown in Fig. 4.20. Head loss 
coefficients for smooth pipe bends provide another example of the 
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dependence of K L on shape of passage ( determined by $ and RJd) and 
Reynolds number; the research of Hofmann20 provides information on 
the ( expected ) dependence of loss coefficient on relative roughness as 
well. To the engineer the most significant feature of head loss in bends 
is the minimum value of K L occurring at certain values of RJd which 
allows selection of bend shapes for maximum efficiency in pipeline 
design. 



Fig. 9.19 Gradual contraction. 



Tests made on bends with R/d = 0 have given values of K L around 1.1. 
Such bends are known as miter bends (Fig. 4.21) and are used widely 
in large ducts such as wind and water tunnels, where space does not 
permit a bend of large radius. In these bends, installation of guide 
vanes materially reduces the head loss and at the same time breaks up 
the spiral motion and improves the velocity distribution downstream. 

The losses of head caused by commercial pipe fittings occur because of 
their rough and irregular shapes which produce excessive large-scale 
turbulence. The shapes of commercial pipe fittings are determined 
more by structural properties ease in handling, and production 
methods than by head loss considerations, and it is, therefore, not 
feasible or economically justifiable to build pipe fittings having 
completely streamlined interiors in order to minimize head loss. The 
loss of head in commercial pipe fittings is usually expressed by Eqn. 
4.46 with V the mean velocity in the pipe and K L a constant (at high 
Reynolds numbers), the magnitude of which depends on the shape of 
the fitting. Values ofK L for various common fittings are available in the 
Engineering Data Book of the Hydraulic Institute; typical values are 
presented in Table 3. 
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It is generally recognized that when fittings are placed in close 
proximity the total head loss caused by them is less than their 
numerical sum obtained by the foregoing methods. Systematic tests 
have not been made on this subject because a simple numerical sum of 
losses gives a result in excess of the actual losses, and thus, produces 
an error on the conservative side when design calculations of pressures 
and flowrate are to be made. 




k l - 0.2 


Fig. 9.21 Miter bends. 


4.7- PIPELINE PROBLEMS— SINGLE PIPES 

All steady-flow pipe problems may be solved by application of the 
work-energy and continuity equations, and the most effective method of 
doing this is the construction of energy and hydraulic grade lines. 
From such lines the variations of pressure, velocity, and unit energy 
can be clearly seen for the whole problem; thus, the construction of 
these lines becomes equivalent to writing numerous equations, but the 
lines lend a clarity to the solution of the problem which equations alone 
never can. 

In engineering offices, tables, charts, monograms, computer software, 
and so forth, are employed where numerous pipe-flow problems are to 
be solved. Although all these methods are different, they have their 
foundations in the work-energy principle, usually with certain 
approximations; no attempt is made here to cover these many 
methods — the following discussion will be confined to the application 
of the work-energy and continuity principles and the use of certain 
approximations. 

Engineering pipe-flow problems usually consist of (1) calculation of 
head loss and pressure variation from flowrate and pipeline 
characteristics, (2) calculation of flowrate from pipeline characteristics 
and the head which produces flow, and (3) calculation of required pipe 
diameter to pass a given flowrate between two regions of known 
pressure difference. The first of these problems can be solved directly, 
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but solution by trial is required for the other two. 22 Trial-and-error 
solutions are necessitated by the fact that the friction factor, /, and loss 
coefficients, KL, depend on the Reynolds number, which in turn 
depends on flowrate and pipe diameter, the unknowns of problem types 
2 and 3, respectively. However, many engineering pipeline problems 
involve flow in rough pipes at high Reynolds numbers. Here trial 
solutions are seldom required (1) because of the tendency off and K L 
toward constancy in this region, (2) because of the inevitable error23 
in selecting f from Fig. 9.10, and (3) because engineering answers are 
usually not needed to a precision which warrants trial-and-error 
solution in the light of the foregoing facts. Construction of energy and 
hydraulic grade lines for some typical pipeline problems will indicate 
further approximations which may frequently be used in the solution of 
engineering problems. 

Consider the calculation of flowrate in a pipeline laid between two 
tanks or reservoirs having a difference of surface elevation H (Fig. 
4.22). The energy line must start in one reservoir surface and end in 
the other; using a gradual drop to represent head loss due to 




Fig. 9.23 

Unless special plots are devised for circumventing this. Due to the 
inexactness of definition of the roughness. 
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pipe friction, h L f, and abrupt drops to represent entrance and exit 
losses, h Lg and h L the energy line is constructed as shown. It is 
apparent from the energy line that 

hr “I - Hj Hr ~ FT 

U e L^J L. x 

which is the work-energy equation written between the reservoir 
surfaces. When the appropriate expressions for the head losses are 
substituted. 


(° 5 + + l ) Tt " " 

If turbulent flow is assumed and a nominal value for f of 0.03 selected 
(Fig. 4.10), the quantity in parentheses becomes 4.5, 31.5, and 301.5 
for -values of 100, 1 000, and 10 000, respectively. The quantities 0.5 
and 1.0, which result from inclusion of local losses, have a decreasing 
effect on the solution with increasing if these terms were omitted 
entirely, errors of about 18, 2, and 0.3%, respectively, would be 
produced in the velocity and flowrate. Evidently, the effect of local 
losses in pipelines of common length is so small that they may often be 
neglected entirely, appreciably simplifying calculations. Another 
convenient approximation accompanies the above; increasing l/d also 
decreases V2/2gn and thus brings the energy line and the hydraulic 
grade line closer together; since V /2g is of the order of the local 
losses, it is consistent to neglect this also, thus making energy and 
hydraulic grade lines coincident ( except near the entrance ) and 
necessitating the construction of only one line. When the single line is 
drawn for this pipeline problem (Fig. 4.23), the equation may be 
written, and the velocity and flowrate may be obtained by trial-and- 
error procedure. 


a *. 


f 


L 


H 


The foregoing approximations are convenient in engineering 
problems but, of course, cannot be applied blindly and without some 
experience; preliminary calculations similar to those above will usually 
indicate the effect of such approximations on the accuracy of the result. 

EXAMPLE.! 5 


A clean cast iron pipeline 0.30 m in diameter and 300 m long connects two reservoirs 
having surface elevations 60 m and 75 m. Calculate the flowrate through this line, assuming 
water at 10°C and a square-edged entrance. 
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SOLUTION 

Before using the work-energy equation to calculate the flowrate, we need to make some 
preliminary calculations. The Reynolds number depends on the unknown velocity but we 
will introduce all the known quantities into the expression for R. From Appendix 2, for 
water at 20°C, v = 1.306 X 10“ 6 m 2 /s. 



v 


V X 0.30 m 
1.306 X 10~ 6 m 2 /s 


= 229 000 V 


From Fig. 9.1 1, e/d = 0.000 83 for clean cast iron pipe. To arrive at a reasonable estimate 
of the Reynolds number so that a good first estimate of the /-value can be found, we 
assume V ^ 2 m/s. Inserting this number into the expression for R, 


R = 229 000 X 2 = 458 000 


From Fig. 9.10 with e/d = 0.000 83 and R = 458 000, we find / — 0.020. Now, we use 
this value along with the local loss coefficients in the work-energy equation. Those local 
loss coefficients are: 


K l = 0.5 for a square-edged entrance 
K L = 1.0 for exit into a reservoir 


We are now prepared to utilize the work-energy Eq. 7.35. 


, + Ex . Yl _ _ . £2 

Z\ + H- 0 — <-2 + 

y 2 s„ y 


vi 

2g„ 


(7.35) 


Note 

that the choice of 

points 1 and 2 at 

the surface of the two 

reservoirs 

P\ = 

/? 2 = 0 and V v = 

= O. 





75 + O + 

O = 60 -h O + O 

-+- ^0.5 

+ / i + 1-0 

a 

) — 
J 2g„ 


15 = ^0.5 

0.020 X 300 m 

4- 1.0^ 

V 2 

1.096 V 2 


0.30 m 

2 X 9.81 


V = 3.70 m/s 





This is considerably above our estimate so we must re-compute R, find a new /"-value and 
find V again through the work-energy equation. With V = 3.70 m/s, R = 847 250. From 
Fig. 9.10, f = 0.019 3. Substituting this value into the work-energy equation yields 

V = 3.76 m/s 

This is clearly close enough to the first calculated valtie that we can be confident another 
iteration is unnecessary. The flowrate is 

Q — VA = 3.76 m/s X (-77-/4X0.30 m) 2 = 0.266m 3 /s «* 

EXAMPLE.! 6 

A smooth PVC pipeline 200 ft long is to carry a flowrate of 0.1 ft 3 /s between two water 
tanks whose difference in surface elevation is 5 ft. If a square-edged entrance and water 
at 50°F are assumed, what diameter of pipe is required. 
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SOLUTION 

This type of problem represents the most difficult to solve because the unknown diameter 
renders both the Reynolds number and the relative roughness as unknown. The only sim- 
plification in the case of smooth pipes is that the relative roughness is irrelevant, i.e., we 
are on the smooth pipe line on the Moody diagram. 

We again turn to Appendix 2 to obtain the v - value of 1.41 X 10“ 5 ft 2 /s for water at 
50°F. We can also use the technique of representing the Reynolds number in terms of the 
flowrate as 

_ Vd _ Qd _ Q _ 0.1 ft 3 /s _ 9 020 

v Av ('7r/4) d X v (-77-/4) d X 1.41 X 10“ 5 ft 2 /s d 

Following the previous illustrative problem and applying the work-energy equation, 
we get 


5 


0.5 + / 


200 


1.0 


XL 

Ign 


We are now involved in a trial-and-error solution for V. The best procedure is to set up a 
table for the solution process. 


Assume d R f V V 2 /2g n ( ) Right-Hand Side 


0.25 

36 000 

0.022 

2.04 

0.064 4 

19.1 

1.23 

0.20 

45 100 

0.021 2 

3.18 

0.157 

22.7 

3.56 

0.18 

50 100 

0.020 8 

3.93 

0.240 

24.6 

5.90 

0.187 

48 200 

0.021 0 

3.64 

0.206 

24.0 

4.94 


The value of d = 0.187 ft produces a reasonable match to the left-hand side of the equation 
(**1% error). Hence, we will have a diameter of d — 2.24 inches, which means we will 
specify the next larger nominal diameter d = 2.5 inches. © 


Another interesting application of the work-energy principle occurs when a pipeline 
extending from a reservoir terminates in a nozzle. This situation for cases where the main 
line velocity head is both significant and negligible is shown in Fig. 9.24. However, in 
either case, the velocity head at the nozzle exit cannot be neglected. In the most general 
case where we will consider local losses, the work-energy equation is 


7 2 g n y 


% + „ L 
2 S, 1 
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1 ) Crude oil at 20 °c flows in riveted steel pipe 1.00 m in diameter at a mean 
velocity of 2.0 m/s. What range of head loss is to be expected in 1 km of 
pipelines? For Crude oil at 20°c: v=8.391 *10' 6 m 2 /s 

2 ) Calculate the loss of head and the pressure drop when air at an absolute 
pressure of 101.3 kPa and 15 °c (p= 1.225 kg/m 3 ) flows through 600 m of 450 
mm by 300 mm smooth rectangular duct with a mean velocity of 3 m/s. 

For air at 15°c: v=1.46 *10' 5 m 2 /s 

3) A clean cast iron pipeline 0.30 m in diameter and 300 m long connects two 
reservoirs having surface elevation 60 m and 75 m. calculate the flow rate 
through this line, assuming water at 10° c and a square - edged entrance. 

For water at 10°c: v=l. 306 *10~ 6 m 2 /s 

4) A smooth PVC pipeline 200 m long is to carry a flow rate 0.05 m /s between 
two tanks whose difference in surface elevation is 8 m. If a square-edged 
entrance and water at 10 °c are assumed, what diameter of pipe is required? 
For water at 10 °c : v= 1.306 *10~ 6 m 2 /s 

5) A clean cast iron pipeline 0.26 m in diameter and 200 m long connects two 
reservoirs having surface elevation 60 m and 160 m. if given: e/d=0.001, 
v w =1.3*10' 6 m 2 /s, Re>2*10 6 , moody diagram, Ki=0.5 for a square edge 
entrance, Ki -1 for exit into reservoir , 

- Calculate the flow rate through this pipeline. 

- Check the flow type in the above pipeline. 
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